NASA CR-1 32691 


MCR-75-143 

NAS1-13611 

Volume III 

PROGRAMMER'S MANIJAL 


Final 

Report April 1975 


PROGRAM TO OPTIMIZE 
SIM! I L.ATFD TR AJECTOFs ! ES 
(POST) 


(>i AS? -C ; - 1 32 601) P TO G 

ST KUI S T T D T ?■ ?;<7 B CT r ~ ITS 
FROG Ti A >•; ?•' 2 ' ' S M f. T U ? T, I' 
tfarivtta Corp.) 2i-2 


7; ah TO OPTI v ITT 

( n C S T ) . VOLUrli 3: 
inal report ("a rt ir 
' ' CSC I 


1175-32148 


n ticlas 

G.3/1.3 4081b 


Prflj: .a a-:* by 

G. I . B*-; ‘uer. 0 Corn!*, A. R. Rabagor, 


Approved 


ti 9P 

flC ui* 



tv 


D- E . Cornice 
Program Manager 


REPRODUCED BY 

NATIONAL TECHNICAL 
INFORMATION SERVICE 

U S DEPARTMENT OF COMMERCE 
SPRINGFIELD, VA. 2216! 


fcJA.TTiM V'A. COEPCRATiGW 

P’,0, *879 

Denver, Colorado £0201 


1 


m* 



FOREWORD 


This final report describing the formulation of the Program 
to Optimize Simulated Trajectories (POST) is provided in accord- 
ance with Part 3.0 of NASA Contract NAS1-13611. The report is 
presented in three volumes as follows; 

Volume I - POST - Formulation Manual; 

Volume II - POST - Utilization Manual; 

Volume III - POST - Programmers Manual. 

This work was conducted under the direction of Mr. Joseph Rehder 
of the Space Systems Division, National Aeronautics and Space Ad- 
ministration, Langley Research Center* 
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FINAL REPORT 

PROGRAM TO OPTIMIZE SIMULATED TRAJECTORIES (POST) 

VOLUME III - PROGRAMERS MANUAL 

By G. L. Brauer, D. E. Comick, A. R. Habeger, 

F. M. Petersen, and R. Stevenson 
Martin Marietta Corporation 


SUMMARY 


This report documents the program structure and logic, 
subroutine descriptions, and other pertinent programing 
information. 

POST, a generalized point mass, discrete parameter targeting 
and optimization program, provides the capability to target and 
optimize point mass trajectories for a powered or unpowered 
vehicle operating near a rotating oblate planet. POST has been 
used successfully to solve a wide variety of atmospheric flight 
mechanics and orbital transfer problems. The generality of the 
program is evidenced by its N-phase simulation capability, which 
features generalized planet and vehicle models. This flexible 
simulation capability is augmented by an efficient discrete 
parameter optimization capability that Includes equality and 
Inequality constraints. 

POST was originally written in Fortran IV for the CDC 6000 
series computers. However, it is also operational on the IBM 
370 and UNIVAC 1108 computers. 

Other volumes In the final report are: 

Volume I - Formulation Manual - Documents the equations 
and numerical techniques used in POST. 

Volume II - Utilization Manual - Documents Information 
pertinent to users of the program. It describes the Input 
required and output available for each of the trajectory 
and targeting/optimization options. 
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I , INTRODUCTION 


The program was written according to guidelines designed to 
provide complete generality wherever possible without sacrificing 
computational speed or computer storage* The guidelines adhered 
to are: 

1) Computer core size of approximately 104 000 octal; 

2) Fortran IV programming language; 

3) Minimum program execution language; 

4) Modular program construction; 

5) Generalize routines to allow simulation of various types 
of vehicles; 

6) Generality of input, output, targeting, and stopping 
variables ; 

7) Compatibility of operation on both 6500 and 6600 CDC 
computers using either MACE or SCOPE operating systems* 

Information pertinent to the programmer is presented in the 
following sections of this report. Included are descriptions of 
the program logic, flow diagrams, Fortran symbols, and subroutine 
listings . 
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II* PROGRAM STRUCTURE AND LOGIC 


POST is coded exclusively in FORTRAN IV. Overlays are used 
to minimize computer core requirements. The program requires ap- 
proximately 104 000 octal cells of computer storage (see table 
II-l) • More than 104 000 cells may be required to load the pro- 
gram, depending on the operating system used. 

Executive programs are used throughout. These control the 
program flow -by calling subroutines containing the actual mathe- 
matical formulations. This procedure allows the program to be 
modified quickly and easily* 


TABLE II-l* - SUMMARY OF POST OCTAL CORE REQUIREMENTS 


Overlay 

Basic 

program 

Scope 

3.4.1 

system 

Subtotal 

(absolute) 

Blank 

common 

Operating 

total 

(absolute) 

(Q,Q) 

23 141 

3 375 

26 536 

3 720 

32 456 

(1,0) 

31 766 

2 076 

62 623 


66 543 

(2,0) 

24 777 

1 006 

54 544 


60 464 

(2,1) 

5 103 

0 

61 647 


65 567 

(2,2) 

15 205 

74 

72 046 


75 766 

(2,3) 

21 147 

405 

77 674 


103 614 

(2,5) 

6 103 

0 

62 651 


66 576 

(2,6) 

1 676 

0 

56 442 


62 362 

Total required 


77 674 


103 614 
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Overall Program Logic 


POST is structured in three overlay levels, as shown in fig- 
ure II-l. The first overlay (0,0) is the master executive over- 
lay, which controls the overall program. This overlay controls 
the read-in of input data and determines which trajectory compu- 
tations are to be performed. 

Overlay (0,0) first calls overlay (1,0), which reads the 
namelist input data from cards and stores the processed data on 
disc for later use. 

Overlay (2,0) is called by (0,0) after (1,0) has completed 
the input processing tasks. The first decision in overlay (2,0) 
concerns the type of simulation; i-e., single trajectory or search/ 
optimization mode. If a single trajectory is to be run, the pro- 
gram calls overlays (2,1), (2,2), and (2,3) sequentially, then 
returns to the master overlay (0,0). If the search/optimization 
mode is to be used, the program control is turned over to subrou- 
tine MINMYS, which calls overlays (2,1), (2,2), (2,3), (2,5), and 
(2,6) as required to perform the search/optimization function. 

When convergence has been achieved or the maximum number of itera- 
tions has been exceeded, control reverts back to the master over- 
lay (0,0) for the next problem. 

An outline of the approximate calling sequence for each rou- 
tine is presented in the following section of this report. 

This outline shows which subroutines are called by a given 
routine, thereby allowing the detailed logic flow to be followed 
easily. The overall program logic described by the overlays is 
as follows: 

1) Overlay (2,1) reads the previously processed input data 
from tape, locates the data for the current phase (event), 
and initializes the program values based on this input; 

2) Overlay (2,2) initializes the equations of motion for 
the current phase; 

3) Overlay (2,3) integrates the equations of motion from 
time t^ to a specified stopping condition for the cur- 
rent phase; 

4) Overlay (2,5) calculates the control corrections based 
on the search/optimization algorithm being used, limits 
the control parameters that violate the control param- 
eter constraints, and tests for convergence; 
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MINMYS ! 

- Determine optimization 
option 

I 

DINPT OVERLAY (2,1) 

- Reads previously stored 
input data for current 
phase 

i 

PHZXMI OVERLAY (2,2) 

- Initializes equations 
of motion for current 
phase 



Figure II-l. - Program Macrologic 
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CALE 

- Calculates targeting 
errors 


PHZXM OVERLAY (2,3) 

- Integrates the equations 
of motion for the current 
plase from t^ to the 

specified stopping condition 


TEST 

- Tests for convergence 


ITERO OVERLAY (2,6) 

- Prints iteration diag- 
nostics 


GRAD 

- Calculates cost gradient 

- Calculates sensitivity 
matrix 


DELTU OVERLAY (2,5) 

- Calculates control 
corrections 

- Updates control para- 
meters 


Figure II-l.- Concluded 









5) Overlay (2,6) prints out an iteration summary at the end 
of each iteration. It also performs any other informa- 
tion output tasks required by search/optimization algo- 
rithm, such as printing trial step summaries. 

The program dictionary (subroutine DICT) performs a one-to- 
one mapping of variables in common and the Hollerith names by 
which the user can select the variables for a variety of uses, 
including output, stopping conditions, control variables, and 
targeting variables. 

All variables in the dictionary are located In common with 
respect to two labeled commons, IV and END. The first of these 
commons defines the starting reference; the last defines the end- 
ing reference. These commons must bracket all commons required 
by the dictionary. 

POST uses a generalized table storing and look-up procedure 
whereby the size of tables is limited only by the total data 
storage allocation of 1500 cells. Each table has its own multi- 
plier. This is accomplished by dimensioning the table by (2). 

The first location contains the address of the table and the sec- 
ond location contains the table function multiplier. The gener- 
alized table lookup (GENTAB) is set up to handle all allowable 
types of tables, namely, constant-value, monovariant, bivariant, 
and trivariant. 


Outline of Program Logic 

This outline shows the calling sequence for a single itera- 
tion for the trajectory and optimization logic. Certain routines 
are called only if certain options are requested. These routines 
are presented in their approximate calling sequence. The outline 
allows the user to follow the program flow either forward or back- 
ward from a given routine to aid in understanding the logic flow. 
Each subroutine that is called by another routine is listed im- 
mediately below and to the right of the calling routine. 
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1.1 MASTER (OVERLAY! POSTOO, 1, 0 ) ),BLKDAT 

2.1 READAT (OVER LAY ( P0ST10»l»0)) f DICT 

3.1 R SEARCH 

4.1 NMLTER 

3.2 RGENDAT 

1.1 NMLTER 

3.3 RTBLMLT 

4.1 NMLTER 

3.4 RTAE 

4.1 NMLTER 

3.5 INPUTX 

2.2 TSPXM(0VERLAY(P0ST20,2,0) ) ,DATA 

3.1 TRAJX 

4.1 CLSPFL 

3.2 MINMYS 

4.1 NOMINL 

5.1 SETIC 

5.2 TRAJ 

6.1 D I NPT (OVERLAY ( PCST21* 2 » 1 J 

7.1 SETESN 

6.2 SAVIC 

6.3 SETIV 

6.4 PHZXMI (OVERLAY ( POST22 » 2»2 ) 

7.1 CYCXMI 

8.1 DYNXMI 

9.1 DYNXA (OR BYNXB ) 

10.1 DERVI 

11.1 MOTIAL 

12.1 CRBTR 

12.2 WINDS 

12.3 GUIDI 

13.1 QUAT1 

13.2 QUAT2 

12.4 GUIDI 

13.1 PLGM 

13.2 CLGOM 

12.5 IBMTRX 

12.6 WGTINI 

12.7 INTGRL 
12. 6 ATMOS 

10.2 DYSI1 

11.1 DLOOK 

5 . 2 AUXFMI 

9.1 EPHEM 

10.1 FORMN 

10.2 SUN 

8.3 TGOEMI 

8.4 INFXMI 

9.1 PAGER 

9.2 PRNTIC 
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6.5 CALE 

6.6 PHZXM ( OVERLAY (P0ST23t2»3) 

7.1 TGOEM 

7.2 CLGM 


8.1 G5ENSR 

5 . 2 GNAV 

6.3 GGUID 

6.4 GCNTRL 

7.3 CYCXM, GR CYC YM 

7.4 DYNXM 

8.1 DYNS1, DYNS2 » OR DYNS3 

9.1 RUK, SVDQ, OR TWOEDY 

9.2 DERIV 
1G.1 GSA 

11.1 UPNGMS 
11-2 MOTION 

11.3 AUXFM 

11.4 CALSPE 

11.5 CALES 

11.6 GRADS 


10.2 


7.5 AUXFM 

1 8.1 GAMLAM 

8. 2 DGAMLAM 

5.3 6ACKCI 

; £ 6.4 BACKOR 

i - 8.5 XRNGE1 

i t 8.6 XRNGE2 

I 6.7 DPRNG 

L.C CONIC 

8.9 MONITR 

8.10 hSWGT 

8.11 TRACKER 

* 8.12 ANMPT 

6.13 CALSPEC 


MOTION, GR MOTENC 

11.1 ATMOS 

12.1 ATMOS 1, OR 
ATMCS2 , OR 
ATM0S3 

11.2 WINDS 

11.3 GUID1 

12.1 GLGM 

12.2 OLGOM 

11.4 IBMTRX 

11.5 AERO 

11.6 PROP 
12.1 TRIM 

13.1 XITER 

11.7 TMOTM 
12.1 GR AV 

11.3 AEROHI 
11.9 GUID2 

12.1 OLGM 
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7.6 INFXM 

8.1 CGNVG 
b.2 PAGER 

8.2 P6LGCK 

7.7 CLSPFL 

5.3 TEST 

5.4 1TERO (OVERLAY (P0ST26»2»6) 

5.5 SE TIC 

5.6 TRAJ 

4.2 GRAD 

5.1 SETIC 

5.2 TRAJ 

5.3 PAD 

4.3 DELTU(OVERLAY,POST25,2,5) 

5.1 WUCAL 

5.2 GMAG 

5.3 SDM , OR 

5.3 CGM , OR 

5.3 OGMP2, GR 

6.1 DGM 

5.3 PGM 

6.1 UPDATE 

6.2 DROP 

7.1 REVISE 

6.1 UPDATS 

7.2 COMB IN 

7.3 UPDATS 

6.3 DGM 

5.4 UNITDU 

5.5 GABDD 

4.4 TRY II 1 

5.1 GENMIN 

6.1 FGAMA 

7.1 SETIC 

7.2 TRAJ 

7.3 1TERO(OVFRLAY,POST26,2,6) 

6.2 TPOSM 

6.3 THPCSM 

6.4 BUCKET 

6.5 THPM 

6.6 FCPMIN 

5.2 DELTU( OVERLAY, PD ET2 5,2,5) 

5.3 UPiNiCM 

5.4 DELTU(OVERLAY,POST25,2,5J 

4.5 TRY IT2 

5.1 FGAMA 

5.2 GENMIN 

4.6 UPNOM 


(RETURN TO 3.2 AND REPFAT UNTIL CONVERGED) 



Tape or File Designations 

The program uses several Tape (File) designations internally 
to perform the simulation tasks. These files are normally stored 
on discs, but tapes can be used by assigning them the proper file 
designations. The file designations are as follows: 

Tape (or File) Definition 

1 Contains the general data and table multipliers 
for the problem 

2 Contains the initial conditions for each event 
that has a control parameter 


3 and 4 

5 (INPUT) 

6 (OUTPUT) 
8 (PROFIL) 


Store input data for multiple runs 

Stores input data 

Stores output data 

Contains the simulation profile 


Common Designations 

POST uses several labeled commons to provide communication 
between subroutines. In addition, a blank common is used to act 
as a data buffer for the table input data and the event criteria. 
The blank common could be labeled, if desired, without adversely 
affecting the operation of the program. 

The labeled commons are briefly described below in alphabeti- 
cal order. The variables are listed in the following section 
alphabetically to provide an easy cross-reference. 

AUXVC : Common AUXVC contains the variables that are com- 

puted as auxiliaries at the end of each integration step. 

CYCVC : Common CYCVC contains variables and flags used to 

perform cycling functions during forward integration. 

DPGVC; Common DPGVC contains the variables and flags as- 
sociated with the guidance (steering) options. 
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DYNVC : Common DYNVC contains variables and flags required 

to perform dynamics functions during the forward integration. 
Primarily this includes time references and discontinuity flags* 

DYTEM : Common DYTEM contains variables and storage used by 

the integration algorithms to integrate the equations of motion 
forward. No variable in this common. may be input or output. 

END: Common END is used to define the end of the dictionary. 

Any variable defined in a common after common END cannot be input, 
output, or used as a search parameter. This common contains only 
one variable, namely, END. 

GENIC: Common GENIC contains variables of a general nature 

that are required in overlay (0,0) . 

GUIDIC : Common GUIDIC contains the input variables for the 

generalized guidance routines . 

GUIDVC : Common GUIDVC contains the computed variables for 

the generalized guidance routines. 

H0LINC : Common H0LINC contains all of the Hollerith input 

variables. 

INFVCi Common INFVC contains variables and flags that may 
be used in the information output routines at any phase. 

r£: Common IV is used to define a reference to the diction- 

ary region. All variables that are to be input, output, or used 
as search parameters must be defined in a common between common 
IV and common END. IV contains the size of this region. This 
common contains one variable, namely, IV(2) . 

KR0IC: Common KR0IC contains input variables for the vari- 

able step/order predictor-corrector integrator. 

KRgSVC: Common KR0VC contains variables and flags calculated 

in the variable step/order predictor-corrector integrator. 

L0CAL ; Common L0CAL contains parameters used in computing 
the equations of motion and the auxiliary equations that are not 

required to be input or output. If a common variable is to be 
added and it is not needed as an output or an input, it should 
be added to this common. 
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MNMMLT : Common MNMMLT contains a list of mnemonic multi- 

pliers associated with the aerodynamic tables* The first cell 
contains the value 1*0. The remaining cells contain the address 
of any input variable within the dictionary. 

M0TBL : Common M0TBL defines all tables to be interpolated 

by the general table lookup routine GENTAB* Each table requires 
two consecutive storage locations. The first is the table ad- 
dress and the second is the value of the table multiplier. When- 
ever a table is added to this common, subroutines DICT and DATA 
must be modified accordingly. 

M0TIC: Common M0TIC contains all parameters that are re- 

quired as inputs to the equations of motion. Input parameters 
do not have to be defined in this common; however, when such a 
parameter is defined in a lower common (e.g., M0TVC) , the pro- 
gram must search for the dictionary and, hence, runs longer. 

H0TVC : Common M0TVC contains all variables used in the equa- 

tions of motion. These are generally not input or constant par- 
ameters. They are available for output, table arguments, or 
search parameters through the dictionary. 

MU LTRC s Common MULTRC contains the variables associated 
with the multiple-run capability. 

0VRLY25 : Common 0VRLY25 contains the variables required by 

overlay (2,5), which contains the direction-of-search logic. 

PHZVC : Common PHZVC contains flags and constants required 

to perform the phasing fun-tions. 

REDAT: Common REDAT is defined in overlay (1,0) by READAT 

and contains variables and storage data required to build the 
general and table data buffers. 

SEARC : Common SEARC is defined in BLKDAT and, in general, 

contains all parameters required by the iteration algorithms. 
Variables in common SEARC can be input only once per run through 
namelist SEARCH. They cannot be changed through input at a 
phase. Since SEARC is defined in overlay (0,0), it is available 
to every routine in the program. 

SERVC: Common SERVC is a service common available to all 

routines in the program. This common contains 50 cells of tempo- 
rary storage, 5 commonly used index parameters, and a list of the 
most frequently used fixed- and floating-point constants. This 
common should be used whenever possible in order to conserve 
storage. 
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SPECAL; Common SPECAL contains the variables associated with 
the special calculation routine CALSPEC. 

TARGVC : Common TARGVC contains parameters calculated for the 

target vehicle. 

TOOVC; Common TGOVC contains variables and flags required 
to perform the time-to-go functions. 

TRACKC; Common TRACKC contains the variables associated with 
the tracking station option. 



Program Additions 


The guidance, navigation, and flight control routines will 
generally be coded by the user. The coding of these routines may 
require the user to make minor program additions. The most fre- 
quently requested types of program additions are: addition of new 

general variables, addition of generally new tables, and addition 
of new integrals. Instructions for making these additions ' are 
presented in this section. Other types of additions will gener- 
ally require in-depth knowledge of the program code, and a pro- 
grammer familiar with the program should be consulted. 

Addition of new general variables .- General variables are 
any variable that are computed in the simulation portion of the 
program and are to be input, output, used as table arguments, 
search parameters, integrals, or derivatives. The program Execu- 
tive processing algorithm expects to find all general variables 
defined in a labeled common which is loaded between the labeled 
common /IV/ and the labeled common /END/ . The labeled commons 
/IV/ and /END/ are defined in subroutine DICT for overlay (1,0) 
and in subroutine DATA for overlay (2.0). The labeled commons 
defined between /IV/ and /END/ must be in the same order and must 
be the same size in subroutine DICT [overlay (1,0)] and in sub- 
routine DATA [overlay (2,0)]. That is, a one-to-one mapping of 
parameters in subroutine DICT to subroutine DATA relative to 
labeled common /IV/ must be maintained; /END/-/IV/ in overlay 
(1,0) must equal /END/-/IV/ in overlay (2,0). This is absolutely 
required for proper program operation. 

The following steps should be followed to add a general vari- 
able : 

1) Add the new variable(s) to an appropriate labeled common. 
For example, if the new variable is an auxiliary param- 
eter, it should be added to labeled common /AUXVC/ . New 
variables should be added on to the end of an existing 
common. Only the length of the common should change, 

NOT the structure. 

If the user does not want to add to an existing common, 
a new labeled common may be defined, and the new vari- 
ables included in it. However, this is not generally 
necessary. 

2) The labeled common to which the new variables have been 
added, or the new labeled common, is replaced or added 
into subroutine PICT in overlay (1,0) . If a new labeled 
common is being added it must be placed after common 
/IV/ and before common /END/; but not between common 
/MOTBL/ and common /MOTEND/. The locations from common 
/MOTBL/ through common /MOTEND/ are reserved for tables. 
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3) For every new variable added, the Hollerith name by which 
it is to be known must be set into its location for use 
during input processing. This is done by a DATA state- 
ment in subroutine DICT. For example, if a new variable 
called AROANG is added, then the data statement DATA 
AR0ANG/6HAR0ANG/ must be added in the subroutine DICT. 

4) If the new variable Is going to be on input quantity it 
must be added to the NAMELIST/GENDAT/. The input 
NAMELIST/GENDAT/ is defined in subroutine RGENDA in over- 
lay (1,0). Include the new or updated labeled common 
and add the new Input variable to the namelist. 

5) Subroutine DATA, in overlay (2,0), establishes the ini- 
tial or nominal values of the variables to be used in 
the simulation. Every new variable must have a nominal 
value set, even If it is zero. Add the new common, or 
update existing common with new variables, in the sub- 
routine DATA. 

6) Add data statement in subroutine DATA to set nominal 
value of new variable. 

7) Add or change common for new variable In routines where 
it is to be used. Add necessary coding to perform com- 
putations involving new variable. 

Adding new tables ,- The program has a generalized table ac- 
cessing feature that allows new tables to be added without add- 
ing dimensional arrays, hard-wired table arguments, table types, 
table dimensions, etc. The program input processor packs all 
tables input by the user, into an array in blank common. At 
execution time, the table interpolation routine, GENTAB, is 
directed to a particular table in the blank common array by a 
pointer, which is set at data initialization at the beginning of 
each phase. Thus, each table has a pointer associated with it. 
Each table also has a multiplier associated with it, by which the 
table Is scaled during execution. To add a new table the user 
need only add the table pointer and the table multiplier. Be- 
cause the pointer and multiplier can change from phase to phase, 
they are Included In the general data area of program. That is, 
they are defined between labeled common /IV/ and labeled /END/ 
as they are specified in Subroutine DICT and DATA for overlays 
(1,0) and (2,0), respectively. The table input processor expects 
to find all table pointers and multipliers together, and in pairs. 
The pairs must be defined between labeled common /MOTBL/ and 
labeled common /MOTEND/, as declared in Subroutine DICT and Sub- 
routine DATA. A new table pointer and multiplier should be 
added to labeled common /MOTBL/ . A new labeled common could be 
declared between common /MOTBL/ and /MOTEND/ , but this is gen- 
erally not necessary. 


To add a new table the following procedure is to be followed 

1) Add two locations to labeled common /MOTBL/ for the 
table pointer and multiplier. 

2) Replace labeled common /MOTBL/ in Subroutine DICT in 
overlay (1,0). 

3) Add data statement in subroutine DICT to set table name 
in pointer, and Hollerith name of multiplier into table 
multiplier. For example, if a new table called EMFT is 
to be added, then EMFT (2) is added to common /MOTBL/. 

In Subroutine DICT the data statement DATA EMFT/4HEMFT, 
6HEMFTM/ is added. This sets the table name, EMFT, and 
the table multiplier EMFTM for the input processor. 

4) Table multipliers are input through namelist /TBLMLT/* 
Thus the table multiplier must be added to namelist 
/TBLMLT/ in subroutine RTBLML in overlay (1,0). Replace 
the labeled common /MOTBL/. Include equivalence state- 
ment to equate table multiplier with desired input name. 
Add input name to namelist /TBLMLT/. For example, 
EQUIVALENCE (EMFT (2), EMFTM); add EMFTM to namelist 
/TBLMLT/. 

5) Replace labeled common /MOTBL/ In subroutine DATA in 
overlay (2,0). 

6) Add data statement in subroutine DATA to set table 

pointer to zero and table multiplier to desired nominal 
value. Generally the table multiplier will be set to 
1.0. For example, DATA EMFT/0 ,1.0 /. 

7) To reference the new table add or replace labeled common 
/MOTBL/ in routine where table look-up is to be per- 
formed. To perform the interpolation, the interpola- 
tion routine GENTAB is called with the table pointer as 
an argument. For example: 

VOLT = GENTAB (EMFT) 

If the table is not Input, GENTAB will return as zero. 

Adding new integrals .- Any general variable computed in the 
simulation model can be integrated provided it satisfies the neces 
sar^ conditions of differentiable and continuity as required by 
the integration algorithms. The variables must be computed In 
the inner loop of the simulation, and be defined as a general 
variable in Subroutine DICT and Subroutine DATA. The program 
determines which variables are to be Integrated during any phas r 
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from an integration list, which is defined in BLKDAT . The inte- 
gration list contains three entries for each integral. These 
entries are the integral name, the derivative name, and a flag to 
indicate whether this integral is to be integrated or not. Dur- 
ing phase initialization this flag can be set to turn the inte- 
gration on or off. The integration list is defined in labeled 
common /DYNIL/ in BLKDAT. The first location in labeled common 
/DYNIL/ contains the total size of the list including itself. 

Thus, to add an integral the common /DYNIL/ must be increased 
by three, and the contents of DYNIL (1) increased by three. 

To add a new integral the following procedure should be 
followed : 

1) In BLKDAT, overlay (0,0), increase the dimension of 
DYNIL in labeled common /DYNIL/ by three for each new 
integral. 

2) In the associated data statement increase the number 
prestored into DYNIL(l) by three for each integral added. 

3) Add the DATA statment to set the Hollerith name of the 
integral, of the derivities and a nominal value of 0 or 
1, depending when the integral is to be nominally off 
or on, into three new locations defined in DYNIL. 

For example, to add DYNPI as the integral of DYNP , the 
following DATA statement should appear. 

DATA DYNIL/M, 6HTIME, 6HDTIME , 1 


6HDYNPI , 6HDYNP , 0 / 

4) Add integral and derivative, if required, to simulation 
as described previously under addition of general vari- 
ables . 

5) If the user desires to turn the integral on or off as a 
function of input, or model selected, then, the asso- 
ciated flag must be set in the integration list in Sub- 
routine MOTIAL (overlay 2,2). The utility routine INTGRL 
can be used. The first argument is the position of the 
integral in the list, the second is the number of inte- 
grals to be set, and the third the flag zero or one. 


11-16 



III. 


« 


SUBROUTINE DESCRIPTIONS, FLOW CHARTS, AND SELECTED LISTINGS 


This section describes the subroutines used in the program. 
Flow charts and/or listings are also presented in order to show 
the detailed operation of the subroutines. 

Note that the routines are presented alphabetically, rather 
than in the order shown in the previous outline of program logic. 
The outline enables the user to follow the program logic flow from 
one subroutine to another with a minimum of searching to find the 
next routine, but this alphabetical listing makes it easier to 
find subroutines at random. 
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AER0 * This routine calculates the aerodynamic forces and moments 
in the body coordinate system. 
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AER0HI : This routine calculates aeroheating indicators 

that are functions of angle of attack, sideslip, and Mach number. 



I 


APPLY CORRECTIONS 
FOR ALPHA AND BETA 
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AER04: This routine calculates the corrections to the lift 
and drag coefficients (DCLV and DCDV) to be applied to CL and CD 
to account for viscous interaction effects. 



COMPUTE VMU, 
REYN0, AND 
VINV 



I 


ADD CORRECTIONS 
TO CL AND CD 




ANMPT: This routine calculates the vacuum impact point of 

the vehicle and the oblate planet at the specified altitude (ALTIP) . 
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ATM0S1: This routine computes the atmospheric parameters 

using generalized table lookups. 
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ATM0S3: This routine computes the atmospheric parameters 

based on the 1963 Patrick AFB atmosphere model. 



CONVERT ALTITUDE 
TO METERS 


♦ 

IF ALTITUDE ^NEGATIVE, 
SET IT EQUAL TO ZERO 

■ 

DETERMINE THE 
ALTITUDE REGION 
AND EVALUATE 
THE PROPER 
POLYNOMIALS 


l 

COMPUTE ATEM, 
PRES, AND DENS 

i — * 

CONVERT FROM 
METRIC TO ENGLISH 
UNITS IF REQUIRED 




AUXFM: This routine calculates the auxiliary variables at 

the end of each integration step. 
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AUXFMI : This routine initializes the auxiliary calculations. 
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* 


BACK0I: This routine computes the inertial Euler angles, 

given the LB matrix. 
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BACK0R: This routine computes the relative Euler angles, 

given the GB matrix. 








BLKDAT : This routine is not called explicity but performs 

its function of presetting all internal program values at load 
time. Certain of these values can be overridden later by input 
if desired. 



♦DECK, BLKDAT 

BLOCK DATA 
C*** BLKDAT 

C SFT INITIAL PROGRAM VALUES 

C *** THE DATA STATEMENTS IN THIS ROUTINE ARE STANDARDIZED 

C 

C 

C STANDARD FORTRAN DECLARATIONS 

C 
C 
C 
C 
C 
C 

c 
c 
c 
c 

♦CALL MULTRC 
C 

♦CALL SERVC 
♦CALL GENIC 
♦CALL SEARC 

COMMON /INF IV/ INFIV(l) 

COMMON/INFIC/ PF ( 396 ) ,HF ADFR ( 1 0 ) 

COMMON/INFND/ INF ND 
COMMON /DYNIL/ DYNI L { 184 ) 

COMM ON /I NP VC/ INPCF 

♦ IF -IDENT »SMALLC» 1 

COMMON /GENRL/ IGEN12500) 

♦IF -IDENT ,LARC»1 

COMMON IBKTt 2500) 

♦ IF IDENT , LARC ,1 

COMMON IBKT ( 1500) 


COMMON 

DIMENSION 

EQUIVALENCE 

TYPE (CAN BE ANYWHFRE ABOVE) 

DATA (CAN PE ANYWHERE) 

NAMELIST 

- WARNING — 

♦♦♦ COMMON MULTRC MUST BE FIRST COMMON IN BLKDAT ♦*♦ 


C 

c 

DIMENSION DYNILK52) 

DIMENSION DYNIL2( 48 ) 

DIMENSION DYNI L3 ( A8 ) 

DIMENSION DYNILM18) 

DIMENSION DYNIL5 ( 36 ) 

DIMENSION SEARCK782) 

DIMENSION SEARC2(478) 

DIMENSION SE ARC3( 53 ) 

DIMENSION SEARCM53) 

DIMENSION SE ARC5 ( 35 ) 

DIMENSION SEARCM172) 

DIMENSION S£ARC7( 31 ) 

DIMENSION SE ARC8 ( 63 ) 

C 

EQUIVALENCE (DYNIL1 (1 ), DYNIL ( 1)) 

EQUIVALENCE (DYNIL2 ( 1 ) »DYNIL( 53)) 
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c 

c 

c 


c 

c 

c 


EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 


CDYNIL3 Cl ) tDVNILClOl) 1 
(DYNIL4C1 ),DYNIL(I49) ) 
(DYNIL5 ( 1 ) »DYNILC167) ) 
(SEARC1 (1 ),DFDCC1» ) 

(SEARC2Cn,GP(l) ) 

( SEA&C3 ( 1 ) » ITC ( 1 ) ) 

( SEARCH C 1 ) tNAC ) 

C SEARC5 ( 1 ) ,PG1 ( 1 ) ) 

(SEARC6 (1 } ,P2NOM J 

( SEARC7C 1 ) f YES (1 ) ) 

( SE AR C8 ( 1 JfTIMIN ) 


COMMON MULTPC 


DATA 

IN 

/3 

DATA 

10 

/ 3 

DATA 

MULTRF/1 

DATA 

ICASE 

/O 

DATA 

IGSIZ 

/O 

COMMON SEE VC 

DATA 

TEMP 

/50*0 

DATA 

STEMP 

/25+0.0 

DATA 

IR1 

/O 

DATA 

I 

/O 

DATA 

J 

/O 

DATA 

K 

/o 

DATA 

L 

/o 

DATA 

M 

/o 

DATA 

NULL 

/l HU 

DATA 

NOD 

/o 

DATA 

N01 

/I 

DATA 

N02 

/2 

DATA 

N03 

/3 

DATA 

N04 

/4 

DATA 

N05 

/5 

data 

N06 

/6 

DATA 

N07 

/7 

DATA 

N08 

/8 

DATA 

N09 

/9 

DATA 

NIC 

/10 

DATA 

Nil 

/II 

DATA 

N12 

/12 

DATA 

N 13 

/ 13 

DATA 

N14 

/14 

DATA 

N15 

/1 5 

DATA 

FPOO 

/0.0 

DATA 

FPP5 

/ .5 

DATA 

FP1 

/1.0 

DATA 

FP2 

/2.0 

DATA 

FP3 

/3.0 


/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 
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DATA 

FP4 

/4.0 

/ 

DATA 

FP5 

/5.0 

/ 

DATA 

FP6 

/6.0 

/ 

DATA 

FP7 

/7.0 

/ 

DATA 

FF8 

/e.o 

/ 

DATA 

FP9 

/9.0 

/ 

DATA 

F PI 0 

/10.0 

/ 

DATA 

FP11 

/ll.O 

/ 

DATA 

F P I 2 

/12.0 

/ 

DATA 

FP 1 3 

/13.0 

/ 

DATA 

FP14 

/ 14 . 0 

/ 

DATA 

FP1 5 

/1 5 .0 

/ 

DATA 

FP60 

/60.0 

/ 

DATA 

FP90 

/90.0 

/ 

DATA 

FP 1 80 

/180.0 

/ 

DATA 

FP270 

/27G.0 

/ 

DATA 

FP360 

/360.0 

/ 

DATA 

PI02 

/I .5707963267948965/ 

DATA 

PI 

/3. 141 592653 589703/ 

DATA 

P PD 

/O. 01745329251994329/ 

DATA 

DPR 

/57. 2957795 

1308232/ 

DATA 

TWO PI 

/6. 283 185307179586/ 

DATA 

FTPNM 

/6076.1155 

/ 

DATA 

CMPFT 

/ . 3 04 8 

/ 

DATA 

I OF LAG /O/ 



C*** CFOPCE = NEWTONS PER POUND 

DATA CFORCE /4. 4482216152605/ 

C*** CPRES = LB/FT** 2 FER NEWTONS/ME TEP$**2 
DATA CPP.ES /. 02 08 8543 47/ 

C*** CTEMP = DEGREES F PER DEGREES K 
DATA C7EMP /I .8 / 

C*** CDENS = SLUGS /FT** 3 PER KI LOGRAM /METER **3 
DATA CDENS / .00194031965 / 

C*** CHEA T= JOULES PER BTU 


DATA 

CHEAT 

/1054. 350264488888 / 


CMAS 

S = KILOGRAMS PEP : 

SLUG 


DATA 

CMASS 

/14. 5939029 

/ 


DATA 

CVD I ST/6076. 1155 

/ 


DA TA 

IDENT 

/1.0 t 0.0 

,0.0 

f 

1 


0.0 ,1.0 

,0.0 

f 

2 


0.0 ,0.0 

,1.0 

/ 

DATA 

IVSZ 

/O 

/ 


DATA 

XINF 

/l.OE+11 

/ 


DAfA 

INF IV 

/O 

/ 


DATA 

INPCF 

/O 

/ 



C*** 

C 


DATA PE 

1/0 ,6HTIME ,0,6HTIMES ,0, 6HTDUB P ,0,6HnENS ,0,6HPPES » 0, 6HAT EM 
2,0,6HALTIT0,0,6HGCRAD »0»6HGDLAT ,0,6HGCLAT , 0*6HLGNG ,0,6HL0NGI 
3 *0 » 6HVEL I »0,6HGAMMA I , 0, 6HA ZVE LI ,0 ,6HXI ,'0,6HVXI ,0 t 6HAXI 
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oooooooooo 


c 


4.0. 6HVELR 

5 *0 tfcHVTlA 

6 *0 *6HGAMAD 
T,0,6HTHRUST 
8*0* 6HE TA 

9.0. 6HFTXP 
0 ,0 *6HFTYB 
A ,0 *f*HF T2E 
B *0 * 6HC A 
C,0,6HCN 

D ,0 , 6HCY 
F ,0 *6HDYN P 
F,?lf*0 
G/ 


0*6HGAMMAR 
0* 6HGAMMA A 
0»6HAZVAD 
0,6HWEIGHT 
0 »6HETAL 
0 * 6HF AXP 
0 »6HF AYB 
0*6FF AZB 
0*6HCD 
0 *6HC L 
0*6HHEATRT 
0*OHMACH 


DATA HFADFR /10*1H / 


0.6HAZVELR *0,6HYI 
0»6HAZVE LA,0,6HZI 
0,6HDWNRNG,0,fcHCRRNG 
0* 6HWD0T ,0,6HWEIC0N 
0 * 6HI PNULL *0 *6HIYNULL 
0* 6HAXB »0,6HALPHA 
0,6HAYB ,0 ,6HP ETA 

0* 6HAZB ,0,6HBNKANG 
0*6HDFAG ,0,6HP0LI 
0* 6HL I FT *0 ,6HYAWI 
0,6HTLHEAT,0,6HPITI 
0, 6HRE YNO *0*6HASXI 


0,6HVYI *0»6HAYI 
0,6HVZI *0,6HAZ I 
0,6HDPRNG1* 0*6HDPRNG2 
0,6HWPR0P ,0,6HASMG 
0,6HINC PCH, 0, 6H INC YAW 
0 *6HALPD0T *0 » 6HALPT0T 
0,6HPETDOT,0,6HQ ALPHA 
0 »6HBNKD0T *0 *6HQAl TOT 
0,6HYAWR ,0,6HROLBD 
0»6HP ITR ,0,6HPITBD 
0*6HR0L R ,0 » 6HYA WBD 
0 »6HASYI *0,6HASZ I 



1/6HDLR * 6HDLRD ,0 

2.6HTVLR ,6HTVLRD ,0 

3 *6HATL *6HATLD *0 

4*6HGLR *6HGLRD ,0 

5,6HVIDFAL,6HVIDLD ,0 
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6»6HPWPP0P , 6HPWD0T ,0 
7,6HAHI ,6H^HID ,0 
8,6HHTPT »6HHTETD ,0 
9,6HHTTP ,6HHTTPD ,0 
0,6HHTRT »f HHTRTD ,0 
At6HHTLF , 6HHTLFD ,0 

B, 6HHTURP * 6HHTURB D, 0 

C, 6HTIMRF 1.6HDTIMR 1,0 

D, 6HTIMPF?,6HDTIMP 2,0 

E , 6HTIMPF 3 ,6HDTI MR 3,0 

F , 6HTIMRF4 »6HDTIMR 4,0 
G/ 


DATA PYNIL2 


1 

/6HDX 

,6HDVX 

»o 

2 

, 6HDY 

, 6HDVY 

»0 

3 

, 6HCZ 

,6HCVZ 

,0 

4 

, 6HDVX 

,6 HD AX 

»o 

5 

, 6HDVY 

,6HDAY 

»o 

6 

,6HDVZ 

,6HDAZ 

»o 


7 ,6HG INTI , 6HGDEP V 1 , 0 
8 ,6HG INT? ,6WGDFRV2,0 
9»6HG INT3 ,fcHGDEPV3,0 
0,6HGIN'T4 ,6HGDEFV4,C 

A , 6HG INT5 ,6HGDERVE ,0 

B , 6HGINT6 ,6HGDFRV6,0 

C, 6HGINT7 ,6HGDERV7,0 

D, 6HGINT8 ,6HGDERV8,0 

E, 6HGINT9 ,6HGDERV9,0 

F , 6HG INTI 0, 6HGDERV0 , 0 
G / 

DATA DYNIL4 

1 /6HDI AMP1 ,6HDIAPP1 ,0 

2 ,6HDTAMP2,6HDI ARP2,0 

3 , 6HDTAMP3 ,6HD I AR P3 , 0 


L 

, 6HDL I 

,6 HD LID 

t0 

r. 

, 6 W T VL I 

.6HTVLID 

»o 

6 

,6HGLT 

t6HGL ID 

tO 

I 

/ 



DATA PYNIL 5 


1 

/6HDXI 

,6HDVX I 

,0 

2 

,6HDYI 

, 6 HDW 1 

,0 

3 

, fcHDZ I 

♦6HDVZ I 

tO 

4 

, 6HDVX I 

,6HDAX I 

»0 

5 

, 6HDVYI 

,6HDAYI 

tO 

6 

, 6HDVZI 

♦fcHDAZI 

.0 

7 

, 6HDXPT 

,6HDVXPT 

»0 

8 

, 6HDYPT 

,6HDVYPT 

»0 

9 

,6HCZPT 

, 6HDVZP T 

»0 

0 

, 6 H 0 VXPT 

,6HDAXPT 

»o 

A 

, 6HDVYPT 

,6HD AYP T 

,0 

B 

, 6HDVZPT 

,6HDAZPT 

♦ 0 

C 

/ 
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non 


COMMON SEARC 


DATA A COB /625*0 
DATA CONEPS/ B9.9 

1 5*.l 

DATA CONSEX/ .000001 

1 t .001 

DATA CTHA / 

DATA CTHAT /0 
DATA DFPPH / 75*900.0 
DATA DFPTL / 25*1 .0 

DATA DEPVAL/ 25*0.0 
DATA DEPVR / 25*0.0 

DATA SFAPC1 /782*0 
DATA FITEPR/ .C00001 
1 » .001 
DATA GAMA 5 /O 


DATA 

GAMASS/O 


DATA 

GAMAX 

/10/ 


DATA 

SE ARC?/478*0 

DATA 

ICGM 

/63 


DATA 

IDAV 

/O 


DATA 

IDEB 

/ 

0 

DATA 

IDE PVR/ 

25*0 

DATA 

IFDEG 

/ 

25*0 

DATA 

IHADIT/O 


DATA 

IMAX 

/O 


DATA 

IMIN 

/O 


DATA 

INDPH 

/ 

25*0 

DATA 

1N0VR 

/ 

25*0 

DATA 

INTRPL/O 


DATA 

I NTRY1/0 


DATA 

IOPT 

/o 


DATA 

IPRO 

/ 

0 

DATA 

1ST ART/63 


DATA 

searc: 

?/53*0 


DATA 

LIMIT 

/5000B 


data 

MAXITP/ 

10 

DATA 

MODEW 

/ 

1 

DATA 

SFARC4/53*0 


DATA 

NDEPV 

/ 

0 

DATA 

NEQC 

/O 


DATA 

NETF 

/o 


DATA 

NFLAG 

/o 


DATA 

NINDV 

/ 

0 

DATA 

NOMF 

/o 


DATA 

NSTEP 

/o 


DATA 

NTC 

/o 


DATA 

0LDG2 

/o 


DATA 

OLD PI 

/o 


DATA 

0LDP2 

/o 



/ 


» 

/ 


/ 

/ 

/ 

/ 

/ 

/ 


/ 
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DATA 

OIDU /O 



/ 

DATA 

OPT / 

0 

/ 


DATA 

OPTPH /900 

• 

/ 


DATA 

OPTVAP/ 

0 

/ 


DATA 

PCTCC / 

.3 

/ 


DATA 

PCTOLD/ 

.3 

/ 


DATA 

PERT /25*1 . E— A 


/ 

DATA 

PGFPS / 

1.0 

/ 


DATA 

SEARC5/35*C 


/ 

DATA 

P2MIN / 

1.0 

/ 


DATA 

SEARC6/172*0 


/ 

DATA 

SRCHM / 

0 

/ 


DATA 

STEP / 6*0 



/ 

DATA 

STMINP/ 

.1 



1 

» 

.1 

/ 


DATA 

STPMAX/ 1 

. 0 E ♦ 1 0 

/ 


DATA 

TGFAD /O 



/ 

DATA 

TTRIAl/O 



/ 

DATA 

U / 

25*0.0 

/ 


DATA 

OMAG /O 



/ 

DATA 

WCON / 

100.0 

/ 


DATA 

WOPT / 

1 .0 

/ 


DATA 

wu / 

25*1.0 

/ 


DATA 

SEARC7/31*0 


/ 

DATA 

ISFLG /63 



/ 

DATA 

IWTFLG/O 



/ 

DATA 

NPAD /10t 

5,10 / 



DATA 

SEARC8 /83*0/ 




c 

C COMMON GENIC 

C 

DATA IPRT / 63 / 

DATA L IS TIN/2 / 

C 

END 




ORIGINAL PAGE IS 
OE POOR QUALITY 
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CALE : This routine calculates the performance index (PI) 

and the error in the target variables (E(I)). 
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GALES: This rountine calculates the steering errors# 





CALSPEC: This routine is a blank routine to be used when 

special calculations of a temporary nature are required. This 
routine is called at the end of each integration step from AUXFM. 



III-26 



CGM: This routine contains the logic for the conjugate 

gradient method. It is a second-generation unconstrained 
optimization technique that has the stability of the steepest- 
descent method and the convergence properties of the second-order 
techniques . 
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CLGM: This routine contains the executive logic for the 

closed-loop guidance routines. 









CLSPFL (NENDFL) : This routine closes out profile records and 

writes a physical end of file on the profile tape for each trajectory. 
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C0MBIN: This routine determines all combinations of the 

indices of the active constraints. These combinations of con- 
straints are used to determine if any constraints can be dropped, 
and are used when the number of tight constraints exceeds the 
number of independent variables. 
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C0NIC: This routine calculates the Keplerian conic for 

either elliptic or hyperbolic orbits, based on the value of the 
orbital energy. 
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C0NV0:: This routine converts the output variables from 

metric to English, or vice-versa, as required. 



I CO’tfUTE 

1 COWERS I ON FACTORS 



CONVERT m 
AXB (I) , AS 
VIM, WGT, 

1, R, VAM, 
M, VRM, 
AND MASS 


T 

^NPC(l): 


CONVERT C0NIC 
PARAMETERS 


I 


] 


f NPC(5): oV 

j[~ 


1 CONVERT Am 
1 AEROHEATING 

3 SPHERIC AND 
PARAMETERS 



m ■ ■■■■- - 1 

■ w 





f t 

I CONVERT AERODYNAMIC 
1 PARAMETERS 




' 

| CONVERT PROPELLANT 
| WEIGHT PARAMETERS 




^NPC(9)7^^- 


CONVERT THRUST 
AND FLOWRATE 


^NPC(27): Oy 


CONVERT SPECIAL 
FLOWRATE PARAMETERS 


Q NPC(12) :0 ^ 







1 

QjPC( 26 )^ 0 ^- 




CONVERT 

PARAMETI 

HEATING 

HEATING 
:rs FOR 
MODEL 
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CYCXM: This routine performs program-cycling functions, 

which includes checking for any new phase. 



III-33 
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CYCXMI: This routine initializes the program at the beginning 

of each new cycle. 






c 


IS THIS THE FIRST 
INTEGRATION PASS? 


> 


YES 



QUANTIZE THE 
INTEGRATION TIME 


(|) 


SET THE TIME 


SET THE TIME 

REFERENCES TO THE 


TO THE INITIAL 

CURRENT TIME 


TIME 


CALL THE ROUTINES 
THAT INITIALIZE 
THE EQUATIONS OF MOTIONl 


J- — 
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DATA: This routine is not called explicitly, but presets 

the program variables to their stored values at the time overlay 
(2,0) is called. 


i 
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♦DECK, DATA 

SUBROUTINE DATA 
C*** DATA 

C DATA - DEFINES COMPUTATIONAL COMMONS 

C (IV - END) INITIAL DATA VALUES 

C *** THE DATA STATEMENTS IN THIS ROUTINE ARE STANDARDIZED 
C 
C 

COMMON/I V/ I V( 2 ) 

C 

♦CALL AUXVC 
♦CALL CYCVC 
♦CALL DPGVC 
♦CALL DYNVC 
♦CALL INFVC 
♦CALL MNMMLT 
♦CALL HOLINC 
♦CALL MOTBL 
♦CALL MOTIC 
♦CALL MOTVC 
♦CALL PHZVC 
♦CALL SPECAL 
♦CALL TGOVC 
♦CALL TRACKC 
♦CALL GUIDIC 
♦CALL GUIDVC 
♦CALL TARGVC 
C 

COMMON/END/ END 
C 
C 
C 

♦CALL DYTEM 
♦CALL LOCAL 
C 

♦CALL KROIC 
C 

♦CALL KROVC 
C 

COMMON/PRO/ 

1 PRO(502) 

C 

COMMON /GFLAGS/ IGFLAG(IO) 

COMMON /GVARS / GVARS(IOO) 


C 


DIMENSION 

DIMENSION 

DIMENSION 

DIMENSION 

DIMENSION 

DIMENSION 

DIMENSION 


AUXVC 1(82) 
AUXVC2 ( 40 ) 
AUXVC3 ( 23 ) 
AUXVC4C11 ) 
DPGVC 1(33) 
DPGVC2( 9) 
DPGVC3 (77) 
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noon no 


c 


DIMENSION 

DIMENSION 

DIMENSION 

DIMENSION 

DIMENSION 

DIMENSION 

DIMENSION 

DIMENSION 

DIMENSION 

DIMENSION 

DIMENSION 

DIMENSION 

DIMENSION 

DIMENSION 

DIMENSION 

DIMENSION 

DIMENSION 


DPGVC4C 27) 
DYNVC1 { 14 > 
HOLIN1 C 51 ) 
MNMMLK 16) 
MOTE LI (304) 
M0TVC1C1AI) 
M0TVC2C 9) 
MOT VC 3 ( 148 ) 
M0TVC4( 26 ) 
SPECA1 ( 24) 
GUIDV1 ( 35 ) 
TARGV1 (133) 
LOCAL 1C 145) 
KR0IC1C 13) 
KP 0VC1C 27) 
KR0VC2 ( 8 ) 
KR0VC3C 59 ) 


EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 


(AUXVC1C1) 
(AUXVC2 ( 1 ) 
( AUXVC2 ( I ) 
(AUXVCA(l) 
(DPGVC1C1 J 
(DPGVC2C 1 ) 
(DPGVC3C1 ) 
(DPGVC4C 1 ) 
(DYNVCI C 1 ) 
(H0LIN1 ( 1 ) 
(MNMML1C1) 
(M0TBL1C1) 
(MOTVCICI ) 
(M0TVC2C1 ) 
(M0TVC3C 1 ) 
(M0TVC4 ( 1 ) 
CSPECA1 ( 1 ) 
CGUIDV1C1) 
(TARGV1C1) 
(L0CAL1 (1 ) 
(KR0IC1C1) 
(KR0VC1 (1 ) 
(KR0VC2C1 ) 
(KR0VC3C 1 ) 


.ALPTOT 

tCRRNG 

♦URXC1) 

tDPRGI 1 

♦ALPHA 

♦HARG(l) 

♦ YAWI 

♦ EN0MSC1 ) 
♦DTIMRC1) 
♦ALPARG 
♦CADPNM 

♦ CSTC1) 

♦ AH I 
♦LONGI 

♦ DENS 
♦DLID) 
♦SPECI(l) 
♦Gpxim 

♦ALTAT 

♦ AC 1 ) 

♦TINT ) 

♦ DDINTCl) ) 

♦ KQQ) 
♦RNDC) 


DATA IVC2) /O 
DATA PRO /502*0/ 
DATA IGFLAG /10*0/ 
DATA GVARS /100*0/ 


AUXVC 
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no non non 


DATA AUXVC 1/82*0 / 

DATA XMAX /10*-1 .0E10/ 

DATA XMIN /10* 1.0E10/ 

DATA AUXVC2 /40*0/ 

DATA AUXVC3 /23*0/ 

DATA NTRK /I / 

DATA XIVE /0. 0,0. 0,0. 0,0. 0,0. 0,0. 0,0. 0,0. 0,1.0/ 
DATA AUXVC4 /11*0/ 

CYCVC 


DATA 

DELT 

/0.0 


/ 

DATA 

DT 

/1.0 

/ 


DATA 

DTIME 

/1.0 

/ 


DATA 

DTM 

/1.0 


/ 

DATA 

DTO 

/0.0 


/ 

DATA 

ENOIS 

/l.E-8 

/ 


DATA 

TREF 

/0.0 


/ 

DATA 

IDTAB 

/6*0 


/ 

DATA 

IFLG 

/O 


/ 

DATA 

CYCF 

/O 


/ 

DATA 

DELTT 

/O 

/ 



DPGVC 


DATA DPGVCl/33*0 


DATA 

AP 

/l. 0,0.0, 

DATA 

GB 

/l. 0,0.0, 

DATA 

IA 

/l. 0,0.0, 

DATA 

IB 

/l. 0,0.0, 

DATA 

IG 

/l. 0,0.0, 

DATA 

IL 

f\ »0f 0 * 0 t 

DATA 

LB 

/\ m 0 f 0*0 f 1 

DATA 

DPGVC2/9*0 

1/0 

,0 

*0 ,0 

2,0 

,2 

,1 ,0 

3,10*0 


4/ 



DATA 

IVCRT 

/3*0 

DATA 

IVETA 

/O 

DATA 

KDG 

/3*1 .0 

DATA 

KRG 

/3*1 .0 

DATA 

DPGVC3 

/77*0/ 

DATA 

DEPTLS 

/4*1.0/ 

DATA 

MAXITS 

/O/ 

DATA 

NITS 

/O/ 

DATA 

PERTS 

/4*1 • E— 4/ 

DATA 

0PGVC4 

/27*0/ 

DYNVC 



.0*0. 0*1 .0*0. 0*0. 0*0.0* 1 .0/ 
. 0 * 0 . 0 * 1 . 0 * 0 . 0 * 0 . 0 * 0 . 0 * 1 . 0 / 
► 0 * 0 . 0 * 1 . 0 * 0 . 0 * 0 . 0 * 0 . 0 * 1 . 0 / 
► 0 * 0.0 * 1 . 0 * 0 . 0 * 0 . 0 * 0 . 0 * 1 . 0 / 
> 0 * 0 . 0 * 1 . 0 * 0 . 0 * 0 . 0 * 0 . 0 * 1 . 0 / 
> 0 * 0 . 0 * 1 . 0 * 0 . 0 * 0 . 0 * 0 . 0 * 1 . 0 / 
> 0 * 0 .0 * 1 . 0 * 0 . 0 * 0 . 0 * 0 . 0 * 1 . 0 / 
/ 

DATA IGUID 


*1 

*0 


*0 


/ 

/ 


/ 

/ 
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n o o no n non non non 


DATA DYNVC1/14*0.0 / 

DATA DEBINT / 0 / 

DATA DLTMIN /0/ 

DATA DLTMAX /1.0EI0/ 

DATA EPS INT /I.O/ 

HOLINC 

DATA H0LIN1 /51*0/ 

INFVC 

DATA E SNPRT/O / 

DATA EXTRAP, LPRNT/2*0/ 

DATA PINC /0.0 / 

DATA PRNC /I HU / 

DATA FID /IOHUD265 FILE,10H I.D. 0000/ 
DATA INFF»IPRNTP* IPRNTR/3*0/ 

DATA PSTOP /6HPST0P / 

DATA TITLE /10*10H / 

DATA SFID /O ,0 / 

DATA MPINC /O / 

MNMMLT 

DATA ONE /l.O/ 

DATA MNMML1/16*0.0 / 

MOTEL 






DATA 

M0TPL1 






1/0,1 

• t 

0,1., 

0,1., 

0,1., 

0,1 

• Y 

0, 1 . , 

0,1. 


0,1 ., 

0,1 ., 

0,1 . 

2,0,1 

• Y 

0,1., 

0,1., 

0,1., 

0,1 

• Y 

0, 1. , 

0,1. 


0,1., 

0,1 ., 

0,1. 

3,0,1 

• Y 

0,1., 

0,1., 

0,1., 

0,1 

• Y 

0,1., 

0,1. 


0,1., 

0,1., 

0,1. 

4,0,1 

• f 

0,1., 

0,1., 

0,1., 

0,1 

• Y 

0,1., 

0,1. 


0,1s, 

0,1., 

0,1. 

5,0,1 

• Y 

0,1., 

0,1., 

0,1., 

0,1 

* Y 

0,1., 

0,1. 


0,1., 

0,1., 

0,1. 

6,0,1 

• Y 

0,1., 

0,1., 

0,1., 

0,1 

• Y 

0,1., 

0,1. 


0,1., 

0,1 ., 

0,1. 

T»0 » 1 

• Y 

0,1., 

0,1., 

0,1., 

0,1 

• Y 

0,1., 

0,1. 


0,1., 

0,1., 

0,1. 

8,0,1 

* Y 

0,1., 

0,1., 

0,1., 

0,1 

• Y 

0,1., 

0,1. 


0,1., 

0,1., 

0,1. 

9,0,1 

* Y 

0,1., 

0,1., 

0,1., 

0,1 

* Y 

0,1., 

0,1. 


0,1., 

0,1 ., 

0,1. 

0,0,1 

• t 

0,1., 

0,1., 

0,1., 

0,1 

• Y 

0,1., 

0,1. 


0,1,, 

0,1., 

0,1. 

A, 0,1 

- Y 

0,1., 

0,1., 

0,1., 

0,1 

• Y 

0,1., 

0,1. 


0,1., 

0,1., 

0,1. 

B ,0 , 1 

• t 

0,1., 

0,1., 

0,1., 

0,1 

• Y 

0,1., 

0,1. 


0,1., 

0,1., 

0,1. 

C,0,1 

• Y 

0,1., 

0,1., 

0,1., 

0,1 


0, 1 • , 

0,1. 


0,1., 

0,1., 

0,1. 

D,0, 1 

• Y 

0,1., 

0,1., 

0,1., 

0,1 

• Y 

0,1., 

0,1. 


0,1., 

0,1., 

0,1. 

E ,0 , 1 

• Y 

0 Y 1 • Y 

0,1., 

0,1., 

0,1 

• Y 

0,1., 

0,1. 


0,1., 

0,1., 

0,1. 


F ,0 ,1a, 0,1* 

1 / 

MOTIC 
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DATA ALTIP /0-0 / 

DATA ALTREF/IOO. / 

DATA LATREF/1HU / 

DATA LONREF/O.O / 

DATA AZREF /0.0 / 

DATA TIMREF/0.0 / 

DATA ARP /10*0.0 / 

DATA ASH AX /3.0 / 

DATA ATMOSK /l.O ,1.0 / 

DATA AZWB /180. / 

DATA AZL /0.0 / 

DATA LATl /1HU / 

DATA LONL /1HU / 

DATA CLCOMX/O.O / 

DATA DETA /0.0 / 

DATA DESN /3*0.0 / 

DATA DESNE /0.0 / 

DATA DVIMAG,DVMAR,DVPCT/3*0.0/ 

DATA ETAPC /1.0,0.,0. ,0. / 

DATA ETA /l.O/ 

DATA GINT /10*0.0 / 

DATA GO /32.174 / 

DATA GXP »GYP ,GZP/45*0 .0/ 

DATA HEATK /l., 17600. ,26000./ 

DATA ALTITO/O.O / 

DATA HRAT /10*0.0 / 

DATA ISPV /15*1.E11 / 

1960 FISCHFR EARTH MODEL 
DATA J2 /1.0823E-3 / 

DATA J3 /0.0 / 

DATA J4 /0.0 / 

DATA LREF /l.O / 

DATA LREFY /0.0 / 

1960 FISCHER EARTH MODEL 
DATA MU /1.4076539E+16 / 

DATA OMEGA /7.29211E-5 / 

DATA PGCLAT/0.0 / 

DATA PSL /1HU / 

DATA PWPROP/O.O / 

DATA RHOSL /.0023769 / 

DATA RN /l.O / 

1960 FISCHER EARTH MODEL 
DATA RE /20925741 . / 

DATA RP /20855590. / 

DATA SREF /0.0 / 

DATA TSL /1HU / 

DATA WGTSG /l.E-10 / 

DATA WJETT ,WPLD,WPROPI ,WEIC0N/4*0.0/ 
DATA XREF /3*0.0 / 

DATA A EXP /*64 / 

DATA CINF /l.O / 

DATA VINFI /.007 / 



o o o non 


DATA 

IENGA 

/15*1 

/ 

DATA 

IENGT 

/15*1 

/ 

DATA 

ITAP 

/10*0.0 

/ 

DATA 

IWPF 

/15*0.0 

/ 

DATA 

NENG 

/I 

/ 

DATA 

NEQS 

/3*0 *0 

/ 

DATA 

NEWSTG/0.0 

/ 

DATA 

1/0 

,1 

t4 ,2 

»2 

2,0 ,0 

3,15*0 

4/ 

tO 1 0 

tO 

DATA 

GHA 

/0.0 

/ 

DATA 

GHAS 

/180.0 

/ 

DATA 

DECL 

/0.0 

/ 

DATA 

TR PM 

/0.0 

/ 

DATA 

DVMAPR/0.0 

/ 

DATA 

DATE 

/3*1HU 

/ 


tO 

tO 


tl 

»0 


tO 

,1 


,0 

,0 


MOT VC 


DATA MOTVC 1/141*0 / 

DATA GCLAT /1HU / 

DATA GDLAT /O / 

DATA LONG /1HU / 

DATA MOT VC 2/9*0 / 

DATA GCRAD /1HU / 

DATA MOTVC3/148*0 / 

DATA ISV / 1 / 

DATA I SV3 /O / 

DATA DIARP /3*0.0 / 

DATA DR AGP /3*0.0 / 

DATA DRAGPT/0.0 / 

DATA FAXBP /3+0.0 / 

DATA CDP /3*0 .0 / 

DATA 01AMP /3*0 .0 / 

DATA DRGPK /3*0.0 / 

DATA DRGPP /3*0 .0 / 

DATA VELAP /0.0 / 

DATA DRGPS /3*0.0 / 

DATA IDRGP /3*0 / 

DATA PARIF /3*0.0 / 

DATA M0TVC4 /26*0/ 


PHZVC 


DATA ALTMAX/1.E20 / 

DATA ALTMIN/— 5000 • / 

DATA MAXTIM/1.0E10 / 

DATA EVTF /O / 

DATA FESN /100 / 

DATA IESNfPHZFfPIF* 14/4*0/ 
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c 

c 

c 


DATA SAVESN /O/ 


SPECAL 


C 

c 

c 


c 

c 

c 

c 


c 

c 

c 


DATA SPECAl/24*0 / 

TGOVC 


DATA 

FUXN 

/10*0 


DA TA 

PCTGO 

/ .9 

/ 

DATA 

SAVE 

/80*0/ 


DATA 

TCO 

/10.0E10/ 

DATA 

TIMX, 

ESN/2 *0/ 


DATA 

IEVNT 

/10*0 


DATA 

ISZEV 

/O 


DATA 

NXEVT 

/3 

/ 

DATA 

15 

/O 


DATA 

GUXN 

/10*0/ 


DATA 

TIMY 

/o/ 



/ 


/ 

/ 

/ 


TRACKC 


DATA 

CTKLAT 

/10*0.0 

/ 



DATA 

CTKLON 

/10*0.0 

/ 



DATA 

ELEV 

/10*0.0 

/ 



DATA 

LKA 

/10*0.0 

/ 



DATA 

LKB 

/10*0.0 

/ 



DATA 

PGT /90*0.0 

/ 



DATA 

SLNTRG 

710*0.0 

/ 



DATA 

SL0S1 

/10*0 

/ 



DATA 

SL0S2 

/10*0 

/ 



DATA 

SL0S3 

/10*0 

/ 



DATA 

STKLAT 

/10*0.0 

/ 



DATA 

STKLON 

/10*0.0 

/ 



DATA 

TKLATC 

/10*0.0 

/ 



DATA 

TKRAD 

/10*0.0 

/ 



DATA 

TRKAZM 

/10*0.0 

/ 



DATA 

TRKGLT 

/28. 22655 ,28. 

41338 

,26.62278 ,34.8155 ,6^0.0 

/ 

DATA 

TRKHIT 

/49.0 ,34. 

0 

,45.0 ,22.0 ,6*0 

/ 

DATA 

TRKLON 

/279. 40002, 279 

.40723 

,281.65167,283.64607,6*0.0/ 


DATA 

TRKXRI 

/30*0 .0 

/ 



DATA 

JTKFLG 

/6*0 

/ 



DATA 

TRKNAH 

/10HP ATRICK 

, 1 OHC APE KEN , 10HGRAN0 BAHA, 


1 


10HALTANT FLD 

,6*10H 

/ 


DATA 

NTRKS 

/10 

/ 



DATA 

NDUM 

/10 

/ 



DATA 

ELEMIN/0.0 

/ 




GUIDIC 

DATA DTG /1.0 / 
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non ooo ooonooo non non 


DATA GVRI /10*0 
DATA IGF /6*0 


/ 

/ 


GUIDVC 

DATA GTIME /1.E10 / 

DATA GUIDV1 /35*0 / 

TARGVC 

DATA TARGV1 /133*0/ 

MVEHF(l) TURNS ON THE INTEGRATION OF TARGET VEHICLE 
MVEHF ( 2 ) SPECIFIES TARGET INITIALIZATION 


DYTEM 

DATA NMAX, /AO / 

LOCAL 

DATA LOC ALI /145*0 .0 / 

DATA MSKFLG /400000000000000000008 , 377777777777777777776/ 

DATA ACOBS /16*0/ 

DATA ALTIO /O/ 

KROIC 

DATA KROIC I /13*0/ 

DATA GAM / 1.0, .5, .416666666666667,. 375, 

* .348611111111111, .329861111111111, .31 5591931216931 , 

* .304224537037037, .294868000440917, .286975446428571, 

* .280189596443937, .274265540031599, .269028846773649, 

* .5, .166666666666667, .125, .105555555555556, 

* 0.375E-2 f 8 .56150793650793E— 2 , 7.95717592592593E-2 , 

* 7.48522<>27689594E-2, 7. 103298611 1 11 IlE-2 , 6.78584 998463470E-2 , 

* 6 .5 1646 2 053571 43E— 2 » 6.28403190954034E-2, 6. 080 7479291 54 <>4E -2/ 

DATA GAS/1 .0 . 5 , -8 .3 33333 33 333 333E -2 »— 4. 16666666666667E-2 , 

* — 2 .63888888 88 888 9E— 2 , -1.R75E-2, 

* -1.42691798941799E-2, -1.13673941798942E-2, 

* —9 . 3565365961 1993E— 3 , -7.89255401234568E-3, 

* —6.7858499846347 IE— 3 , -5.92405641233767E-3, 

* -5.2366O325795029E-3/ 

DATA 

1 PT / 1., 2., 4., 8. 

2 , 16., 32., 64., 128. 

3 , 256., 512. ,1024. ,2048. 

4 ,4096. ,8192. / 

DATA CETPdl ,1 = 1,13)/ 3 .33333330E-01 , 2. 5000000 OE -01, 

1 1 • 13636360E— 0 1 , 6.73076^30E-02 , 4 .60526330E-02 , 3.43749980E-02 , 

2 2.71381590E— 02, 2.22547310E-02, 1.87484580E-02, 1.61 123220E-02 , 
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3 1.406030006-02, 1.24197060E-02, I .108021706-02/ 

DATA (ETP(I) *1=14,26)/ 2 .000000006-01 t 4.00000000E-01 » 

1 3.40909090E-01, 2.01923080E-01 * 1 .381579006-01 * 1 . 03124990E-0 1 * 

2 8 • 14 144780 E— 02* 6.67641930E-02, 5 .624537306-02 t 4.833696706-02, 

3 4.21809010E-02, 3.72591170E-02, 3.324065166-62/ 

DATA (6TP(I) *1=27,39)/ 1 .42857140E-01 , 2.85714280E -01 , 

1 3.428571406-01, 3 .46153840E-0 1, 2 .456140406-01 * 1. 87500000F-01 , 

2 1.503036506-01* 1.246264906-01* 1.058736406-01* 9.158583206-02* 

3 8.03445710E-02, 7.1 2783 130E-02, 6.382205106-02/ 

DATA (ETPf I) ,1=40,52)/ 1.111111106-01* 2.22 222220E-01, 

1 2.65714280E-01* 2 .53968250E-01 * 3.070175406-01* 2. 50000000E-0 1 , 

2 2 .08 75 5060 E— 01 * 1.78037850E-01, 1.543990606-01, 1 .35682710F-01 , 

3 1 • 20516850E— 01 , 1 .07997450E-0 1, 9.750590806-02/ 

DATA (ETP(I) ,1=53,65)/ 9.090909106-02, 1 .818 18 180E -01 , 

1 2.424242406-01, 2.424242406-01, 1.731601706-01, 2.500000006-01, 

2 2.27732800E-01, 2 .05428290E-01, 1 .85278 880E-01 , 1 .676080506-01 , 

3 1 .52231820E-01, 1 .38853860E-01 , 1 .271816206-01/ 

DATA (ETPCI) ,1=66,78)/ 7.69230760E-02 , 1.538461506-01, 

1 2.09790210E— 01, 2.23776220E-01, 1 .86480190E-01 , 1. 118881 10E-01 , 

2 1.91295550E— 01, 1 .91733070E-01 , 1.852788806-01, 1.75988460E-01 , 

3 1.657635306-01, 1.555163306-01, 1.456807706-01/ 

DATA (ETPCI) ,1=79,91)/ 6.666666606-02, 1.333333306-01, 

1 1.846153806-01, 2.05128200E-01, 1 .86480190F-01 , 1.34265730E-01 , 

2 6.961926906-02, 1.394422306-01, 1.520236906-01, 1.564341906-01, 

3 1.56012730E— 01, 1.52787970E-01 , 1 .47993 160E -01/ 

DATA (ETP(I) *1=92,104)/ 5. 882352906-02, 1. 176470606-01 , 

1 1.6470 58 80E-01, 1 .88235290E-0 1 , 1 .809954706-01 , 1 . 44796380E-0 1 , 

2 9.21431500E— 02, 4.21225830E-02, 9.77295190E-0 2, 1. 149312406-01 , 

3 1.253673806-01, 1 .30961120E-01, 1.331938506-01/ 

DATA (ETPCI), 1*105, 117)/ 5.263157906-02, 1.05263160E-01, 

1 1.486068106-01, 1 .7337461Q6-01 , 1.733746106-^01, 1.48606810E-01 , 

2 1.066920706-01, 6.096689^06-02, 2 .49410030E-02 , 6.63064840F-02 , 

3 8.357825306-02, 9.629493906-02, 1.051530406-01/ 

DATA (6TP( I) *1=118,130)/ 4.761904806-02, 9. 52380950E-02, 

1 1.353383506-01, 1.604010006-01, 1.651186806-01, 1 .48606810E-01 , 

2 1.155830806-01, 7.548282406-02, 3 .919300506-02 , 1.451592806-02* 

3 4.377908606-02, 5.884690806-02, 7.14002090F-02/ 

DATA (ETP(I) *1=131,143)/ 4.347B2610E-02 « 8.69565210E-02, 

1 1.242236006-01, 1.490683206-01, 1 .56914020E-01 , 1 . 46453090E-01 , 

2 1 .206084306—01, 8 .61488760E-02, 5 .16893250E-02 , 2.46139640 E-02, 

3 8.330880306-03* 2.82465160E-02* 4.032011806-02/ 

DATA (6TPU ) ,1 = 144,156 )/ 4.000000006-02, 7. 9999999 0E -02, 

1 1.14782610E-01, 1 .39130430^-0 1 , 1 .49068 320E-01 , 1.43105590F-01 , 

2 1.230205906-01, 9 .3729977OE-02, 6 .20271 9006-02 , 3.44595500E-02 , 

3 1.51622020E— 02, 4.725881 20E-03, 1 .78691 420E -02/ 

DATA (ETPCI) ,1=157,169)/ 3.703703706-02, 7.40740740E-02 , 

1 1.066666706-01, 1.303703706-01, 1.417069206-01, 1 . 391 30430E-01 , 

2 1.236715006-01, 9.893719806-02* 7.029748206-02* 4. 33935080F-02 , 

3 2. 246252206-02, 9.18921340E-03, 2.654661706-03/ 

DATA FAC /I . ,.5/ 

DATA FPP1 /.I / 

DATA FPP01 /.Ol / 
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non 


DATA FPP25 / .25 / 

KROVC 


C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


DATA KROVC1 /27*0/ 
DATA KBIT2 /98/ 

DATA KPC /0/ 

DATA KDD /0/ 

DATA KDS /O/ 

DATA KMAXO /?/ 

DATA KOM /O/ 

DATA KGMAX /12/ 

DATA KROVC 2 /e*0/ 

DATA RNDINT /2.8E-14/ 
DATA KPCVC3 /5o*0/ 


NPC(l) 
NPC 1 2 ) 
NPC ( 3 ) 
NPC(4) 
NPC (5) 
NPC(6) 
NPC ( 7) 
NPC (8) 
NPC ( 9 ) 
NPC ( 10 ) 
NPCU1 ) 
NPC( 12 } 
NPC ( 13 ) 
NPC ( 14 ) 
NPC ( 15) 
NP C ( 1 6 ) 
NPC ( 17 ) 
NPC ( 18 ) 
NPC ( 19 ) 
NPC ( 20 ) 
NPC (21 ) 
NPC ( 22 ) 
NPC ( 23 ) 
NPC ( 24 ) 
NPC (25) 
NPC ( 26 ) 
NPC ( 27 ) 
NPC ( 28 ) 
NPC (29) 
NPC (30) 
NPC ( 31 ) 
NPC ( 32 ) 
NPC ( 33 ) 
NPC ( 34 ) 
NPC ( 35 ) 


FLAGS CONIC CALCULATION 
FLAGS INTEGRATION SCHEME 
FLAGS INITIAL VELOCITY INPUT OPTION 
FLAGS INITIAL POSITION INPUT OPTION 
FLAGS ATMOSPHERE MODEL 
FLAGS ATMOSPHERE WINDS 
FLAGS ACCELERATION LIMIT FY THROTTLING 
FLAGS AERODYNAMIC COEFFICIENT TYPE 
FLAGS ENGINE TYPE 
FLAGS STATIC TRIM 

FLAGS ENVIRONMENTAL INEQUALITY CONSTRAINTS 

FLAGS CROSS AND DOWNRANGE CALCULATION 

FLAGS PROPELLANT JETTISON OPTION 

FLAGS HOLD DOWN OPTION 

FLAGS HEATING RATE 

FLAGS EARTH MODEL 

FLAGS MASS FRACTION JETTISON OPTION 
FLAGS TRAJECTORY TERMINATION 
FLAGS INPUT CONDITIONS PRINTOUT 
FLAGS DT MODEL 

FLAGS FLOWRATE METHOD FOR ROCKFT ENGINES 

FLAGS THROTTLING PARAMETER 

FLAGS DVMAR CALCULATION 

FLAGS GENERAL INTEGRATION 

FLAGS VELOCITY LOSSES 

FLAGS AERODYNAMIC HEATING INDICATORS 

FLAGS INDIVIDUAL ENGINE FLOWRATFS 

FLAGS TRACKING STATIONS 

FLAGS ANALYTIC IMPACT 

FLAGS WEIGHT AS FUNCTION OF TABLF LOOKUP 

flags vernal equanox calculations 

FLAGS PARACHUTE DRAG OPTION 
FLAGS BTL COORDINATES COMPUTATIONS 

- NOT USED 

- NOT USED 


END 
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DELTU: This is the main program of overlay (2,5). This 

routine determines the direction of search, based on the search/ 
optimization mode selected by user input. 



CALL WUCAL 
FOR PROBLEM 
WEIGHTING 


CALL GMAG 
TO CALCULATE 
GRADIENT 
MAGNITUDES 


C 


BRANCH ON SEARCH 
MODE (SRCHM) 


) 


/ CALL SDM, \ 

#201 / IS THE y_ 

y \ STEEPEST-DESCENT 
\ METHOD / 

/ WHICH IS THE \ 
\. 40 y \ CONJUGATE f 

^ \ GRADIENT METHOD / 

/ CALL DGMP2, \ 

>160 1,/ WICH IS THE V 

V L# \ DAVIDON METHOD / 

\ ° N P2 / 

/ CALL PGM, \ 
/ 80 1 r / WHICH IS THE \ 
\ PROJECTED f 
\ GRADIENT / 



CALL UNITDU 
TO UNITIZE THEl 
DIRECTION OF 
SEARCH 



CALL GABDD TO 

DETERMINE 

INEQUALITY 

STEP-SIZE 

BOUNDARIES 


> 
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DERVI: This routine initializes the time references and calls 

the routines to initialize the equations of motion. 



i 

i 
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DGAMLAM: This routine computes the rate of change in path 

angle and azimuth relative to the atmosphere. 



COMPUTE INERTIAL 
VELOCITY COMPONENTS 
IN THE G FRAME 

7 ~ * 

COMPUTE RATE OF 
CHANGE IN 
LATITUDE, 

ALTITUDE, AND 
INERTIAL VELOCITY 


COMPUTE RATE OF CHANGE 
IN THE ATMOSPHERIC 
RELATIVE VELOCITY 

7=r= 

I DGAMA, AND DLAMD 




DGM (G) : This routine computes the direction of search based 

on the Davidon deflected gradient method and is used to minimize 
a scalar-valued function whose gradient vector is given by G(l). 



COMPUTE A NEW DEFLECTION 
MATRIX BASED ON G(I), 

GP ( I ) , AND HESS (I) 



COMPUTE TIIE NEW DIRECTION 
OF SEARCH BASED ON THE 
NEW DEFLECTION MATRIX 
AND THE CURRENT GRADIENT 
VECTOR 

* 

STORE THE COMPUTED 
DIRECTION OF SEARCH IN 
PP(J) AND THE CURRENT 
GRADIENT VECTOR IN 
GP (J) 
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DGMP2: This routine computes the direction of search for 

minimizing P2 via the Davidon variable metric method. 
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DICT: This routine is not call explicity, but maps the 

variable names to core locations at the time overlay (1,0) is 
called. 



♦DECKfDICT 

SUBROUTINE DICT 
C*** DICT 

C DICT - DEFINES COMPUTATIONAL COMMONS 

C (IV - END) DICTIONARY VALUES DATA 

C*** THE DATA STATEMENTS IN THIS ROUTINE ARE STANDARDIZED 
C*** THERE MUST BE AN EVEN NUMBER OF INTFGERS IN COMMONS BETWEEN 
C IV AND END FOR THE UNIVAC 1108 DOUBLE PRECISION VERSION 
0 

COMMON/IV/ IV<2) 

C 

C COMPUTATIONAL DATA REGION 

C 

♦CALL AUXVC 
♦CALL CYCVC 
♦CALL DPGVC 
♦CALL DYNVC 
♦CALL INFVC 
♦CALL MNMMLT 
♦CALL HOLINC 
♦CALL MOTBL 
♦CALL MOTIC 
♦CALL MOTVC 
♦CALL PHZVC 
♦CALL SPECAL 
♦CALL TGOVC 
♦CALL TRACKC 
♦CALL GUIDIC 
♦CALL GUIDVC 
♦CALL TARGVC 

COMMON/END/END 

C 

£♦♦* COMMONS NOT INCLUDED IN DICTIONARY 

♦CALL DYTEM 
C 

£♦♦* COMPUTATIO N A L DATA DICTIONARY 

DIMENSION AUXVCM 126) 

DIMENSION AUXVC2(60) 

DIMENSION CYCVC 1( 16) 

DIMENSION DPGVC1 ( 126 ) 

DIMENSION DPGVC2(119) 

DIMENSION DYNVC 1 ( 18 ) 

DIMENSION INFVC1 ( 24) 

DIMENSION MNMML1 ( 17 ) 

DIMENSION H0LIN1 ( 51 ) 

DIMENSION MOTICK 126) 

DIMENSION M0TIC2 ( 126 ) 

DIMENSION M0TIC3< 1) 

DIMENSION MOTVC 1 ( 126 ) 

DIMENSION M0TVC2C 126) 
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c 


DIMENSION 

DIMENSION 

DIMENSION 

DIMENSION 

DIMENSION 

DIMENSION 

DIMENSION 

DIMENSION 

DIMENSION 

DIMENSION 

DIMENSION 


M0TVC3( 110) 
PHZVC1U0) 
SPECA1 ( 24) 
TG0VC1 ( 118 ) 
TRACKM126) 
TR ACK2 ( 126 ) 
TRACK3 ( 61) 
GUIDI1 ( 17) 
GUIDVl ( 37) 
TARGVK126) 
TARGV2 (7) 


EQUIVALENCE ( AUXVC1 ( 1 ) , ALPTOT ) 
EQUIVALENCE CAUXVC2 (1)»XR(3) ) 

EQUIVALENCE (CYCVCim.DELT ) 

EQUIVALENCE (DPGVC1 ( 1 ) , AL PHA ) 
EdUI VAlENCT TbPGVC2 ( 1 ) »KDG ( 2 ) ) 

EQUIVALENCE (DYNVC1 ( 1 ) ,DTIMR ( 1 ) ) 
EQUIVALENCE (INFVC1 (1 ) »ESNPRT ) 
EQUIVALENCE (MNMML1 (1 ) rONE ) 

EQUIVALENCE (H0LIN1 ( 1 ) , ALPARG ) 
EQUIVALENCE (M0TIC1 (1 ) , ALTIP ) 
EQUIVALENCE (M0TIC2 ( 1 ) ,LREFY ) 
EQUIVALENCE (M0TIC3 ( 1 ) *OATE ( 3 ) ) 

EQUIVALENCE (M0TVC1 ( 1 ) ,AHI ) 

EQUIVALENCE (M0TVC2 ( 1 ) , GXI (3 ) ) 

EQUIVALENCE CM0TVC3 ( 1 ) ,XLO (3 ) ) 

EQUIVALENCE (PH2VC1 ( 1 ) » ALTMAX ) 
EQUIVALENCE (SPEC A1 (1 ) ,SPECI (1) ) 
EQUIVALENCE (TG0VC1 ( 1 ) * FUXNU ) ) 

EQUIVALENCE (TRACK1 (1 ),CTKLAT(1 )) 
EQUIVALENCE (TRACK2U) ,PGT(77) ) 

EQUIVALENCE (TRACK3 ( 1 ) ,TRKL0N(3 ) ) 
EQUIVALENCE (GUIDI1 < 1 ) ,DTG ) 

EQUIVALENCE ( GUIDVl ( 1 ) t GTIME ) 
EQUIVALENCE ( TARGV1 ( 1 ) , ALTAT ) 
EQUIVALENCE (TARGV2 (1)»DAXPT(3) ) 


1 

2 

3 

4 

5 

6 

7 

8 
9 
0 
A 


/6HALPT0T 

♦6HARGV 

♦6HDWNRNG 

♦6HHYPVEL 

*6HlSAGRE 

»6HPGVEL 

,6HSEMJAX 

♦6HTIMTP 

♦6HTLHT7 

,6HTTLISP 

» 6HWTP2 


♦6HALTA 

.6HDECLIN 

t 6HDVCIR 

,6HlNt 

•6HLAN 

♦6HOALTOT 

,6HSIDEAL 

♦6HTLHT1 

♦6HTLHT8 

,6HUBAR 

t 6HWTP3 


DATA 

♦6HALTP 

♦6HDRAG 

.6HDVEXS 

t 6HlNCPCH 

♦6HL0NGIP 

♦6HREYN0 

,6HSSVIDL 

,6HTLHT2 

♦6HTLHT9 

♦6HVIP1 

*6HWTP4 


AUXVC1 

♦6HANGM0M 

♦6HLIFT 

,6HECCAN 

*6H INCYAW 

,6HMEAAN 

.6HRTP1 

,6HTHTP 

♦ 6HTLHT3 
♦6HTLHT10 
1 6 H VI P 2 

♦ 6HWTP5 


.6HAP0RAD 

,6HDGENV 

.6HECCEN 

♦6HIPNULL 

♦6HPERI0D 

♦6HRIP2 

♦6HTHTPL 

*6HTLHT4 

,6HTLPWT 

♦6HVIP3 

,6HWTP6 


♦6HAPVEL 

.6HDPRNG1 

♦6HENERGY 

,6HIYNULL 

,6HPGERAD 

,6HRIP3 

*6HTIM IP 

,6HTLHT5 

,6HTRUAN 

»6HVMU 

»6HWTP 7 


,6HARGP 

,6HDPRNG? 

♦6HGDLTIP 

♦6HIRANGE 

♦6HPGL0N 

♦6HRTASC 

,6HTIMSP 

»6HTLHT6 

,6HTRUNMX 

♦ 6HWTP1 

t 6HWTP8 
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B *6HWTP9 .6HWTP10 ,6HXISAV1 ,6HXISAV2 ,6HXISAV3,6HXMAX1 ,6HXMAX2 
C ,6HXMAX3 .6HXMAX4 ,6HXMAX5 ,6HXMAX6 ,6HXMAX7 ,6HXMAX8 ,6HXMAX9 

0 ,6HXMAX10,6HXMIN1 .6HXMIN2 ,6HXMIN3 .6HXMIN4 ,6HXMIN5 *6HXMIN6 

E *6HXMIN7 ,6HXMIN8 .6HXMIN9 ,6HXMIN10 ,6HCRRNG .6HYXMN1 ,6HYXMN2 
F » 6HYXMN3 ,6HYXMN4 ,6HYXMN5 ,6HYXMN6 ,6HYXMN7 ,6HYXMN8 .6HYXMN9 

G *6HYXMN10»6HYXMX 1 ,6HYXMX2 ,6HYXMX3 ,6HYXMX4 .6HYXMX5 ,6HYXMX6 

H »6HYXMX7 ,6HYXMX8 ,6HYXMX9 ,6HYXMX10,6HMACHDT,6HXR ,6HYR 

1 / 

C 

DATA AUXVC2 

1 /6HZR *6HXVE ,6HYVE f 6HZVE ,6HVXVE ,6HVYVE .6HVZVE 

2 »6HXS VE *6HYSVE .6HZSVE ,6HXSI ,6HYSI ,6HZSI ,6HSHADF 

3 » 6HSC0NF ♦6HSCL0CK T 6HXIVEIl,6HXIVE12 f 6HXIVE13»6HXI VE21 t 6HXIVE22 

4 y 6HXIVE 23 *6HX I VE 31 *6HXIVE32 *6H XI VE33 »6HURX ,6HURY ,6HURZ 

5 ,6HUTX ,6HUTY ,6HUTZ ,6HUNX ,6HUNY *6HUNZ ,6HRAS 

6 t 6HLANVE ,6HVCIRC ,6HGULIE .6HRASGM ,6HBTLX .6HBTLY » 6HBTLZ 

7 ,6HBTLXD ,6HBTLYD ,6HBTLZD »6HZETA 1 ,6HZETA2 ,6HZETA3 »6HNTRK 

8 »6HDPRGI1 »6HDPRGI2 »6HGCLTIP ? 6HGENV1 »6HGENV2 *6HPGENV »6HRGENV 

9 ,6HSGENV ♦6HDTRUAN ,6HDYNPD »6HBETAS 
I / 


DATA CYCVC1 

1 /6HDELT *6HDT »6HDTIME »6HDTH *6HDT0 »6HEN0IS f 6HTREF 

2 *6HIDTABl»6HIDTAB2r6HIDTAB3t6HIDTAB4*6HIDTAB5*6HIDTAB6f 6HIFLG 

3 ,6HCYCF »6HDELTT 

4 / 


DATA DPGVC1 

1 /6HALPHA »6HBETA » 6HBNKANG *6H ALPPC1 »6HALPPC2 »6HALPPC3 « 6HALPPC4 

2 t6HPETPCl »6HBETPC2 »6HBETPC3 ,6HBETPC4, 6HBNKPC1 *6HBNKPC2»6HBNKPC3 

3 *6HBNKPC4f 6HALPD OT »6HEETD0T » 6HPNKD0T t6HDALPHA»6HDBETA ,6HDBANK 

4 t6HDYAW *6HDPITCH»6HDR0LL *6HAB1 «6HAB2 »6HAB3 ,6HAB4 

5 *6HAB5 » 6MAB6 *6HAB7 »6HAB8 *6HAB9 »6HGB1 »6HGB2 

6 »6HGB3 »6HGB4 t 6HGB5 »6HGB6 »6HGB7 «6HGB8 ,6HGB9 

7 t 6HIA1 t 6HIA2 *6HIA3 *6HIA4 »6HIA5 »6HIA6 *6HIA7 

8 ,6HIA8 t 6HIA9 *6HIE11 ,6HIB12 *6HIB13 »6HIB21 »6HIB22 

9 .6HIB23 .6HIB31 ,6HIB32 ,6HIB33 ,6HIG1 ,6HIG2 ,6HIG3 

0 v 6HIG4 *6HIG5 ,6HIG6 *6HIG7 ,6HIG8 ,6HIG9 ,6HIL1 

A »6HIL2 *6HIL3 ,6HIL4 »6HIL5 ,6HIL6 »6HIL7 »6HIL8 

B *6HIL9 ,6HLB1 ,6HLB2 ,6HLB3 ,6HLB4 ,6HLB5 ,6HLB6 

C ? 6HLB7 *6HLB8 ,6HLB9 ,6HHAPG1 ,6HHARG2 .6HHARG3 ,6HHERR0R 
D ,6HHERRC1R*6HHERR0R,6HIDGF1 ,6HIDGF2 ,6HIDGF3 ,6HIGUID1 ,6HIGUID2 
E *6HIGWID3*6RIGUID4t6HIGUlD5 »6H IGUID6»6HIGUID7*6HIGUID8t 6HIGUID9 
F ,6HIGUIDl f 6HIGUIDl »6HIGUIDl,6HIGUIDl t 6HIGUIDl,6HIGUIDl*6HIGUIDl 
G ,6HIGUIDI,6HIGUID1,6HIGUID1,6HIGUID2,6HIGUID2,6HIGUID2,6HIGUID2 

H 1 6HIGUID2 ?6HIGUID2 » 6HIVCRT1 »6H IVCRT2 »8HI VCRT3 »6HIVE TA ,6HKDG1 

DATA DPGVC2 
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1 /6HKDG2 ,6HKDG3 *6HKRG1 » 6HKRG2 »GHKRG3 ,6HYAWI ,6HP1TI 

2 ,6HR0LI »6HYAWR ,6HPITR ,6HR0LR ,6HR0LPC1,6HR0LPC2,6HR0LPC3 

3 ,6HR0LPC4,6HYAWPCl,6HYAWPC2 f 6HYAWPC3,6HYAWPC4,6HPITPCl»6HPITPC2 

4 ,6HPITPC3,6HPITPC4,6HR0LBD ,6HPITB0 ,6HYAWBD ,6HYAWID ,6HPITID 

5 » 6HR0LI0 »6HYAWRD ,6HPITRD ,6HR0LRD ,6HXY0MA1 ,6HXY0MA2,6HXY0MA3 

6 ,6HXYOMBl t 6HXYOMB2,6HXYOMB3,6HXYOMSl»6HXYOM52,6HXYOMS3,6HXYOMH 

7 *6HXY0M12 *6HXY0M 13 ,6HXY0M21 ,6HXY0M22 ,6HXY0M23 ,6HXY0M3 1 , 6HXY0M32 

8 ,6HXY0M33,6HXY0MElt6HXY0ME2,6HXY0ME3,6HGUI0XH*6HGUIDYH,6HGUIDZH 

9 t 6HALPHI »6HBETA I ,6HBANKI *6HIIV11 *6HIIV12 .6HIIV13 ,6HIIV21 

0 *6HIIV22 ,6HIIV23 ,6HIIV31 ,6HIIV32 ,6HIIV33 ,6HIVB11 .6HIVE12 

A ,6HIVB13 »6HI VB21 ,6HIVB22 .6HIVB23 ,6HIVB31 ,6HIVB32 ,6HIVB33 

B tBHNOEPVS *6HDPVLS1 »6HDPVLS2 * 6HDPVLS3 »6HDPVLS4»6HDPTLS 1 *6HDPTL S2 
C ,6HDPTLS3,6HDPTLS4,6HMAXITSt6HNITS ,6HPERTS1 ,6HPERTS2,6HPERTS3 
D ,6HPERTS4,6HEN0MSl,6HEN0MS2,6HEN0MS3,6HEN0MS4t6HESl ,6HES2 

E *6HES3 » 6HES4 ,6HUN0MS1 ,6HUN0MS2 .6HUN0MS3 ,6HUN0MS4,6HU$ 1 
F *6HUS2 ,6HUS3 ,6HUS4 ,6HEN0RMS*6HISTEPS,6HNAMVSl,6HNAMVS2 
G ,6HNAMVS3 ,6HNAMVS4,6HNAMVS5 *6HNAMVS6 ,6HNAMVS7,6HNAMVS8 , 6HGS AI TS 

1 / 


DATA DYNVC1 

1 /6HDTIMR 1 »6HDTIMR2 »6HDTIMR3»6HDTIMR4»6HTDURP .6HTIMES *6HTIMRF 1 

2 ,6HTIMRF2,6HTIMRF3,6HTIMRF4,6HTREFP ,6HTREFS »6HNDISC .6HNPASS 

3 *6HERRINT »6HDLTM IN »6HDLTMAX»6HEPS INT 

4 / 


DATA INF VC 1 

1 /6HESNPRT »6HEXTRAP »6HLPRNT ,6HPINC ,6HPRNC t 6HFIDl ,6HFID2 

2 .6HINFF ,6HIPRNTe,6HlPRNTR,6HPST0P ,6HTITLE 1 .6HTITLE2 , 6HTITLE 3 

3 ,6HTITLE4,6HTITLE5,6HTITLE6,6HTITLE7,6HTITLE8,6HTITLE9,6HTITLE1 

4 »6HSFID1 *6HSF ID2 ,6HNPINC 
5/ 


DATA MNMML1 

1 /6H0NE ,6HCADPNM,6HCADYNM,6HCANM f 6HCDDPNM ,6HCDDYNM , 6HCDNM 

2 »6HCLDPNM*6HCLNM ,6HCMANM ,6HCMDPNM .6HCNANM ,6HCNDPNM ,6HCWBNM 

3 ,6HCWDYNM*6HCYBNM ,6HCYDYNM 

4 / 


DATA HOL INI 

1 /6HALPARG»6HBETARG *6HBNKARG*6HEYAARG»6HDGF1 *6HDGF2 ,6HDGF3 

2 * 6HGD ER VI , £WGtS)§RV2 , 6HGtftMV3 , V4 , 6 WdYRV 5 ,6H6i>E RV6 , 6HGDFR V7 

3 *6HGDERV8 ♦6HG0ERV9»6HG>DERV0 »6HM0NF1 ,6HM0NF2 .6HM0NF3 »6HM0NX! 

4 .6HM0NX2 ,6HM0NX3 ,6HM0NX4 ,6HM0NX5 •6HM0NX6 ,6HM0NX7 f 6HM0NX8 

5 *6HM0NX9 ,6HM0NX10,6HM0NYI ,6HMONY2 .6HM0NY3 ,6HM0NY4 » 6HM0NY5 
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6 ,6HM0NY6 .6HM0NY7 ,6HM0NY8 .6HM0NY9 ,6HM0NY10,6HYAWARG,6HPITARG 

7 y6HR0LARG*6HDPVRSI y6HDPVRS2 »6HDP VRS3»6HDPVRS4,6HINVRS 1 y6HINVRS2 

8 1 6 HINVR S3 y 6HINVRS4 
8 / 


DATA 

DATA 

I V( 2 ) 
CST 

/6H* 

/6HCST 

/ 

» 6HCSH 

/ 

DATA 

ATEMT 

/6HATEMT 

♦6HATEMM 

/ 

DATA 

PREST 

/6HPREST 

« 6HPRE SM 

/ 

DATA 

VWUT 

/6HVWUT 

♦6HVWUM 

/ 

DATA 

VWVT 

/6HVWVT 

» 6HVWVM 

/ 

DATA 

VWWT 

/6HVWWT 

1 6 HVWWM 

/ 

DATA 

AZWT 

/6HAZWT 

y6HAZWM 

/ 

DATA 

VWT 

/6HVWT 

y 6HVWM 

/ 

DATA 

ALPHAT/6HALP HAT » 6HAL PHAM/ 

DATA 

BETAT 

/6HBETAT 

y 6HBE TAM 

/ 

DATA 

BANKT 

/6HBANKT 

y 6HBANKM 

/ 

DATA 

YAWT 

/6HYAWT 

y6HYAWM 

/ 

DATA 

PITT 

/6HPITT 

y6HPITM 

/ 

DATA 

ROLT 

/6HR0LT 

y 6HR0LM 

/ 

DATA 

CDT 

/6HCDT 

y 6HCDM 

/ 

DATA 

CLT 

/6HCLT 

yfeHCLM 

/ 

DATA 

CAT 

/6HCAT 

y 6HCAM 

/ 

DATA 

CNAT 

/6HCNAT 

,6HCNAM 

/ 

DATA 

CYBT 

/6HCYBT 

t 6HCYBM 

/ 

DATA 

CMAT 

/6HCMAT 

y 6HCMAM 

/ 

DATA 

CWBT 

/6HCWBT 

, 6HCWBM 

/ 

DATA 

CMDPT 

/6HCMDPT 

y 6HCMDPM 

/ 

DATA 

CADPT 

/6HCADPT 

y 6HCADPM 

/ 

DATA 

CNDPT 

/6HCNDPT 

y 6HCNDPM 

/ 

DATA 

CWDYT 

/6HCWDYT 

y 6HCWDYM 

/ 

DATA 

CADYT 

/6HCA0YT 

y 6HCADYM 

/ 

DATA 

CYDYT 

/6HCYDYT 

y 6HCYDYM 

/ 

DATA 

XCGT 

/6HXCGT 

y 6HXCGM 

/ 

DATA 

YCGT 

/6HYCGT 

t 6HYCGM 

/ 

DATA 

ZCGT 

/6HZCGT 

y6HZCGM 

/ 


DATA TVCIT 

1/6HTVC1T y 6HTVC1M , 6HTVC2T ,6HTVC2M ,6HTVC3T ,6HTVC3M 
2 *6HTVC4T ,6HTVC4M t 6HTVC5T ,6HTVC5M ,6HTVC6T ,6HTVC6M 
3.6HTVC7T ,6HTVC7M *6HTVC8T ,6HTVC8M ,6HTVC9T ,6HTVC9M 

4y6HTVC10T,6HTVCl0M,6HTVCllT,6HTVCllM,6HTVC12T,6HTVC12M 
5y6HTVC13T,6HTVC13M,6HTVC14Ty6HTVC14M,6HTVC15T*6HTVC15M 
6 / 


1/6HW01 T ,6HWD1M 
2.6HWD4T ,6HWD4M 
3 t 6MWD7T t 6HWD7M 
4*6HWD10T , 6HWD10M 


DATA WDIT 

♦6HWD2T t6HWD2M »6HHD3T 
♦ 6HWD5T T 6HWD5M ,6HWD6T 
»6HWD8T t 6HWD8M ,6HWP9T 
»6HWD1 1 T .6HWD11M *6HWD12T 


t 6HWD3M 
♦6HWD6M 
.6HWD9M 
t 6HWD 12M 
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5 *6HWD13T ,6HWD13M ,6HWD14T ,6HWD14M *6HW015T 
6 / 


DATA AEIT 

1/6HAE1T *6HAE1M .6HAE2T y6HAE2M ,6HAE3T 
2 »6HAE4T ,6HAE4M ,6HAE5T ,6HAE5M ,6HAE6T 
3y6HAE7T ,6HAE7M ,6HAE8T ,6HAE8M ,6HAE9T 
7,6HAE10T ,6HAE10M ,6HAE11T ,6HAEI1M ,6HAE12T 
5.6HAE13T ,6HAE13M *6HAF1AT »6HAE1AM .6HAE15T 
6 / 


DATA 

FLIT 

/6HFL1T 

y 6HFL1M 

/ 

DATA 

FL2T 

/6HFL2T 

,6HFL2M 

/ 

DATA 

FL3T 

/6HFL3T 

y 6HFL3M 

/ 

DATA 

XREFT 

/6HXREFT 

» 6HXREFM 

/ 

DATA 

YREFT 

/6HYREFT 

, 6HYREFM 

/ 

DATA 

ZREFT 

/6HZREFT 

1 6 HZREFM 

/ 

DATA 

CLDPT 

/6HCLDPT 

f 6HCLDPM 

/ 

DATA 

CDDPT 

/6HCDDPT 

1 6 HCDDPM 

/ 

DATA 

CDDYT 

/6HCDDYT 

y 6HCDDYM 

/ 

DATA 

DENST 

/6HDENST 

♦6HDENSM 

/ 

DATA 

HTRTT 

/6KHTRTT 

1 6 HHTRTM 

/ 

DATA 

ETAT 

/6HETAT 

♦6HETAM 

/ 

DATA 

CAOT 

/6HCA0T 

1 6HCA0M 

/ 

DATA 

CNOT 

/6HCN0T 

1 6HCN0M 

/ 

DATA 

CYOT 

/6HCY0T 

v 6HCY0M 

/ 




DATA 

PIITyYIIT 

1/6HPI1T 

t6HP I I M 

,6HPI2T 

,6HP I2M 

,6HPI3T 

2,6HPIAT 

yfcHPIAM 

,6HPI5T 

*6HP I5M 

,6HPI6T 

3,6HPI7T 

t 6HPI7M 

,6HP I8T 

,6HP IPM 

y 6HP 1 9T 

4.6HPI10T 

,6HPI10M 

y6HP HIT 

♦6HPI11M 

,6HPI12T 

5 y6HP I13T 

,6HPI13M 

,6HPI14T 

♦6HPI14M 

,6HPI 1 5T 

6,6HYI1T 

♦6HYI1M 

,6HYI2T 

♦6HYI2M 

y 6HYI3T 

7,6HYIAT 

»6HYIAM 

♦6HYI5T 

,6HYI5M 

,6HYI6T 

8,6HYI7T 

,6HYI7M 

»6HYI8T 

,6HYI8M 

♦6HYI9T 

9 *6HYI 10T 

y 6HYI I OM 

,6HYI 1 1 T 

,6HYII1M 

,6HYI12T 

0,6HYI13T 

,6HYI 13M 

1 6 HYI 1AT 

,6HYI14M 

,6HYI15T 


DATA CDOT /6HCD0T ,6HCD0M / 
DATA CLOT /6HCL0T ,6HCL0M / 
DATA CWOT /6HCW0T ,6HCW0M / 
DATA CMOT /6HCM0T ,6HCM0M / 
DATA WUA1T /6HWUA1T * 6HWUA1M / 
DATA WUA2T /6HWUA2T t 6HWUA2M / 
DATA DENKT /6HDENKT ,6HDENKM / 
DATA GDF1T /6HGDF1T ,6HGDF1M / 
DATA GDF2T /6HGDF2T »6HGDF2M / 
DATA GDF3T /6HGDF3T , 6HGDF3M / 
DATA GN0M1T/6HGN0M1T» 6HGN0M1M/ 
DATA GN0M2T/6HGN0M2T y 6HGN0M2M/ 
DATA GNOM3T/6HGNOM3T * 6HGN0M3M/ 
DATA GNMXlT/6HGNMXlTy6HGNMXlM/ 
DATA GNMX2T/6HGNMX2T * 6HGNMX2M/ 
DATA GNMX3T/6HGNMX3T 1 6HGNMX3M/ 
DATA GNMN1T/6HGNMN1T »6HGNMN1M/ 


»6HWD15M 


*6HAE3M 
,6HAE6M 
♦6HAE9H 
t 6HAE12M 
»6HAE 15M 


»6HPI3M 

y6HP I 6M 

,6HPI9M 
»6HPI 12M 
♦6HPI15M 
.6HYI3M 
♦6HYI6M 
♦6HYI9M 
,6HYI12M 
•6HYI15M / 
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DATA GNMN2T/6HGNMN2T »6HGNMN2M/ 

DATA GNMN3T/6HGNMN3T y 6HGNMN3M/ 

DATA GENV1T/6HGENV1T.6HGENV1M/ 

DATA GENV2T/6HGENV2T »6HGENV2M/ 

OATA FMASST/6HFMASST »6H FMASSH/ 

DATA ZLALPT/6HZLALPT,6HZtALPM/ 

DATA CAIOT /6HCAI0T ,6HCAI0M / 

DATA MGT1T /6HWGT1T .6HWGT1M / 

DATA WGT2T /6HWGT2T ,6HWGT2M / 

DATA WGTD1T/6HWGTD1T,6HWGTD1M/ 

DATA WGTD2T/6HWGTD2T « 6HWGTD2M/ 

DATA CDPIT /6HCDP1T ,6HCDP1M ,6HCDP2T .6HCDP2M ,6HCDP3T .6HCDP3M / 
DATA HTRT1T/6HHTR TIT *6HHTRT1 M/ 


DATA MOTIC1 

1 /6HALTIP ,6HALTREF,6HLATREF,6HL0NREF f 6HAZREF ,6HTIMREF ,6HARP1 

2 * 6HARP2 f 6HAR P3 ,6HARP4 ,6HARP5 ,6HARP6 ,6HARP7 ,6HARP8 

3 f 6HARP9 »6MARP10 ,6HASMAX »6H ATMSK1 ,6HATMSK2,6HAZWB ? 6HAZL 

4 » 6HLATL »6HLONL »6HCLCDMX *6HDFTA ,6HDESN1 ,6HDESN2 y6HDESN3 

5 ,6HDESNE »6HDVIMAG,6HDVMAR ,6HDVPCT ,6HETAPCly6HETAPC2y6HETAPC3 

6 »6HETAPC4,6HETA ,6HGINT1 ,6HGINT2 ,6HGINT3 ,6HGINT4 ,6HGINT5 

7 ,6HGINT6 »6HGINT7 ,6HGINT8 ,6HGINT9 »6HG INTIO »6HG0 ,6HGXP1 

8 y 6HGXP2 ,6HGXP3 ,6HGXP4 ,6>HGXP5 ,6HGXP6 ,6HGXP7 r 6HGXP8 

9 ,6HGXP9 .6HGXP10 ,6HGXP11 »6HGXP12 »4HGXP13 t 6HGXP14 ,6HGXP15 

0 « 6HGYP1 t 6HGYP2 .6HGYP3 ,6HGYP4 ,6HGYP5 .6HGYP6 ,6HGYP7 

A » 6HGYP8 »6HGYP9 ,6HGYP10 »6HGYP11 ,6HGYP12 ,6HGYP13 ,6HGYP14 

B *6HGYP15 t 6HGZPl ? 6HGZP2 ,6HGZP3 ,6HGZP4 ,6HGZP5 f 6HGZP6 

C »6HGZP7 ,6HGZP8 ,6HGZP9 ,6HGZP10 yAHGZPll ,6HGZP12 ,6HGZP13 

D » 6HGZP14 .6HGZPI5 ,6HHEATK1 , 6HHEATK2 *6HHEATK3 y6HALTIT0 ,6HHR ATI 
E » 6HHRAT2 ,6HHRAT3 » 6HHRAT4 ,6HHRAT5 ,6HHRAT6 ,6HHRAT7 .6HHRAT8 
F ,6HHRAT9 ,6HHR AT10,6HISPV1 t 6HISPV2 ,6HISPV3 ,6HISPV4 ,6HISPV5 
G »6HISPV6 t 6HI SPV7 ,6HISPV8 ,6HISPV9 ,6HISPV10,6HISPVlly6HISPV12 
H *6HISPV13»6HISPV14,6HISPV15 t 6HJ2 ,6HJ3 *6HJ4 f 6HLREF 

1 / 

DATA MOTIC2 

1 /6HLREFY f fcHMU » 6H0MEGA ,6HPGCL AT,6HPSL .6HPWPR0P, 6HRH0SL 

2 * 6HPN ,6HRE ,6HRP ,6HSREF ,<»HTSL ,6HWGTSG »6HWJETT 

3 y 6HWPLD y6HWPR0PI f 6HWEIC0Ny6HXREF t 6HYREF y6HZREF y6HAEXP 

4 y 6HCINF f 6HVINF I y6HIENGAl ,6HIENGA2,6HIENGA3y6HIENGA4y6HIENGA5 

5 »6HIFNGA6,6HIENGA7y6HIENGA8 y6HIENGA9 ,6HI ENGA 1 ,6HIENGA1 ,6HIENGA1 

6 »6HIENGA1,6HIENGA1*6HIENGA1 y6H IENGT1 y6HIENGT2,6HIENGT3 »6HIENGT4 

7 ,6HIENGT5,6HIENGT6,6HIENGT7,6HIENGT8,6HIENGT9*6HIENGT1,6HIENGT1 

8 y6HIENGTl,6HIEK'GTl,6HIENGTl,6HIENGTl,6HITAPl .6HITAP2 t 6HITAP3 

9 .6HITAP4 ,6HITAP5 t 6HITAP6 t 6HITAP7 *6HTTAP8 y6HITAP9 t 6HITAP10 

0 y 6HIWPF1 y 6HIWPF 2 y6HIWPF3 y6HIWPF4 .6HIWPF5 y6HIWPF6 t 6HIWPF7 

A t 6HIWPF8 y6H IWPF 9 ,6HIWPF10»6HIWPFn y6HIWPFI 2 y6HIWP FI 3 , 6HIWPF14 
B t 6HIWPF15y6HNENG f 6HNE0Sl y6HNEQS2 t 6HNECS3 y6HNEWSTG,6HNPCl 

C y&HNPC2 y6HNPC3 t 6HNPC4 ,6HNPC5 ,6HNPC6 ,6HNPC7 »6HNPC8 

D y 6HNPC9 y 6HNPC1 0 t 6HNPCll t 6HNPC12 ,6HNPC13 ,6HNPC14 y6HNPC15 

E *6HNPC16 ,6HNPC17 ,6HNPC18 t 6HNPC19 t 6HNPC20 # 6HNPC21 ,6HNPC22 
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non 


F 

» 6HNPC23 

,6HNPC24 

.6HNPC25 

» 6HNPC26 

f 6HNPC27 

,6HNPC28 

t6HNPC29 

G 

,6HNPC30 

,6HNPC31 

♦6HNPC32 

,6HNPC33 

»6HNPC34 

*6HNPC35 

t 6HNPC36 

H 

T 

*6HGHA 

/ 

,6HGHAS 

»6HDECL 

1 6 HTRPM 

,6HDVMARR,6HDATE1 

♦6HDATE2 

X 



DATA 

MOTIC3 




1 

I 

/6HDATE3 

/ 


DATA 

MOT VC 1 




1 

/6HAHI 

»6HAHID 

,6HAMXB 

»6HAMYB 

»6HAMZB 

♦6HASM 

,6HASMG 

2 

» 6HATEM 

,6HAXB 

»6HAYB 

,6HAZB 

t 6HASXI 

»6HASYI 

* 6HASZ I 

3 

» 6HAXI 

,6HAYI 

»6HAZ I 

,6HAXL1 

♦feHAXtl 

♦6HAXL3 

♦ 6HCA 

4 

1 6 HCN 

t 6HCD 

,6HCL 

,6HCDDP 

,6HC($DY 

t6HCLDP 

,6HCM 

5 

» 6HCMDP 

♦6HCADP 

,6HCNDP 

,6HCWDY 

♦6HCADY 

*6HCYDY 

»6HCY 

6 

,6HCW 

,6HCIP1 

»6HCIP2 

.6HCIP3 

,6HCIP4 

,6HCIP5 

,6HCIP6 

7 

t 6HCIP7 

,6HCIP8 

,6HCIP9 

t 6HCIP10 

,6HCIP11 

♦6HCIP12 

,6HCIP13 

6 

,6HCIP14 

*6HCIP1 5 

,6HCIY1 

»6HCI Y2 

,6HCIY3 

,6HCIY4 

» 6HC I Y5 

9 

♦6HCIY6 

»6HC IY7 

♦6HCIY8 

*6HCIY9 

.6HCIY10 

,6HCIY11 

*6HCI Y1 2 

0 

♦6HCIY13 

»6HC IY1 4 

♦6HCIY15 

,6HCS 

♦6HDAX 

♦ 6HDAY 

,6HDAZ 

A 

1 6 HDE 0 

,6HDE1 

t 6HDE2 

*6HDE3 

*6HDFVAL1 

♦6HDFVAL2 

♦6HDFVAL3 

B 

»6HDFLP 

♦6HDFLY 

♦6HDFVLH1 

,6HDFVLH2 

t 6HDFVLH3 

,6HGAMAD 

♦6HAZVAD 

C 

? 6HDMASS 

.6HHEATRT 

,6HDUA 

,6HDVA 

ff»HDWA 

.6HVELAD 

,6HDVWH1 

D 

,6HDVWH2 

♦6HDVWH3 

,6HDLRD 

♦6HTVLRD 

«6HATLD 

♦6HGLRD 

,6HVIDLD 

E 

♦6HDVX 

,6HDVY 

»6HDVZ 

1 6 HDX 

f 6HDY 

,6HDZ 

,6HE0 

F 

» fcHF 1 

,6HE2 

»6HE3 

.6HETAL 

,6HFAXB 

»6HFAYB 

♦6HFAZB 

G 

,6HFMXB 

»6HFMYB 

♦6HFMZB 

,6HFTXB 

,6HFTYB 

*6HF TZB 

.6HFVAL1 

H 

I 

♦6HFVAL2 

/ 

♦6HFVAL3 

♦6HGAMMAA 

♦6HGAMMAI 

»6HG AMM A R 

,6HGXI 

,6HGYI 


DATA MOT VC? 


1 

/6HGZI 

,6HH 

,6HHTBT 

*6HHTBTD 

♦6HHTLF 

,6HHTLFD 

f 6HHTTP 

2 

.6HHTTPD 

,6HHTRT 

,6HHTRTD 

,6HHTURB 

,6HHTURBD 

'.6HAZVELA 

,,6HAZVEL 

3 

♦6HAZVELR 

♦6HGCLAT 

♦6HGDLAT 

,6HL0NG 

t 6HL0NGI 

♦6HMACH 

,6HMASS 

4 

*6HP JETTS 

,6HPRES 

,6HPWD0T 

♦6HDYNP 

.6HOALPHA 

*6HTLHEAT 

» 6HGCRAD 

5 

» 6HDENS 

,6HRS0 

,6HRS 

,6HSIP1 

♦6HSIP2 

,6HSIP3 

,6HSIP4 

6 

,6HSIP5 

,6HS IP6 

»6HS IP7 

♦6HSIP8 

,4HSIP9 

» 6HS 1 P 10 

♦6HSIP11 

7 

♦6HSIP1? 

,6HSIP13 

»6HS IP 14 

t6HSIPl5 

,6H$IY1 

*6HS IY2 

,6HSIY3 

8 

,6HSIY4 

»6HS IY5 

»6HS IY6 

,6HSIY7 

»6HS IY8 

,6HSIY9 

.6HSIY10 

9 

t 6HSIY1 1 

.6HSIY12 

» 6HS IY1 3 

,6HSIY14 

*6HS IY15 

t6HTHRUST 

.6HTIME 

0 

*6HTMXB 

»6HTMYB 

♦6HTMZB 

♦ 6HTTMXB 

♦6HTTMYB 

,6HTTMZB 

♦ 6HTVAC 

A 

,6HU 

,6HV 

,6HW 

.6HUA 

,6HVA 

,6HWA 

♦6HVELA 

B 

,6HVAXI 

t 6HVAYI 

♦6HVAZI 

*6HUB 

»6HVB 

t6HWB 

♦6HVFLT 

C 

*6HDLR 

♦6HTVLR 

♦6HATL 

,6HGLR 

♦6HVIDEAL 

»6HUR 

» 6HVR 

D 

*6HWR 

♦6HVELR 

» 6HVRXI 

,6HVRYI 

,6HVRZI 

,6HUW 

,6HVW 

E 

♦ 6HWW 

.6HVWXI 

,6HVWYI 

,6HVWZI 

,6HVXI 

,6HVYI 

♦6HVZI 

F 

•6HVXL 

t6HVYL 

*6HVZL 

♦6HVXL0 

,6HVYL0 

♦ 6HVZL0 

,6HWD0T 

G 

,6HWEIGHT 

,6HWJETTM 

.6HWPR0P 

♦6HXCG 

♦6HYCG 

♦6HZCG 

,6HXI 

H 

T 

♦ 6HYI 

/ 

*6HZI 

*6HXL 

*6HYL 

,6HZL 

,6HXL0 

*6HYL0 

X 



DATA 

MOTVC3 




1 

/6HZL0 

,6HYAWRH 

♦ 6HDCLV 

♦6HDCDV 

♦6HVINV 

♦6HAE1 

♦ 6HAE2 
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2 

♦6HAE3 

,6HAE4 

,6HAE5 

,6H AF6 

,6HAE7 

,6HAE8 

t 6HAE9 

3 

♦6HAE10 

y6HAE 11 

,6HAE12 

,6HAE13 

,6HAE 14 

yfcHAEl 3 

y6HWDl 

4 

y 6HWD2 

,6HWD3 

,6HWD4 

»fcHWD5 

y6HWD6 

y6HWD7 

y6HWD8 

5 

y 6HWD9 

y6HWD10 

,6HHD1 1 

,6HWD12 

f 6HWD13 

»6HWD14 

y6HWD15 

6 

,6HTHR1 

y 6HTHR2 

,6HTHR3 

,6HTHR^ 

,6HTHP5 

♦ 6HTHR6 

,6HTHR7 

7 

y 6HTHP 8 

»6HTHR9 

,6HTHR 1 0 

.6HTHR11 

♦6HTHR12 

.6HTHR13 

,6HTHR14 

8 

y 6HTHR 1 5 

,6HISV 

,6HISV3 

,6HDIARP1 

.6HDIARP2 

,6HDIARP3 

,6HDRAGP1 


9 ,6HDRAGP2y6HDRAGP3y6HDRAGPT,6HFAXBPly6HFAXBP2y6HFAXBP3,6HCDPl 

0 * 6HCDP2 ,6HCDP3 ,6HD I AMP1 ,6HDI AMP2 y6HP IAMP3 y6HDRGPKl y 6HDRGPK2 

A ,6HDP&PK3,6HDRGPPl,6HDRGPP2t6KDRGPP3,6HVFLAP ,6HDRGPS ly 6HDRGPS2 
B *6HDRGPS3t6HIDRGPl,6HIDRGP2,6HIDRGP3,6HPARIFl,6HPARIF2t6HPARIF3 
C y 6HPL 10 y6HTVLID y6HGLID f 6HDLI y6HTVLI t 6HGLI y6HHTRTl I 

D ,6HHTRT1 yfeHTIME 0 ,6HTTIME ,6HAXG ,6HA YG ,6HAZG ,6HAH0RIZ 

E y 6HAVERT y6H IAER0H ,6HR0VET1 ,6HR0VET2y6HR0VET3,6HR0VET4,6HR0VET5 
F y 6HR0VET6 y6HR0VET7,6HR0VET8 ,6HR0VFT9,6HR OVETO 

1 / 


DATA PHZVC1 

1 /6HALTMAX t 6HALTM IN y6HMAXTIM » 6HEVTF ,6HFESN »6HIESN ,6HPHZF 

2 y6HPIF ,6HI4 ,6HSAVE£N 

3 / 


DATA SPECA1 

1 /6HSPECIl,6HSPECI2,6HSPECI3,6H5PECI4,6HSPECI5,6HSPECI6y6HSPECI7 

2 y 6HSPEC 18 ,6HSPEC I9,6HSPECV1 t 6HSPEC V2 f 6HSPFCV3 y6HSPE CV4 t 6HSPECV5 

3 ,6HSPECV6,6HSPECV7,6HSPECV8,6HSPECV9*6HNSPEC1,6HNSPEC2,6HNSPEC3 

4 y 6HNSPE C4y 6HNSPEC5 ,6HNSPEC6 

5 / 


DATA TG0VC1 

1 / 6HFUXN1 ,6HFUXN2 ,6HFUXN3 »6HFUXN4 ,6HFUXN5 ,6HFUXN6 .6HFUXN7 

2 y 6HFUXN8 f 6HFUXN9 ,6HFUXN10 y6HPCTG0 y6HSAVEl .6HSAVE2 t 6HSAVE3 

3 y 6HSAVE4 ,6HSAVE5 ,6HSAVE6 ,6HSAVE7 ,6HSAVE8 ,6HSAVE9 , 6HSAVE 10 

4 ,6HSAVE11 ,6HSAVE 12y6HSAVE13 f 6HSAVE 14 t 6HSAVEl 5y6HSAVE16 t 6HSAVE 17 

5 ,6HSAVE18,6HSAVE19,6HSAVE20y6HSAVE21 »6HSAVE22 y6HS AVE23y6HSAVE24 

6 y 6HS AVE25 y 6HSAVE26 y6HSAVE27 y6HSAVE28 y6HSAVE29,6HSAVE30y 6HSAVE31 

7 t 6H SAVE 3 2 t 6HSA VE33 y6HS AVE34 y 6H SAVE 35 ,6HS AVE36 y6HS AVE37,6HSAVE 38 

8 y 6HSAVE39y6HSAVE40y6HSAVE41 ,6HSAVE42 y6HSAVE^3y6HSAVE44,6HSAVE45 

9 t 6HS AVE46 y 6HSA VE47 t 6HS AVE48 t 6HSAVE49,6HSAVE50y6HS AVE5 1 » 6HSAVE 52 
0 *6HSAVE53y6HSAVE54,6HSAVE55,6HSAVE56,6HSAVE57y6HSAVE58y6HSAVE59 
A y6HSAVE60y6HSAVE61 t 6HSAVE62y6HSAVF63y6HSAVE64 f 6HSAVE65y6HSAVE66 
B y 6HSAVE67 y 6HSA VE68 y6HSAVE69 r 6HSAVE70 y6HS AVE71 ,6HS AVE72 y 6HS AVE73 
C »6HSAVE 74,6HSAVE75 y6HSAVE76,6HSAVE77y6HS AVE78y6HSAVE79,6HSAVE80 
D f 6HTG0 t 6HTIMX ,6HESN 

C y 6H1EVNT1 y6HIEVNT2 y6H IEVNT3 t 6HTF VNT4 1 6HI EVNT5 y6HI EVNT6 » 6HIEVNT7 
D y 6HIEVNT8 y6HIEVNT9 t 6HIENV10,6HISZEV .6HNXEVT ,6HI5 .6HGUXN1 
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E » 6HGUXN? *6HGUXN3 ,6HGUXN4 ,6HGUXN5 ,6HGUXN6 .6HGUXN7 ,6HGUXN8 
F *6HGUXN9 ,6HGUXN10,6HTIMY 
I / 


C 


c 


DATA TRACK1 

1 /6HCTKLT1 ,6HCTKLT2 ,6HCTKLT3 ,6HCTKLT4 ,6HCTKLT5 t 6HCTKLT6»6HCTKL T7 

2 f 6HCTKLT8,6HCTKLT9,6HGTKLT0,6HCTKLNl ,6HCKTLN2 ,6HCKTLN3 ,6HCKTLN4 


3 

»6HCTKLN5 

» 6HCTKLN6 

♦6HCTKLN7 

»6HCTKLN8 

i ,6HC TKLN9, 6HC TKLNO 

♦6HELEV1 

A 

» 6HELEV2 

t 6HELEV3 

*6HELE VA 

»6HFLEV5 

,6HFLFV6 

»6HELFV7 

»6HFLEVP 

5 

♦ 6HELEV9 

♦6HELEV10 

,6HLKA1 

,6HLKA2 

.6HLKA3 

♦6HLKAA 

♦6HLKA5 

6 

.6HLKA6 

,6HLKA7 

,6HLKA8 

,6HLKA9 

.6HLKA10 

,6HLKF1 

»6HLKB2 

7 

♦6HLKB3 

,6HLKPA 

,6HLKB5 

f 6HLKP6 

♦6HLKB7 

»6HIKE 8 

1 6 HLKP 9 

8 

t 6HLKP 1 0 

»6HPGT1 

» 6HPGT2 

,6HPGT3 

,6HPGTA 

♦6HPGTF 

,6HPGT6 

9 

t 6HPGT7 

,6HPGT8 

♦ 6HPGT 9 

♦6HPGT10 

,6HPGT11 

,6HPGT12 

♦6HPGT13 

0 

»6HPGT1 A 

1 6HPGT1 5 

♦6HPGT16 

,6HPGT17 

♦6HPGT18 

.6HPGT19 

,6HPGT20 

A 

» 6HPGT2 1 

*6HPGT2 2 

,6HPGT23 

»6HPGT2A 

,6HPGT25 

»6HPGT26 

♦6HPGT27 

B 

»6HPGT28 

»6HPGT29 

.6HPGT30 

,6HPGT31 

♦6HPGT32 

,6HPGT33 

,6HPGT3A 

C 

,6HPGT35 

»6HPGTj6 

.6HPGT37 

1 6HPCT3 8 

,6HPGT39 

,6HPGTA0 

,6HPGTA1 

D 

» 6HPGTA2 

»6HPGT43 

,6HPGTAA 

.6HPGTA5 

,6HRGT46 

,6HPGTA7 

»6HPGTA 8 

E 

♦6HPGTA9 

♦6HPGT50 

»6HPGT5 1 

,6HPGT52 

,6HPGT53 

*6HPGT5A 

,6HPGT55 

F 

t 6HPGT56 

» 6HPGT57 

♦6HPGTE8 

,6HPGT5« 

,6HPGT60 

♦6HPGT61 

»6HPGT62 

G 

1 6HPGT63 

,6HPGT6A 

»6HPGT65 

,6HPGT66 

♦6HPGT67 

.6HPGT68 

♦6HPGT69 

H 

I 

* 6HPGT70 
/ 

t 6HPGT71 

♦6HPGT72 

,6HPGT73 

,6HPGT7A 

♦6HPGT75 

,6HPGT76 




DATA 

TRACK2 




1 

/6HPGT77 

,6HPGT78 

.6HPGT79 

,6HPGT80 

♦6HPGT81 

*6HPGT 82 

,6HPGTS3 

2 

♦6HPGT8A 

♦6HPGT85 

♦6HPGT86 

,6HPGT87 

♦6HPGT88 

*6HPGT 89 

f 6HPGT90 


3 ,6HSLTRGl,6HSLTRG2,6HSLTRG3,6HSLTRG4 f 6HSLTR&5,6HSLTRG6,6HSLTRG7 
A ♦6HSLTRG8,6HSLTRG9 t 6HSLTRGO,6HSLOSU,6HSLOS12,6HSLOS13,6HSLOSlA 

5 »6H5L0S1 5 »6h'SL0S 16 »6HSL0S 17» 6HSL0S 18 »6HSL0S19»6HSL0S 10 »6HSL0S21 

6 » 6HSL0S22 » 6HSLOS23 »6HS L0S2A*6HSL0S25 t6HSL0S26 » 6HS L0S27 1 6HSL0S28 

7 ,6HSL0S29*6HSL0S20,6HSL0$31 ,6MSL 0S32 ,6HS L0S33,6HS L0S3A , 6HSL0S35 

8 t 6HSL0S36 * 6HSL0S37»6HSL0S38 t 6HFL0S39 »6HSLOS30 *6HS TKLT1 » 6HSTKLT2 

9 ,6HSTKLT3,6HSTKLTA,6HSTKLT5 t 6HSTKLT6,6HSTKLT7,6HSTKLT8,6H£TKLT4 

0 1 6HSTKLT0 »6HS TKLN1 » 6HSTKLN2 *6HSTKLN3 »6HSTKLNA »6HSTKLN5 »6HSTKLN6 
A ,6HSTKLN7 t 6HSTKLN8,6HSTKLN9,6HSTKLN0*6HTKLTCl,6HTKLTC2,6HTKLTC3 
P ,6HTKLTCA,6HTKLTC5,6HTKLTG6,6HTKLTC7,6MTKLTC8 t 6HTKLTC9,6HTKLTC0 
C ,6HTKRA01 t 6HTKRAD2,6HTKRAD3,6HTKRADA t 6HTKRAD5,6HTKRAD6,6HTKRAD7 
D 1 6 HTKR ADS »6HTK*> AD9 f 6HTKR ADO »6HTK AZMI » 6HTKAZM2 *6HTKAZM3 *6HTKAZM4 
E ,6HTKAZM5*6HTKAZM6,6HTKAZM7,6HTKAZM8,6HTKAZM9,6HTKAZM0,6HTRKLT1 
F ,6HTRKLT2 t 6HTRKLT3,6HTRKLTA t 6HTRKLT5,6HTRKLT6,6HTRKLT7,6HTRKLT8 
G ,6HTRKLT9,6HTRKLT0,6HTkKHTl ,6HTRKHT2 ,6HTRKHT3 ,6HTRKHT4,6HTRKHT5 
H ,6HTRKHT6,6HTRKHT7,6HTRKHT8,6HTRKHT9,6HTRKHT0*6HTRKLN1,6HTRKIN? 

1 / 

DATA TRACK3 

1 /6HTRKLN3 ,6HTRKLNA ,6HTRKLN5 ,6HTRKLN6 »6HTRKLN7 » 6HTRK LN8 ,6HTRKLN9 

2 ,6HTRKLN0,6HTKRX1 *6HTKRX2 ,6HTKRX3 »6HTKRX4 .6HTKRX5 ,6HTKRXfc 

3 ,6HTKRX7 ,6HTKRX8 ,6HTKRX9 ,6HTKRX10 ♦6HTKRX1 1 ,6HTKRX12 ,6HTKRX 12 
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4 ,6HTKRXl?,6HTKRX14,6HTKPXlf> , 6HTKRX 1 6 ,6HTKF XI 7,fcHTKR XI 8 ,6HTKRX 19 

5 ,6«TKRX20t6HTKRX?1 ,6HTKRX22 ,6HTKRX2? ,6HTKRX24,6HTKR X25 ,6HTKRX26 

6 ,6HTKRX27,6MTKRX28 ,fcHTKRX29 t 6HTKRX?0 ? 6HTPKFL 1 »6HTRKFL2 t 6RTPKFL? 

7 »6HTRKFL^» 6HTRKF L6 »6HTRKFL7» 6HTP KF L8 »6HTFKFL9 »6HTRKFL0 »feHTR KNM1 
S « 6HTR KNM2 t 6HTRKNM3 ,6HTPKNM4 , 6HTP KNM5 »6HTRKNM6 »6HTRKNM7»6HTRKNM8 
9 T 6RTP KN'M9 ,6HTP KNMO t 6HNTR KS *6HNPUM »6HF LEMIN 

I / 

C 


C. 

c 

c 


TATA GUIPI1 

1/6HDTC ,fcF'CVRIl 1 6 H G VP I 2 ,6HGV*I3 ,6HGVFI4 ,6HGVRIL y6RGVR 16 
2 »6HGVR 17 ,6RGVH6 ,6KGVRI9 , 6HG VR 1 1 0, 6HIGF 1 »6HIGF? , 6HIGF3 
2 *6MI GF4 » tH I GF c . ,6HIGF6 
4/ 


LATA GUinVl 

1/6HGT1MC ,<HPGTlME,dHGPXl ,6HGPY] ,6HGPZI »6HGPVX I ,6HCPVYI 

2 »6HGP VZ I 1 6 P. C P A X I » £ HGF AY I ,FHGPAZI »6HG ALPHA yGHCEFTA » 6 HG F A NK 

3 »6HGYAWR ,CRGPITR ydHCPOLR ,6HGRPLI ,6R-GYAWI * 6H G P I T I ,6RGA5M 

4 t 6HGPA5yi,6HGPASYI,6HGPASZI t 6HGTHPST t fcHGWGT ,6HGWD0T ,6HGVRC1 
F »6RGVP C? ft RGVR C ? ,6HGVEC4 ,6HGVFCF ,6HGVRC6 t dHGVPC7 ,6HGVRC£ 
6,6HGVPC9 , dHGVRCl 0 

7/ 


DATA TAPGV1 

1/6HAL T AT ,6HALTPT »6H AN'CMOT » 6HA PORT * 6HAR GPT »6HAF GVT »6HDE CLT 
2 » 6HE CCENT »6H F NR GYT »6HHYPVT ,6HINCT ,6HLANT ,6HPFR ] DT,6HPC-FRT 
3»6HPGLPNT , 6HPGC LTT , 6HRTASCT » 6HS FMAXT, 6HTRRANT ,6HTRl!MXT,6HE TP 1 
4, 6 HE TP 2 ,6HETR3 » 6HE TR4 , 6HF TP 3 ,6HETR6 ,fcHETP7 , 6HE T & 8 

5 »6HE TR9 ,6HDXRT »6HPYR T ,6HnZPT ,6HPVXRT ♦6HDVYRT ,6H0VZRT 

6,6HPXI »d RPYI »6HDZI »6HDVXI ,6HPVYT ,6HDVZI *6HDAXI 

7 t 6HDAYl »6HDAZ I t 6HXIT »6HYIT y6MZIT t 6RVXIT ,6HVYIT 

8 tf-HVZ IT » 6HAXI T « 6HA YI T » 6HA Z I T , 6HGCP APT *6HVELI T ,6HTDXR 

9*6HTGYP v 6hTDZR ,6HTDVXR y 6HTPVYR ,6WTPVZP ,6 HDVCI RT t 6HECCANT 
0,6HMEAANT,6hTJM5PT , 6-HT1 MTPT , 6HPC VE LT ,6HAP VE LT , 6HGT AX I »6HGT AYI 
A »6HGTAZ I »6HCTVX I ,6HGTVYI f 6HGTVZI ,6*->GTX I ,6HGTYI ,6HGTZI 

F ,6HHVEHF 1 t 6hMVEHF 2 , 6HHVEHF 3 » 6HM VEHF 4 f 6HMVE HF3 *6HMVEH F6 ,6HHVFHF7 
C 1 6 HMVEHF 8 1 6KMVEHF 9 T 6HMVEHF0 1 6 HA Z VIT ,6HGAKIT ,6HIGT1 ,6HIGT2 
C,6HIGT3 »6HIGT4 ,6HIGT5 »6H1CT6 ,6MIGT7 ,6HIGT6 ,6HIGT9 

E, 6HDRT t6HVCIRCT,6HLANVET,6HGCLATT,6HL0NGT ,6HXLTC » 6HYLTP 

F, 6H?LTC ,6HVXLT0 ,6HVYLTC ytHVZLTP ,6HXLT ,6HYLT ,6HZLT 

GydHVXLT ,6HVYLT ,6HVZLT »6HAXLT , 6HAYLT » 6H A Z L T »6HPXPT 

H,6HPYPT 1 6 HDZ P T ,6HDVXPT t 6HPVYPT ,6HPVZPT t 6HDAXP T ,6HPAYPT 
1/ 

DATA TARGV2 

1/6HPAZ F T * 6HLVVX I ,6HDVVYI ,6HDVV2I »6HPVVXIT »6HPVVYIT »6HLVVZ1 T 
1/ 

END 
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ORIGINAL PAGE IS 
OF POOR QUALITY 



DINPT: This routine reads the previously stored input data 

from the disc and locates the data for the current phase. 





YES 


STORE THE DATA 
FOR THIS EVENT IN 
THE APPROPRIATE 
CORE LOCATIONS 
IN COMMON IV 
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DL00K (NAME) : This routine sets the core addresses of the 

variables (NAME) to be integrated into the integration list. 
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DPRNG: This routine calculates the range based on the dot 

product of the initial position vector of the vehicle and the 
Impact point vector with the oblate planet. 
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DR0P: This subroutine drops tight constraints 


NORMAL 

LOOP 



TAKE ALL COMBINATIONS 
OF ACTIVE CONSTRAINTS 
TO DETERMINE IF THERE 
ARE ANY COMBINATIONS 
FOR WHICH 

NAC .LE. NINDU 
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DTMDL: This routine checks the user-specified tables to 

ensure that the next integration step size is less than or equal 
to the next time point in any of the tables. 



III-69 




DVADDM: This routine adds instantaneous delta inertial veloc 

at the beginning of any phase. 


















DYNS1: This routine integrates the equations of motion using 

the Runge-Kutta method. 
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DYNS2; This routine integrates the equations of motion using 
the predictor-corrector method. 



RESTART INTEGRATION 
IF DELT WAS CHANGED 
IN CYCXM 


RESTART INTEGRATION 
IF A DISCONTINUITY 
OCCURRED 







DYUS3: This routine integrates the equations of motion using 

Laplace’s Method of integrating orbits about a spherical planet. 











DYNXA: 

integration 


This routine initializes the fourth-order Runge-Kutta 
scheme . 



SET T1IE INTEGRATION PASS 
AND DISCONTINUITY 
FLAGS TO ZERO 
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DYNXB: This routine initializes the fourth-order predictor- 

corrector integration scheme. 



! ' SET THE INITIALIZATION FLAG 
AND THE HISTORY VALUE OF 
THE INTEGRATION INTERVAL 







G- ^ 
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DYNXMI: This routine selects the integration scheme to be 

used. 



NOTE: 

DYNXA - FOURTH-ORDER RUNGE-KUTTA INTEGRATION 

- LAPLACE’S METHOD OF INTEGRATION 

- ENCKE'S METHOD OF INTEGRATION 

DYNXB - FOURTH-ORDER PREDICTOR-CORRECTOR INTEGRATION 
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DYSI1: This routine sets the integration list. 



^NPC(24) : 0 J 1 

■i* 


TURN ON THE 
GENERAL 
INTEGRATION 
VARIABLES 



^ i 

3 


d: 



i 

ez_ 


ZERO OUT THE CORE 
ADDRESSES OP THE 
INTEGRALS AND 
DERIVATIVES 



SET THE NUMBER 
OF INTEGRALS TO 
ZERO 


< 


DL00K (DYNIL) 



> 


CALL DL00K TO SET UP 
THE LIST OF INTEGRALS 
FOR THIS PHASE 


NOTE: 

N = THE NUMBER OF INTEGRALS FOR 
THE CURRENT PHASE 
LISTI = INTEGRALS 
LI STD = DERIVATIVES 
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' FGAMA(IS): This subroutine calculates the values of PI and 

P2 associated with a particular GAMA in the direction of search 
by changing the controls according to the equation 


U(I) = U(I) + GAMA * DU (I) 
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F0PC: This routine contains the fourth-order predictor- 

corrector integration algorithm. 
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FXRNG2(T): This routine calculates the crossrange parameter 
used in the iteration scheme required to determine the crossrange 
relative to a reference ground track. 



COMPUTE CENT 
AND 0MEGA 


l 

COMPUTE L01IGR 
AND SIGT 


COMPUTE THE COSINE: 
OF THE RANGE 

T ; 

C0?1PUTE RANGE 
ANGLE TO THE 
VEHICLE 


i 


COMPUTE SIGV 
AND SIGTV 


i 


COMPUTE THE 
CROSS RANGE 
PARAMETER 
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GABDD: This routine computes the maximum step size in the 

direction of search, STPMAX. STPMAX is computed as the minimum 
step needed to make one of the inactive inequality constraints 
become active. This prediction is based on a linearization of 
the inactive constraints. 



ESTIMATE STEP SIZE 


TO INACTIVE 
INEQUALITY 
CONSTRAINT 
BOUNDARIES 


STPIIAX IS TAKEN 
EQUAL TO THE MINIMUM 
OF ALL OF THESE STEP- 
SIZE PARAMETERS 
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Mili: ri 


GAMLAH: This routine calculates the flight path angle and 

azimuth angle of the atmospheric relative velocity vector* 
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GCNTRL: This is a blank routine used for simulating the 

hardware lags and errors associated with implementing the closed- 
loop guidance steering commands. 
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GENMIN; This routine finds the minimum of the function y, 
given the pairs (X, Y) and the slope at X * 0. If the pairs (X, Y) 
are not given, then GENMIN generates these pairs automatically. 

The minimization is based on approximations to the functions 
that are made using quadratic and cubic polynomials* The analyti- 
cally calculated minimums of these polynomials are used as esti- 
mates to the actual minimum values of y. 
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MINIMIZE Y VIA 
3 POINTS AND NO SLOPE 







GGUID: This is a blank routine that is to be used for simu- 

lating the closed-loop guidance equations being analyzed. 
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GMAG: This routir 

vectors G1(I) and G2(I) 














GNAV: This is a blank routine that is to be used for simu- 

lating the navigation equations for the closed-loop guidance system 
being analyzed. 
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GRAD: This routine computes the gradients of the penalty 

functions (G1 and G2) with respect to the control parameters. 
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GRADS: This routine computes the gradients of the steering 

constraints with respect to the steering control parameters. 

i 





CALES 

i 

COMPUTE STEERING SENSITIVITY 
MATRIX AC0BS (I) 











GRAV: This routine calculates the gravity acceleration vector. 
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GSA: This routine (Generalized Steering Algorithm) computes the necessary 

values of steering control variables to yield the desired values of the steering 
constraints . 




SET NUMBER OF ITERATIONS = 

EH 


20 i 



INCREMENT NUMBER OF 
ITERATIONS BY 1 



Gb 


CALCULATE NEW STEERING CONSTRAINT 
VALUES BASED ON CURRENT VALUE OF 
STEERING CONTROLS 


J 



SAVE NOMINAL ERRORS 
EN0MS = ES 



i 


COMPUTE WEIGHTED ERROR 
WES = S0RT (SUM ES(l) ** 2) 







GSENSR: This routine is to be used for simulating the inter- 

face between the trajectory simulator (real world) and the guidance 
sensor (hardware detected) being analyzed. 
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GUIDI: This routine initializes the guidance (steering) par- 

ameters based on the option selected by the user via the IGUID array. 



INITIALIZE THE 
HISTORY VALUES 
FOR LINEAR 
FEEDBACK 
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GUIDX (N) : This function returns as the commanded angle for 

channel N, where N=l,2,3, based on the generalized linear feedback 
guidance (steering) algorithm. 
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GUIDl: This routine calculates the steering angles based on 

the user-specified steering option. 




COMPUTE AERODYNAMIC 
ANGLES BASED ON 
SELECTED STEERING 
OPTION 


COMPUTE EULER 
ANGLES FOR 
EITHER RELATIVE 
OR INERTIAL ANGLES 

j ~ 
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Ill II Rl 


GUID2: This routine calculates the Euler paramenter deriva- 

tives based on the inertial roll, pitch, and yaw, or the angle-of- 
attack, sideslip, and bank angle rate commands. 









ROLL, ALPHA, AND 
BETA RATE COMMANDS 
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IISWGT: This routine calculates the total weight of the heat 

shield by summing the weights of its individual components. 
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IBMTRX: This routine calculates the inertial-to-body (IB) 

matrix, based on the guidance (steering) option selected by user 
input . 



COMPUTE THE I I COMPUTE TIIE I I COMPUTE THE 
INERTIAL I I COMPUTE THE INERTIAL- TO- I I GEOGRAPHIC- 
LAUNCH-TO-BODa I INERTIAL-TO- BODY (LB) I I TO-BODY 


(LB) MATRIX I ATMOSPHERIC I MATRIX USING I I (GB) MATRIX 
USING THE I CIA) MATRIX I R0LI, Y AWI, I I USING YAWR, 
QUATERNIONS I I PITI I I PITR, R0LR 
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INFXM: 
functions . 
calculation 


This routine performs the output data processing 
It also calls subroutine C0NIC, depending on the conic 
option requested by NPC(l). 



CONVERT TO METRIC 
OR ENGLISH UNITS 
IF REQUIRED 


I 
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( TIME TO WRITE 
THE PROFILE TAPE? 

I yes 


UPDATE THE 
PROFILE WRITE 
TIME 


WRITE THE 
REQUESTED 
PARAMETERS ON 
THE PROFILE TAPE 


RECONVERT 
TO ENGLISH 
OR METRIC UNITS 
IF REQUIRED 



















-4 



INPUTX : This routine prints a summary of the input table 

data ' The variable LISTIN' is also checked to determine whether or 
not to rewind the output file, which eliminates the listing of the 
input data deck. 
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INTGRL (LIST, NUM, KEY): This routine initializes the list 
of variables to be integrated. The list of variables to be in- 
tegrated is called DYNIL . It contains 148 cells that are stored 
in subroutine BLKDAT. There are three cells for each integrated 
variable, which corresponds to a total of 49 integration variables. 
The first cell indicates the total size of the array, the second 
cell contains the Hollerith name of the first variable to be 
integrated, the third cell contains the Hollerith name of its 
derivative, and the fourth cell contains a flag to indicate whether 
or not to integrate that variable, etc. If the flag is zero, the 
variable is not integrated; if the flag is equal to 1, the vari- 
able is integrated. 


INTGRL is also used to turn the integration of variables on 
or off as desired. For example, if NPC(ll) = 1, we wish to acti- 
vate the inequality constraint integrations (L.e. , the integration 
o FVAL1, FVAL2, FVAL3) . In this case, subroutine M0TIAL calls 
INTGRL as follows to activate the integration: 

CALL INTGRL (DYNIL (38), N03, KOI) 1 — - 


where : 


DYNIL (38) = the position of FVAL1 in the array DYNIL 

W03 = fixed point 3, which means that the three 

integrals namely, FVAL1, FVAL2, and FVAL3 , 
are to be turned on, since they are in 
sequence. ' 

= fixed point 1; this means turn on the 
integration of the variables. If the 
argument were N00 (fixed-point zero), 
the integration would be turned off. 
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IRTBR : This routine calculates body rates from inertial 

Euler angle rates. 
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ITER0: Main program of overlay (2,6). This routine prints 

out the iteration summary as required. The iteration summary con- 
tains all the information relating to the search/optimization 
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MASTER: This is the main program of overlay (0,0). It 

decides whether to read input data or execute the problem. 
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MINMYS: This routine contains the optimization executive 

logic for the various options that are available. It phases the 
various iteration paths by examining various properties of the 
objective function and the constraint manifolds. 
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M0MENT: This routine calculates the total thrust and 

aerodynamic moments for static trim. 
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M0NITR: This routine determines the maximum and minimum 

values of the user-specified monitor variables. 



FIND THE 
ADDRESS OF 
THE MONITOR 
VARIABLES 

T 

COMPARE THE 
CURRENT VALUE 
WITH THE MAXIMUM 
AND MINIMUM VALUES 


UPDATE THE 
MAXIMUM AND 
MINIMUM VALUES 
AS REQUIRED 






COMPUTE YAW 
MOMENTS FOR 
THE TRIMMING 
ENGINES^ __ 

T 

COMPUTE YAW 
AERODYNAMIC 
MOMENTS 

ZEZ 

COMPUTE YAW 
MOMENT TO BE 
CANCELLED BY 
THE TRIMMING 
ENGINES 



COMPUTE THE 
TOTAL PITCH 
AND YAW MOMENTS 




M0TENC: This routine computes the difference between the 

total acceleration on the vehicle and the two-body accelera- 



GO TO: 

100, 700, 300, 
700, 400, 100; Ll 


^ 00 ^ 


CALCULATE KEPLERIAN 
POSITION XL(I) 

AND VELOCITY VXL(I) 






(£o) 


CALCULATE XL(I) AND 
VXL(I) AT 

TIME = TIME +0.5 DELT 
BY CALLING TW0BDY 



CALCULATE TWO-BODY 
ACCELERATIONS AXL(I) 
AT XL(I) 




COMPUTE DIFFERENCE IN 
ACCELERATION AS 
DAX(I) = AXI(I) - AXL(I) 




UPDATE XL 
VXL(I) AT 
MIDPOINT 
INTEGRATI 

(I) AND 
THE 

OF THE 
ON INTERVAL 




< 


CALL H0TION 
TO COMPUTE AXI(I) 


> 


COMPUTES 
— i TOTAL 

ACCELERATION 


^NPASSj^l^--«» ^0^ 


CALCULATE NEW 
XL0(I) AND VXL0(I) 
AND SET DX(I) AND 
DVXI(I) BACK TO ZERO 


CALCULATE TWO-BODY 
ACCELERATIONS , 
AXL(I) 


T 

(1200) 


in-120 




I'j " 1 ! j !• y • i; W| jl l] III ^ 
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M0TIAL: This routine initializes the equations of motion. 





IS THIS THE FIRST 
INTEGRATI ON PAS S ? 

YES 


DETERMINE THE 
INITIALIZATION 
OPTION AND 
COMPUTE THE 
INITIAL POSITION 
AND VELOCITY 
PARAMETERS 




INITIALIZE THE 
IL MATRIX THRUST 
INCIDENCE ANGLES 
AND THE STATIC 
TRIM PARAMETERS 


COMPUTE THE 
IG MATRIX 


COMPUTE THE 
ATMOSPHERIC WINDS 


INITIALIZE THE 
AUTOPILOT 


r\ 


CHECK FOR INERTIA 
DISCONTINUITIES 


I 


INITIALIZE THE 
VEHICLE WEIGHT 
PARAMETERS AND 
ADD DELTA VELOCITY] 
IF DESIRED 


INITIALIZE THE 

GRAVITY 

PARAMETERS 


TURN ON THE 
INTEGRATION OF 
THE OPTIONAL 
INTEGRATION 
VARIABLES FOR 
THIS PHASE 
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M0TI0N: This routine calculates the 


equations of motion. 







COMPUTE HOLDDOWN 
EQUATIONS 


RETURN 


COMPUTE TOTAL - 

ACCELERATION 

COMPONENTS 


COMPUTE 

AEROHEATING 

DERIVATIVES 
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N0MH0L : This routine initializes the values of all HOLLERITH 

variables to the stored values. 
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N0MINL: This subroutine runs the nominal trajectories (one 

per iteration), saving core at the beginning of each phase. The 
routine also calls TEST to determine if the iteration reference 
has converged or failed to converge. 



WRITE: 

*** ITERA r 

riON LIMIT 




C maxitr • ° M 


9 


P2N0M ; 

c P2MIN 

_J 

\ 




rz 

I> 




-o 


WRITE: 

*** TIME LIMIT 
ABOUT TO BE EXCEEDED 


SOLVED = 77B 


6 





6 
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o 


1 

I NFLAG t 0 M 


IPRT = 77B 


C nsth ' 11 1 x 


I 


REWIND 2 


< H ^ 


REWIND 2 


CALCULATE 
P1N0M, P2N0M, 
AND EN0M(I) 


CALCULATE THE 
NUMBER OF TIGHT 
CONSTRAINTS, NTC 
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0LGM: This is a blank routine that is to be used for special 

open-loop guidance (steering) models. The polynominal coefficients 
or angular values to be used by GUID1 can be calculated in this 
routine and then used by the user— selected option, based on the 
IGUID array. 



T 
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0LG0M: This routine allows the guidance (steering) option 

values determined by the IGUID array to be overridden if a 
specified parameter test has been violated by the commanded 
angle. For example, the commanded angle of attack can be 
overridden if the value of QALPHA exceeds an input limit. This 
allows the program to follow the limit until it is no longer 
violated by the commanded angle . 
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4 PAGER(N): This routine determines when a new page is required 

prior to printing. The argument is the number of lines to be 
printed . 



PRECEDING PAGE BLANK NOT FILMED' 
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PBL0CK (N) : This routine generates a summary print block 

for the option requested by the argument N. 



PRINT PRINT PRINT PRINT 

CONIC TRACKER VELOCITY VACUUM 

BLOCK BLOCK LOSSES IMPACT 

BLOCK BLOCK 
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PGM: This routine determines the direction of search when 

the projected gradient method is used. 





c 


IF: 


0PT = 0 
INTRY1 + 0 
INTRBL 4 0 




— © 



THE NUMBER OF ACTIVEX 
CONSTRAINTS LT THE NUMBER] 
INDEPENDENT VARIABLES# 


COMPUTE LEAST-SQUARES 
DIRECTION OF SEARCH 


©I 


1 COMPUTE 1 

1 G2(I) I 

| 


IF THE 
ANY AC 
CONSTR 
COMPUT 

RE ARE 
TIVE 
AINTS 
E P2 


< 


IF THE CONSTRAINTS 
ARE NOT SATISFIED 


1 

NO 

COMPUTE TH 
THE SINGLE 
G1(I) and 
GIVEN BY C 

E COSINE OF 
BETWEEN 
PG1 (I) 

THA 



MAKE DAVIDO 
TESTS ON THI 
ISTART, NACE 

N RESTART 
l PARAMETERS 
>, AND NAC 


■^YES^^ 

ME) 




COMPUTE 

THE I 



PROJECTED | 



GRADIENT, PG1(I)I 



|g| 




COMPUTE 
RETARGETING 
DIRECTION 
OF SEARCH 
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PHZXM: This routine is the executive routine for overlay 

(2,3). It controls the integration of the equations of motion 
and determines whether the parameters MAXTIM, ALTMIN, or ALTMAX 
have been exceeded; if they have, the trajectory is terminated. 
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PHZXMI: This routine is the executive program of overlay 

(2,2). It performs the executive function for the phase initial- 
ization process. 
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i a 


PRNTIC: This routine prints a summary of the initial con- 

ditions for the current phase. 



{ CHECK FOR THE FIRST 
PASS; IF NOT, BYPASS 
SEARCH/OPTIMIZATION 
INPUT 
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PR0P: This routine calculates the thrust forces in the body- 

coordinate system, based on the use of either rocket or jet engines. 
The routine also determines the value of the throttling parameter 
required to limit the acceleration when NPC(7) =1. 
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QMULT: This routine multiplies quaternion El and quaternion 

E2 to yield quaternion E3. 
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QRTATE: This routine computes the quaternion E that results 

from rotating through the specified angles, ANGLE, about the speci- 
fied axes, NAXIS. 



COMPUTE TEMP (J) , WHICH 
IS THE QUATERNION FOR 
ANGLE (I) , AND STEMP(J), 
WHICH IS THE BACK VALUE 


< QMULT 

(STEMP, TEMP, E) 



NO 


< IS THIS THE"\ 
LAST ANGLE ? 3 


YES 




QUAT1: This routine computes the quaternions, given the 

angle of attack, sideslip, and bank angle. 



TEMP =90.0 
RANGLE (1) = 
RANGLE(2) = 
RANGLE (3) = 
RANGLE (4) = 
RANGLE (5) = 
RANGLE (6) = 
RANGLE (7) = 
RANGLE (8) = 
RANGLE (9) « 


AZL 

LATL 

L0NGI 

-LATC 

LAMDA 

GAMA 

SIGMA 

-BETA 

ALPHA 


- L0NL 

- TEMP 



SET UP THE 
ANGLES OF 
ROTATION 




COMPUTE THE 
QUATERNIONS 
FOR THE GIVEN 
ANGLES 
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QUAT2: This routine computes the quaternions, given the 

local attitude angles in yaw, pitch, and roll. 



SET UP THE 
ANGLES OF 
ROTATION 


COMPUTE THE 
QUATERNIONS 
FOR THE GIVEN 
ANGLES 
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QUAT3: This routine computes the quaternions given the relative 

euler angles. 



^ GAMLAM \ 


RANGLE (1) 
RANGLE (2) 

RANGLE (3) 

RANGLE (4) 
RANGLE (5) 
RANGLE (6) 
RANGLE (7) 


AZL 

LATL 

L0NGI - L0NL 

-GCLAT - 90 

YAWR 

PITR 

R0LR 


< 


I 


QRTATE (RANGLE ,NAXIS ,E0) 



>- 


SET UP THE 
ANGLES OF 
ROTATION 


COMPUTE THE 
QUATERNIONS 
FOR THE 
GIVEN ANGLES 


i 

■ 

a 
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READAT : This routine performs the actual processing of the 

input data. Subroutines RSEARCH, RGENDAT, RTBLMLT, and RTAB are 
called as required to perform the actual reading of namelists 
"SEARCH", "GENDAT", "TBLMLT" , and "TAB", respectively. The name- 
lists are always read in a given sequence that can be terminated 
at any point by setting ENDPHS = 1. The calling sequence for 
reading the namelists is: RSEARCH, RGENDAT, RTBLMLT, and RTAB. 

After reading each namelist, the data for that phase are 
packed into one of two data buffers, depending on the type of 
data being processed. The two data buffers are broken down as 
follows: f 

1) IGEN : 1500 decimal cells of storage for all input vari- 

bles except for event criteria and input tables for all 
phases. The table multipliers are also stored in this 
array. 

2) I BKT : 1500 decimal cells of storage for the event cri- 

teria and tables for all phases. 

The actual storage detail is shown in Table 2, The input 
variables are stored in sequence as they are read and the data 
are grouped according to phase number. 
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TABLE 2 


CONTENTS OF DATA BUFFERS 


Location 

BUFFER "IGEN" 
(1500 decimal cells) 

Location 

BUFFER "IBKT" 
(1500 decimal cells) 

i 

number of cells occupied in 
IGEN 

i 

number of cells occupied in 
IBKT 

2 

number of cells occupied by 
data for the first phase 
(NG1) 

2 

number of cells occupied by 
event sequence data (NB1) 

3 

event sequence number for the 
first phase 

3 

event sequence number for 
the first phase 

4 

address of the first varia- 
ble following the event 
sequence number in IGEN 

4 

event type for the first 
phase 

5 

value of the first varia- 
ble in IGEN 

5 

event criteria for the first 
phase 

6 

unused cell associated with 
the first variable stored 
in IGEN 

6 

criteria value for the first 
phase 

7 thru NG1 

variables stored by repeat- 
ing the sequence shown in 
4, 5, 6 for each variable 

7 

derivative name of the 
event criteria variable 
for the first phase 



8 

tolerance 

NGl+l 

number of cells occupied 
by data for the second 
phase (NG2) 

9 

model number 

NGH-2 

repeat sequence 3 through 
7-NG2 as before for the 
second event 

10 

unused cell associated with 
the event criteria for the 
first event 



11 thru NB1 

repeat sequence 3-10 for 
remaining events 
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TABLE 2 


CONTENTS OF DATA BUFFERS - CONCLUDED 


Loc a Lion 

BUFFER "IGEN" 

Location 

BUFFER "IBKT" 



NBl + 1 

number of cells occupied by 
all tables for the first 
phase (NB2) 



NB1 + 2 

phase number associated with 
the first set of tables 



NBl + 3 

size of the first table to 
be input 



NBl + 4 

name of the first table 
(HOLLERITH) 



NBl + 5 

table pointers and values 




repeat above sequence for 
each table in the first 
phase (NBl -f 3 thru NBl + 5) 




repeat above sequence for 
all phases 

(NBl -f 1 thru NBl + 5) 




* . 
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III-148 








DETERMINE TABLE TYPE 

- CONSTANT 

- MONOVARIANT 

- B I VARIANT 

- TRIVARIANT j 

r 

CONVERT INDICES AND 
POINTERS TO INTEGER 
REPRESENTATION. 

CHECK TABLE DIMENSIONS. 
INSERT INTO IBKT PACKED 
TABLE DATA ARRAY 
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PRINT TABLE INPUT 
SUMMARY FOR CURRENT 
CASE 


RESTORE DICTIONARY 
FROM DISC - UNIT 1 


i 


! 




rC 


IS THIS A 
SEARCH RUN? 

YES 




ISET ADDRESSES FO: 
■SEARCH VARIABLES 


SET RELATIVE ADDRESSES FOR 
INTEGRALS - CALL SREL - 


FIND RELATIVE ADDRESS 
FOR PRINT REQUESTS 

I 

FIND RELATIVE 
ADDRESSES OF TABLE 
ARGUMENTS 

FIND ADDRESS OF 
TABLES RELATIVE TO 
IBKT AND STORE IN IGEN 

t ~ 

SET RELATIVE 
ADDRESS OF 
EVENT CRITERIA 




FIND ADDRESS OF 
EVENT CRITERIA 
USED AS CONTROLS 
RELATIVE TO IV 


WRITE DAI 
RELATIVE 
DISC (UNI 

’A ARRAYS 
TO IV ONTO 
T 1) 




C 


WERE THERE FATAL 
SEARCH ERRORS? 


C 


I 


NO 


CALL ERROR I 

■x^YES I J 


IS THIS THE LAST 
CASE, ENDJ0B f 0? 



rv YES f I 

)? V»-llNPCF = 7 7 b| 
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REVISE: This routine determines the indices of active con- 
straints that are candidates to be dropped. This is done by ex- 
panding Gl(I) in terms of the elements of the sensitivity matrix 
The constraint with the most negative coefficient DFDC(I) is a 


candidate to be dropped. 
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RGENDA : This routine reads namelist "GENDAT" and checks 

for any namelist errors. 
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RSC0RE: This routine checks NXTRUN, IFRUNF(i), and ISC0RE 

to determine if criteria for executing the next case are met. 

A history of success or failure of previous cases are packed into 
ISC0RE by bit position. IRUNF is returned nonzero if the next 
case is to be executed. 
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and checks 


RSEARC : This routine reads namelist "SEARCH” 

for any namelist errors. 



I 


1 

1 

r 

1 ENDPRB 

= 0.0 | 

j 

i 


READ NAMELIST (SEARCH) 


] 


YES 


- ^INPCF < 0? ^ 1 


NO 


CHECK TO SEE IF 
AN INPUT ERROR 
HAS PREVIOUSLY 
OCCURRED 


C 


NMLTER ( NAMER ) 


r 


3 


NO 


1 


RETURN IF NO 
INPUT ERROR HAS 
OCCURRED 


INPCF - -77B 


I 


1 


{ SET THE INPUT ERROR 
FLAG IF AN ERROR HAS 
OCCURRED 


1 


READ THE REMAINING 
INPUT NAMELISTS IF 
AN ERROR HAS OCCURRED 
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RTBLML: This routine reads namelist "TBLMLT" and checks 

for any namelist errors. 
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"TAB" and checks for any 


{ CHECK TO SEE IF AN 
INPUT ERROR HAS 
PREVIOUSLY OCCURRED 



RSEARCH 




RUK: This routine contains the Runge-Kutta integration 

algorithm. 
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SAVIC: This routine buffers out commons IV and IBKT, which 

contain the state conditions, to file 2 at the beginning of each 
phase that contains an independent control variable. This infor- 
mation is used in running perturbed trajectories in the discrete- 
parameter mode. 
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SC0RE: This routine sets the ICASE bit of ISC0RE(1) or 

ISC0RE(2) if S0LVED is + or respectively. No bit is set if 
S0LVED = 0. 



SET ICASE BIT OF ISC0RE (2) 

1 ' I 


| SET ICASE BIT OF ISC0RE (1) 
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SETIC (NPHASE) : This routine resets the initial conditions 

for the specified phase (NPHASE) equal to their nominal values. 
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SUN: This routine calculates the right ascension and declination 

of the sun for any Julian date (from 1900 January). 
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TEST: This routine tests for convergence. 




IS THE PROJECTED 
GRADIENT METHOD 
BEING USED? 

|ncT"*^ 




IS 0PTVAR 
OPTIM IZED ? 

YES 







COMPARE U, PI, P2, and G2; 
WITH THEIR VALDES ON THE 
PREVIOUS ITERATION TO 
DETERMINE IF THE ITERATION 
SHOULD CONTINUE 



IF NONE OF THE KEY 
VARIABLES ARE 
CHANGING WRITE 


|***NO CHANGE IN STATeI 
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TG0EM: This routine Computes the time-to-go with the next 

event for each criterion being monitored during the current phase. 
The smallest value is then selected and returned to CYCXM. 


r 


l 



I 


i 



SET TG0 .TO INFINITY 


© 







TG0EMI: This routine loads array IEVNT with the addresses of 

the events to be monitored during the current phase and specifies 
the order in which they are to be monitored. 
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TRACKER: This routine computes the slant range, angle of 

elevation, azimuth angle, look angles, and space losses between 
specified tracker stations and the vehicle* 
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TRAJ: This routine propagates a trajectory from the 

beginning of the specified phase to the final cutoff condition 
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TRAJX: This routine is called when the single trajectory 

option is requested (SRCHM = 0) . This routine calls overlays 
(2, 1), (2, 2), and (2, 3) in sequence to implement the trajectory 
simulation desired. 
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TRIM: This routine calculates the engine deflection angle 

required to balance the aerodynamic moments when using rocket 
engines, or the flap deflection angle required to balance the 
aerodynamic moments when using jet engines. 







TRYIT1 : This routine minimizes the estimated net cost func- 

tion as a function of the step-size parameter. The principal 
function of the routine is to setup of the data required by 
GEMNIN, where GENMIN is the routine that actually minimizes the 
function. 



COMPUTE CUI 
BOUNDARIES 
S TP MAX, THI 
P1NET, AND 
FIT TOLERAI 

IVE-FIT 
PCTCC AND 
i SLOPE OF 
THE CURVE - 
JCE EPSA(I) 



INITIALIZE 
STEP VARIAI 
COMPUTE THI 
TRIAL STEP, 

THE TRIAL 
iLES AND 
FIRST 
STEP (2) 


L_ 




CHECK TO SEE IF 
PINET IS MINIMIZED 
OF IF THE NUMBER OF 
ACTIVE CONSTRAINTS 
HAS CHANGED 



IF THE MINIMIZATION WAS 
SUCCESSFUL, SAVE ALL 
TRIAL STEP VARIABLES 


UPDATE THE INDEPENDENT 
VARIABLES 


I 
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TRYIT2 : 
function P2. 


This function minimizes the constraint error 



COMPUTE CURVE-FIT 
BOUNDARIES AND 
CURVE-FIT TOLERANCES 


GENMIN 



MINIMIZES P2 




CHECK TC 
P2 IS MI 

SEE IF 
NIMIZED 



^P2 X 

P2MIN^ — 


J RETURN 1 

xy 


CONSTRAINTS HAVE 
NOT BEEN SATISFIED, 
SO SET THE CON- 
VERGENCE FAILURE 
FLAGS IHADIT AND 
INTRBL 





TSPXM: This routine is the main program of overlay (2, 0) 
and controls the overall operation of the trajectory simulation 
routines . 
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TW0BDY : This routine solves the orbital equations of motion 

about a spherical planet using the method of Laplace. 



~ x 

COMPUTE THE SCALAR COEFFICIENTS FOrJ 
LAPLACE'S EQUATIONS AS A FUNCTION 1 
OF THE CHANGE IN ECCENTRIC ANOMALY | 

= I 

UPDATE THE POSITION 
AND VELOCITY VECTORS 





UNITDU: This routine unitizes the control correction 

vector and computes the magnitude of the control correction vector. 



CALCl 

DUMA( 

JLATE 

n 

J 



r ^ 

DUMAG : 0 



| NFLAG * 0 j 






COMPUTE A UNIT 
VECTOR IN THE 
DIRECTION OF 
SEARCH DU(I) 
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UPDATS : This routine updates the sensitivity matrix by 

deleting those rows that can be dropped. The routine also com- 
putes the projection matrix PR0J(I) and the error map matrix in 
subroutine PGM to determine the direction of search. 



! DELETE FROM AC0B(I) 

THE ROWS CORRESPONDING 
TO THE INACTIVE INEQUALITY 
CONSTRAINTS. THIS DETER- 
MINES THE CURRENT 
SENSITIVITY MATRIX SMAT(I) 



62 


CALCULATE PR0J(I) , 
EMAP(I) , AND PG1(I) 


CALCULATE THE MAGNITUDE 
PG1(I) AS PG1MAG 


| — 
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UPN0M: This routine updates the nominal values of the 

independent variables according the direction of search DU(I) 
and the step size GAMAS, 
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I 

! 

i 

i 

i 


WGTINI : This routine initializes the vehicle's weight for the 

current phase and adds the specified delta velocity if requested 
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WINDS: This routine calculates the components of the wind 

velocity vector in the Earth-centered inertial coordinate system. 



■compute 

— 


COMPUTE 

GEOGRAPHIC 


GEOGRAPHIC 

WIND VECTOR 


WIND VECTOR 

FROM AZIMUTH 


FROM COMPONENT 

|AND SPEEI 

) TABLES | 


TABLES 


TRANSFORM WIND 
VECTOR FROM 
GEOGRAPHIC TO 
INERTIAL COORDINATES 
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WUCAL: This routine calculates the weighting matrix for the 

independent variables. 
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XITER (A, B, C, SIP, CIP, TTMXB) : This routine solves the 

quadratic equation to determine the sine and cosine of the, thrust 
incidence angle (SIP and CIP) and the thrust moment required to 
balance the moments (TTMXB) , given the aerodynamic and thrust 
moment to be balanced (A) and the components of thrust (B and C) . 
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XRNGE1: This routine calculates the downrange and crossrange 

distances, based on relative great circles. 



COMPUTE DELTA 


COMPUTE DELTA 

LONGITUDE FOR 


LONGITUDE FOR 

RELATIVE RANGE 


INERTIAL RANGE 



1 


COMPUTE DELTA 
ANGLES FOR 
CROSSRANGE & 
DOWNRANGE 


COMPUTE AVERAGE RADIUS 
BETWEEN INITIAL AND 
PRESENT LOCATIONS 


COMPUTE CROSSRANGE 
AND DOWNRANGE 
DISTANCES 
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XRNGE2: This routine calculates the crossrange distance 

based on a ground track of a reference circular orbit. 



i 


COMPUTE THE INERTIAL 
LONGITUDE FROM THE 
ASCENDING NODE 

* 

COMPUTE AN0DER 
AND 0MEG0 

COMPUTE THE 
DOWNRANGE 
PARAMETERS 

l 

CALL S0LVE TO 
ITERATIVELY 
FIND THE CROSS- 
RANGE ANGLE 

j 

COMPUTE THE 
CROSS RANGE 
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IV. SERVICE ROUTINES 


The following service routines are used to perform arithmetic 
tasks that are used frequently throughout the program. 

ABT (A,B , C ,L ,M,N) : This routine computes matrix C given 

matrices A and B, the dimensions of A (L by M) , and the dimen- 
sions of B (M by N) as matrix A times the transpose of matrix B. 

ADDREL: This routine prints a namelist like statement for 

a single variable name. This output is for use in ODINEX executed 
programs . 

ANGLE2 (X,Y, ALPHA): This routine computes the angle ALPHA, 

in degrees, between vectors X and Y. ALPHA is measured counter- 
clockwise from X to Y and can range from 0 to 360 deg. 

ATANH (A): This function computes the angle, from 0 to 

2 pi radians, whose arc- tangent is A over B. 

BTW (B ,W,D,L ,M.N) : This routine computes matrix D as the 

transpose of matrix B times matrix W, where L, M, and N are the 
dimensions of matrices B and W. 

BUCKET (B,W,D,L,M,N) : This routine rearranges an array X 

and the corresponding elements of Y in ascending order. N Is the 
number of elements in each array. NP is a pointer indicating the 
first element of Y that Is less than the next element. NP is 
zero if it does not exist. The ordered arrays are returned as 
XX and YY. 

C0NICT: This routine calculates the Keplerian conic for the 

target vehicle. Same flow chart as C0NIC. 

C0SH (A) : This function returns as the hyperbolic cosine 

of A. 

CUBMIN (A,XMIN,YMIN) : This routine evaluates a cubic poly- 

nominal coefficients (A) , and returns the minimum values of the 
function (YMIN) and the corresponding value of the argument 
(XMIN) . 

CYCYM: This routine performs the cycling functions for the 

variable step/order predictor-corrector integrator. 

DIFTAB (TABLE, DELT) : This function computes the derivative 

of a monovariant table (TABLE) using finite differences over delta 
X (DELT). 
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DIGDIF (M,N,NDIF) : This routine finds the number of differ- 
ent digits (NDIF) between M and N, where M and N are base-10 num- 
bers. 


DISPLY : This routine displays certain parameters from the 

search/optimization algorithm on the cathode-ray tube of the 6400 
computer console. 

DP0LY3 (A,N) : This function evaluates the derivative of 

polynomial: A(l) + A(2)*X + A(3) x X**2 +A(4)*X**3 where N is 

:he address of X in the dictionary. 

DTAB (XS,YS,N,D): This function performs the table lookup 

tasks for double-valued (bivariant) tables. XS is the value of 
the X argument, YS is the value of the Y argument, N is the 
pointer array, and D is the table array. 

ERR0R (I,J): This routine writes an error message and de- 

termines if it is fatal or nonfatal. I and J are Hollerith words. 
If the first letter in I or J is blank, the error is fatal; other- 
wise it is nonfatal. 

EXPN (IBKT, IX, IL, IXBKT) : This routine moves the contents 

of IBKT from IXBKT down IL locations, and shifts the contents of 
IBKT from IX down IL locations. IXBKT is the last occupied cell 
in IBKT. 

FITER (FPDATA,IDATA,XN,FN,N,FCT) : This routine is called by 

subroutine S0LVE to find the zero of the function by either the 
Regula Falsi or secant methods. The arguments FPDATA and IDATA 
are the same as those defined for subroutine S0LVE. XN is the 
current value of X and FN is the corresponding function value. 

N is the current iteration number and FCT is the external func- 
tion defining the function. 

F0PMIN (X,Y,XMIN,YMIN, IERR) : This routine calculates the 

minimum of a polynomial based on four ordered pairs (X,Y) . The 
abscissa value that minimizes this cubic polynomial is returned 
as XMIN and the corresponding ordinate value is returned as YMIN, 

F0RMN: This routine computes the nutation in right ascen- 

sion. 

GENTAB (TABLE) : This function is the executive table look- 

up routine. It determines whether the table is constant-valued, 
monovariant, bivariant, or trivariant, and calls functions TAB, 
DTAB, and TRITAB as required. 

GSAI: This subroutine sets a flag (ISSEPS=0) to indicate 

no steering sensitivity matrix has been calculated. It is called 
once a phase. 
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INVM (A,N): This routine inverts an N by N matrix A and 

returns A as the inverse. 

IRMTRX (R,V,ETR): This subroutine computes the transforma- 

tion matrix (ETR) from EOT coordinates to target-centered relative 
coordinates where R = ECI position and V = ECI velocity. 

LEVEL: This routine calculates the angle of attack and 
throttling parameter required to maintain level flight. 

MADD (A,SA,B,SB,C,K,L) : This routine adds matris A and 

matrix B to form matrix C, where SA and SB are scalars that are 
applied to A and B before the addition is performed. The dimen- 
sions of A and B are given as K by L. 

MATPY (A,B,C,L,M,N) : This routine multiplies matrix A by 

matrix B to produce matrix C, where A is an L by M matrix and B 
is a M by N matrix. Thus, C is an L by N matrix. 

MTRXM (A,B,C): This routine multiplies matrix A by matrix 

B to produce matrix C, where the maximum size of A or B is 3 by 3. 

MTRXT (A,B,C): This routine multiplies matrix A by the trans- 

pose of matrix B to form matrix C, where A, B, and C are 3 by 3 
matrices . 

MTRXTV (AT,V,W): This routine multiplies the transpose of 

matrix A by vector V to form vector W, where A is a 3 by 3 matrix 
and V and W are 3 by 1 vectors. 

MTRXV (A,B,C): This routine multiplies a 3 by 3 matrix A by 

a vector B, and returns the answer as C. 

0RBTRT: This routine transforms the target orbit parameters 

to earth centered inertial coordinates. The flow chart is the 
same as 0RBTR. 

PAD (A,B,M0D): This routine determines the delta X that 

produces the most precise derivative without losing significance 
in being rounded off. 


Argument 

Calls 

Call 


1 to (n - 1) 

n 

A 

f (X) 

AX 

B 

f(X + AX) 

Dummy 

HOD 

0 

1 


P0LY (N1,C0FI,ARGI) : This function evaluates a polynomial 

of degree N1 with coefficients C0FI as a function of ARGI. 
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P0LY1 (A,X): This function evaluates a polynomial of degree 

1 with coefficients A as a function of X. 

P0LY2 (A,X): This function evaluates a polynomial of degree 

2 with coefficients A as a function of X. 

P0LY3 (A,X): This function evaluates a polynomial of degree 

3 with coefficients A as a function of X. 

PPT (P,C,M,N,S): This routine computes matrix C as matrix 

P times the transpose of matrix P, where P is an M by N matrix 
and S is a scalar that must be set to zero in the calling program. 

RRTBR: This subroutine calculates body rates, given rela- 

tive Euler rates . 

REVOAT: This routine contains date of last program revision. 

SERCH (K,IV) : This routine locates the address of variable 

I (Hollerith) with respect to the beginning of the common refer- 
ence . 


SETESN: This routine resets the valves of all variables in 

thedictionary to their values at the minus side of event ESNI. 

SINH ( A ) : This function returns as the hyperbolic sine of 

A. 


S0LVE (FPDATA, IDATA, XSTAR, FCT) : This routine is a gen- 

eralized one-dimensional iteration scheme that finds the zero of 
the function FCT. The argments required are: 


FPDATA(l) - 
FPDATA(2) « 
FPDATA (3) - 
FPDATA(4) - 
FPDATA(5) = 
FPDATA (6) = 


IDATA (2) - 

IDATA (3) = 


IDATA(4) = 


IDATA (5) - 


initial guess on the value of the argument X 

the minimum value of X 

the maximum value of X 

the tolerance on the function value 

the tolerance on the value of X 

the increment of X in the secant method 

1 (Regula Falsi) 

2 (secant method) 

the maximum number of iterations allowed 
iteration flag 

1 (iteration limit was reached) 

2 (X within tolerance) 

3 (the iteration converged) 
type of solution to accept 

0 (ignore slope of the function) 

1 (positive slope of function only) 
debugging print selector 

0 (no debugging printout) 

1 (debugging printout) 
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XSTAR ■ the value of X that satisfies the iteration 

FCT - the external function whose zero is to be found. 

SP (X,Y,N): This function computes the scalar product of 

two N-dimensional vectors X and Y. 

SEEL (LIST) : This routine locates the addresses of the vari- 

ables in the array LIST, where LIST is constructed as follows: 

NAME1,NAME2,X 1 , 

NAME3 ,NAME4 ,X 2 , 

NAME5,NAME6,X 3 , 

The routine is also used to find the addresses of the derivatives 
and their corresponding integrals. 

SVDQ : This routine performs the variable step/order pre- 

dictor-corrector integrator. 

SVDQI : This routine performs the initialization for SVDQ. 

TAB (XS,N,X,Y): This function is a generalized table look- 

up scheme for single-valued (monovariant) tables. The table 
lookup is performed using either linear or cubic interpolation. 

XS is x* (the value of x for which the table value is desired) . 

N is the number of pairs of points in the table. X and Y are the 
values of x and f (x) , respectively. If the value of x* lies be- 
yond the table values of x, extrapolation will occur. If linear 
interpolation is desired, the x values immediately above and below 
x* are used. If cubic interpolation is desired the two x values 
immediately below x* and the two values immediately above x* are 
used. The array N in the argument list is defined as: 

N(l) = number of pairs of points in the table 

N(2) = type of interpolation desired 
* 1 (linear interpolation) 

■ 3 (cubic interpolation) 

N(3) = type of x's 

= -1 (decreasing) 

= 1 (increasing) 

N(4) = pointer to the value of x that is just below the 
last value of x* 

TRPM (X,Y ,XMIN ,YHIN) : This routine fits a quadratic poly- 

nomial through three points. It returns the minimum of this 
polynomial as YMIN and the minimizing value of X as XMIN . 
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THP0SM (X,Y,DYDX1,XMIN,YMIN): This routine fits a cubic 

polynomial using three points and the slope of the function. 

The routine returns the minimum of the polynomial as YMIN and 
the minimizing value of X as XMIN. 

TP0SM (X,Y ,DYDX1 , XMIN , YMIN) : This routine fits a quadratic 

polynomial using two points and the slope of the function. It 
returns the minimum of the function as YMIN and the minimizing 
value of X as XMIN, 

TRITAB (XS ,YS ,ZS ,N ,D) : This routine is a trivariant table 

interpolator that finds F(X*,Y*,Z*). 


D 

es 

the array containing the table to be inter- 
polated . 

N (1) 

= 

number of x’s 

N(2) 

= 

number of y's 

N(3) 

= 

number of z's 

N(4) 

= 

type of interpolation 

N(5) 

= 

type of x* s 

N(6) 

= 

type of y's 

N(7) 

= 

type of z T s 

N(8) thru (20) 

= 

pointers for the last used point on the 
table. 

X*, Y* , Z* 


values of the table arguments 


UNIT (X,XBAR,XMAG) : This routine computes the magnitude, 

XMAG, and unit components XBAR, of a 3-dimens lonal vector X. 

UPN0MS: This subroutine stores the current value of the steer- 

ing controls (US) into the dictionary. 

VCR0SS (A,B,C): This routine computes the vector cross- 

product, C, between two three-dimensional vectors A and B. 

VD0T (X,Y) : This function computes the dot product between 

two three-dimensional vectors X and Y. 

ZER0 (A,M,N) ; This routine zeroes the M by N vector A. 

EPHEM (DATE, GULIE, GHA, RAS , DECL) : This routine computes 

the Greenwich hour angle, and the right ascension and declination 
of the sun given the month /day /year date and the floating point 
Julian date from 1950.0 . 

FORMIN (CENJUL, GULIE1, RANUT) : This routine computes the 

nutation of the right ascension of the Sun given the true Julian 
Century and the Julian date. 
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SUN (AJP, ARA, DEC): This routine computes the right ascen- 

sion and declination of the Sun. 

RUK8: This routine performs eight-order Runge-Kutta inte- 

gration . 
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V. DEFINITION OF INTERNAL FORTRAN SYMBOLS 


This section presents the symbols used internally in the 
program. Most variables in the program are located in COMMON to 
conserve core locations. Certain variables are local variables to 
a specific routine. These types of variables are not shown in 
this list but are presented along with the flow chart for that 
routine in the flow chart section. Variables that are either in- 
put or output variables are defined in Volume II. 

Mathematical symbols are presented for each variable where 
applicable. The common that contains the variable is also shown, 
along with a definition of the variable and the routine that de- 
fines it. 


FORTRAN MATH 

SYMEOL SYMBOL COMMON 

A ( I ) A LOCAL 

1 = 1,2 

ABU) - DPGVC 

1 = 1.9 

ACOB(I) - SEARC 

1 = 1,625 

ACOBSU) - LOCAL 

1 = 1,16 

AR( I ) MOTVC 

1 = 1,15 

AXLU) - MOTVC 

1 = 1,3 

B (I ) LOCAL 

1 = 1,3 

CCI) C LOCAL 

1 = 1,2 

CADPS - LOCAL 

CADYS - LOCAL 

CODPS - LOCAL 

CDOYS - LOCAL 

CLOPS - LOCAL 

CHOPS - LOCAL 

CHOPS - LOCAL 

V-2 


DEFINITION/SUBROUTINE 


TOTAL AERO AND THRUST MOMENTS IN 
PITCH ANO YAW TO BE CANCELLED BY 
THE TRIMMING ENGINES/ TRIM 

THE TRANSFORMATION MATRIX FROM 
THE ATMOSPHERIC RELATIVE SYSTEM 
TO THE CODY SYSTEM/ IBMTRX 

THE JACOBIAN OF THE CONSTRAINT 
VECTOR/ GRAD 

THE JACOBIAN OF THE STEERING 
CONSTRAINT VECTOR/ GRADS 

THE TABLE LOOK— UP VALUE OF EXIT 
AREA FOR ENGINE 1/ PROP 

THE INERTIAL ACCELERATION VECTOR 
DUE TO GRAVITY AS PROPAGATED 
BY THE LAPLACE TWO-BODY 
EQUATIONS/ MOTENC 

THE ELEMENTS OF THE B VECTOR USED 
TO CLMPUTE THE QUATERNION RATES/ 
GUID2 

ThRUST TIMES THE PITCH AND YAW 
MOMENT ARMS, RESPECTIVELY/ TRIM 

THE TABLE LOCK-UP OF CADPT/ AERO 

THE TABLE LOOK-UP OF CADYT/ AERO 

THE TABLE LOOK-UP OF CDDPT/ AERO 

THE TABLE LOOK-UP OF CDDYT/ AERO 

THE TABLE LOOK-UP OF CLDPT/ AERO 

THE TABLE LOOK-UP OF CMDPT/ AERO 

THE TABLE LOOK-UP OF CNDPT/ AERO 



DEFINITION OF 


FORTRAN MATH 

SYMBOL SYMBOL COMMON 



CNPl(I) 

— 

DYTEM 

: ::r~ 

1=1,30 



1 - 

CTHAT 

- 

SEARC 


CWDYS 


LOCAL 


CYCF 

- 

CYC VC 

i 

CYDYS 

- 

LOCAL 


D ( I ) 


LOCAL 

- ■ — 

1 = 1,3 



I : 

DAX(I) 


MOTVC 


1=1,3 



I . -s : 

DELT 

- 

CYCVC 

r - ' 

DG( I ) 



SEARC 

_ 

1=1,25 



| :§L 

DFVLK(I) 

_ 

MOTVC 

, Si , 

1 = 1,3 




DMI(I) 

mm 

MOTVC 


1=1,15 



> 

DPR 

- 

SERVC 


DTIME 

- 

CYCVC 

[ w 

DTO 

— 

CYCVC 



FORTRAN SYMBOLS ( CONTD) 


DEFINITION/SUBRCUTINE 


CORRECTED SOLUTION IN THE FOURTH 
ORDER PREDICTOR-CORRECTOR 
INTEGRATION FORMULA/ FOPC 

COSINE OF CONEPS ( 1 )/ TEST 

THE TABLE LCOK-UP OF CWDYT/ AERO 

CYCLING FLAG SUCH THAT IF =0, THIS 
IS A DERIVATIVE PASS WITH DELT=0/ 
CYCXMI, CYCXM 


THE TABLE LOOK-UP OF CYDYT/ AERO 

THE ELEMENTS OF THE D VECTOR USED 
TO COMPUTE THE QUATERNION RATES/ 
GUID2 

THE DIFFERENCE BETWEEN THE INERTIAL 
ACCELERATION COMPUTED VIA MOTION 
AND THE INERTIAL ACCELERATION 
COMPUTED VIA TWOBDY/ MOTENC 

CURRENT INTEGRATION STEP SIZE/ 

CYCXM 

THE DIFFERENCE BETWEEN THE COST 
GRADIENT ON TWO SUCCESSIVE ITER- 
ATIONS/ DGM 

THE HISTORY VALUES OF THE FUNCTIONAL 
INEQUALITY CONSTRAINT DERIVATIVES/ 
MOTION 

3fT / 

THE TABLE LOOK-UP VALUE OF FLOWRATE 
FOR ENGINE 1/ PROP 

tffSREES PER RADIAN CONVERSION 
FACTOR/ BLKDAT 

THE DERIVATIVE OF TIME, I.E., 1.0 

THE INTEGRATION STEP SIZE IN THE 
CURRENT PHASE/ CYCXMI 
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DEFINITION OF INTERNAL FORTRAN SYMBOLS (CONTD) 


FORTRAN 

SYMBOL 

MATH 

SYMBOL 

COMMON 

DEFINITION/SUBROUTINE 

DTSAV 

- 

DYTEM 

THE LAST INTEGRATION STEP SIZE. 
USED TO CHECK FOR CHANGES IN STEP 
SIZE WHEN USING THE PREDICTOR- 
CORRECTOR FORMULA/ DYNS2 

DXPII) 
DYP(I) 
DIP (I) 
1=1*15 


LOCAL 

THE THRUST MOMENT ARM OF ENGINE 1/ 
TRIM 

DXR {11 
1=1*3 


LOCAL 

THE AERODYNAMIC MOMENT ARM/ TRIM 

DYNIL(l) 

1=1,148 


DYNIL 

AN ARRAY WHICH CONTAINS THE 
INTEGRALS, DERIVATIVES, AND A 
FLAG FOR EACH VARIABLE TO INDICATE 
WHETHER IT IS TO BE INTEGRATED/ 
BLKDAT 

EACI) 

1=1,25 

- 

SEARC 

THE ERROR VECTOR FOR THE ACTIVE 
CONSTRAINTS/ PGM 

EMAPtl) 

1=1,625 


SEARC 

THE TRANSFORMATION THAT DETER- 
MINES THE DIRECTION OF SEARCH FOR 
CONSTRAINT SATISFACTION/ UPDATS 

ENOIS 

- 

CYC VC 

A SMALL NUMBER USED AS A TOLERANCE 
TEST 

ENOM(I) 

1=1,25 


SEARC 

THE NOMINAL TARGET ERRORS/ NOMINL 

ENOMS(I) 

1=1,4 

- 

DPGVC 

THE NOMINAL STEERING ERRORS/ GRADS 

EPSAII) 

1=1,4 


SEARC 

THE STEPSIZE CONTROL FOR THE ONE- 
DIMENSIONAL MINIMIZATION ROUTINE* 
I.E., UPPER AND LOWER BOUNDS AND 
CURVEFIT ERROR TOLERANCES/ FGAMA, 
TRYIT1 , TRYIT2 

ESN 

— 

TGOVC 

THE CURRENT EVENT SEQUENCE NUMBER/ 
DINPT 

ESNPR1 

- 

INFVC 

THE CURRENT EVENT SEQUENCE NUMBER 
FOR PRINTOUT/ INFXMI 



DEFINITION OF INTERNAL FORTRAN SYMBOLS (CONTD) 


0 

FORTRAN 

MATH 




SYMBOL 

SYMBOL 

COMMON 

DEFINITION/SUBROUTINE 

ILrir 

EVNTtI ) 


REDAT 

THE EVENT SEOUENCE NUMBER AND 

= : — - 

1=1,8 



CRITERIA ARRAY/ READAT 

===== 

EVTF 


PHZVC 

EVENT FLAG/ PHZXMI 

: ~ 




=0, NOT AN EVENT 

1 




=1, ON MINUS SIDE OF AN EVENT 

i 




=2, ON PLUS SIDE OF AN EVENT 

, 

EXTRAP 


INF VC 

THE VALUE OF TIME WhEN THE LAST 

1 




PROFIL TIME SLICE WAS WRITTEN. 

i 




USED TO COMPUTE THE PROFIL WRITE 

! 




INTERVAL/ INFXM 


FPP5 

— 

SERVC 

FLOATING POINT NUMBER = .5 

- 

FPOO 

• 

- 

SERVC 

FLOATING POINT NUMBERS 0,1,. ..,15 


• 

FP15 





FP60 

- 

SERVC 

FLOATING POINT 60 

T 

FP90 

— 

SERVC 

FLOATING POINT 90 


FP160 

- 

SERVC 

FLOATING POINT 180 

L 

FP270 

- 

SERVC 

FLOATING POINT 270 

1 • 

FP360 

- 

SERVC 

FLOATING POINT 360 


FUXN ( I ) 


TGOVC 

THE HISTORY VALUES OF THE CRI- 

■ • rip.-; 

1=1,10 



TERIA VARIABLES CURRENTLY BEING 





MONITORED/ TGOEM 


GAMASS 


SEARC 

THE STEPSIZE FOR THE NON-UNITIZED 

V- 




DIRECTION OF SEARCH/ TRYIT1 

r 

GAMAX 

, 

SEARC 

THE MAXIMUM STEPSIZE ALLOWED 

V 




CONSIDERING THE INACTIVE 





INEQUALITY CONSTRAINTS, AND 

— 



' ' “ ' 

Tfir Maximum allowed percentage 

I 




CHANGE IN THE CONTROL 

I - 




PARAMETERS/ TRYIT1 

r : - ^ 
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DFFINITION 

CF INTERNAL 

FORTRAN SYMBOLS (CONTD) 


FORTRAN 

SYMBOL 

MATH 

SYMBOL 

COMMON 

DEFINITION/SUBRCUTINE 

— 

GP(I) 

1-1*9 

- 

DPG VC 

THE TRANSFORMATION MATRIX FROM 
THE GEOGRAPHIC SYSTEM TO THE EODY 
SYSTEM/ IBMTRX 

w 

GCON(l) 
1 = 1*3 


LOCAL 

CONSTANTS USED BY GRAV IN 
COMPUTING THE GRAVITATIONAL 
POTENTIAL/ MOTIAL 


GPCI ) 
1=1,25 

— 

SEARC 

THE COST GRADIENT ON THE PREVIOUS 
ITERATION/ PGM, DGM, DGMP2 


GSQOLD 


SEARC 

THE SQUARE OF THE GRADIENT MAGNI- 
TUDE ON THE PREVIOUS ITERATION/ 
CGM 


HARGC 1 ) 
1 = 1,3 


DPGVC 

hi STORY VALUE OF THE ARGUMENT 
OF ThE GDF(I) TABLE USED IN 
GENERALIZED LINEAR FEEDBACK 
STEERING/ GUIDX 


HDG( 1 ) 
1=1^25 


SEARC 

THE ESTIMATE OF THE HESSIAN 
MATRIX GENERATED IN THE DAVIDON 
ALGORITHM/ DGM 


HEADER (I) 
1=1,10 

— 

INF 1C 

THE TITLE TC BE PRINTED OUT AT THE 
TOP OF EACH PAGE/ INFXMI 


HEATC 

— 

LOCAL 

CONSTANT USED IN AERO HEATING 
CALCULATION/ MOTIAL 


HERRCR(I) 
1=1 ,3 


DPGVC 

HISTORY VALUE OF ERROR TERM IN 
GENERALIZED LINEAR FEEDBACK 
STEERING EQUATIONS/ GUIDX 


HESS (1 ) 
1=1,325 


SEARC 

THE APPROXIMATION TO THE HESSIAN 
MATRIX GENERATED IN THE DAVIDON 
ALGORITHM/ DGM, DGMP2 , PGM 

\ 

I ,J,K,L,M 

- 

SERVC 

INTEGER VARIABLES FOR TEMPORARY USE 


I A ( I ) 

1 = 1,9 


DPGVC 

TRANSFORMATION MATRIX FROM PLANET 
CENTERED INERTIAL TO ATMOSPHERIC 
RELATIVE FRAME/ IBMTRX 
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DEFINITION 

OF INTERNAL 

FORTRAN SYMBOLS ( CONTD ) 

FORTRAN 

SYMBOL 

MATH 

SYMBOL 

COMMON 

DEFINITION/SUBROUTINE 

IACS(I) 

1=1,25 

- 

SEARC 

THE SAVED VALUES OF THE SUBSCRIPTS 
OF THE ACTIVE CONSTRAINTS/ TRYIT1 

IBKTU ) 
1=1,2000 


£LANK 

A DATA BUFFER WHICH CONTAINS THE 
EVENT CRITERIA AND THE TABLE INPUT 
DATA/ READAT, DINPT 

I FTC ( I ) 
1=1,25 


0VRLY25 

i 

A COMBINATION OF ACTIVE CONSTRAINT 
INDICES USED TO DETERMINE IF SOME 
CONSTRAINTS CAN BE MADE INACTIVE/ 
DROP 

ICASE 

- 

MULTRC 

CURRENT CASE ( PROBLEM ) NUMBER/ 
READAT 

I CD 


OVRLY25 

INDEX OF THE ACTIVE CONSTRAINT 
WHICH WAS DROPPED/ REVISE 

ICGM 


SEARC 

A FLAG TO INDICATE THAT INITIAL- 
IZATION MUST BE DONE IN THE 
CONJUGATE GRADIENT ROUTINE (CGM) 
IF ICGM IS NON-ZERO/ DE LTU 

IDAV 

- 

SEARC 

NOT USED. 

IDENT(l) 
1 = 1,9 

- 

SERVC 

THE IDENTITY MATRIX (3X3)/ BLKDAT 

IDTAB(I) 

1=1,5 


CYCVC 

ADDRESSES OF THE TABLES THAT ARE 
TO BE USED IN COMPUTING THE 
INTEGRATION STEP SIZE DURING 
THE CURRENT PHASE/ CYCXM, DTMDL 

IDX 

- 

REDAT 

INDEX USED IN READAT TO MERGE 
MULTIPLE RUN DATA/ READAT 

IENT 

— 

RED AT 

CURRENT ESN BEING MATCHED DURING 
GENERATION CF MULTIPLE RUN DATA/ 
READAT 

IERRX 


REDAT 

USED TO SAVE FATAL ERROR FLAG 
BEFORE EVALUATING NON-FATAL 
ERRORS/ READAT 


V-7 



DEFINITION OF INTERNAL FORTRAN SYMBOLS (CONTD) 


FORTRAN 

SYMBOL 

MATH 

SYMBOL 

COMMON 

DEFINITION/SUBROUTINE 

IESN 

- 

PHZVC 

THE INITIAL EVENT SEQUENCE NUMBER/ 
DINPT 

IEVNT(I> 

1*1.10 

- 

TGOVC 

THE ARRAY OF EVENT LOCATIONS 
CURRENTLY BEING MONITORED/ TGOEMI 

IFLG 

- 

CYC VC 

A FLAG USED TO INDICATE THAT AN 
EVENT HAS BEEN INITIATED/ CYCXM 

IFRST 

- 

REDAT 

A FLAG TO INDICATE THAT THIS IS 
THE BEGINNING OF A PROBLEM/ READAT 

IGCI ) 
1*1.9 


DPGVC 

THE TRANSFORMATION MATRIX FROM THE 
•PLANET CENTERED INERTIAL TO THE 
GEOGRAPHIC SYSTEM/ MOTION 

IG6N(I) 

I=1,20G0 

- 

GENRL 

A DATA BUFFER WHICH CONTAINS THE 
GENERAL INPUT DATA/ READAT, DINPT 

J HAD IT 


SEARC 

A FLAG WHICH IS SET TO 1 IF THE 
PROGRAM COULD NOT GET TARGETED ON 
THE LAST OPTIMIZATION STEP/ TRYIT1 
TRYIT2 

II 

- 

REDAT 

USED AS INDEX AND COUNTER DURING 
DATA PROCESSING IN READAT/ READAT 

IL( I ) 
1*1.9 

— 

DPGVC 

THE TRANSFORMATION MATRIX FROM THE 
PLANET CENTERED INERTIAL TO THE 
LAUNCH INERTIAL SYSTEM/ MOTIAL 

IMAX 

- 

SEARC 

THE LAST PHASE NUMBER OF DEPPH(I) 
AND OPTPH TO OCCUR IN TIME/ MINMYS 

IMIN 

- 

SEARC 

THE FIRST PHASE NUMBER OF DEPPH(I) 
AND OPTPH TO OCCUR IN TIME/ MINMYS 

IMLT 

- 

REDAT 

POINTS TO THE BEGINNING OF LABELED 
COMMON MNMMLT IN IV/ READAT 

IMULT 

- 

MULTRC 

SIZE OF GENERAL DATA RECORD SAVED 
FOR MULTIPLE RUNS/ READAT 
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DEFINITION OF INTERNAL FORTRAN SYMBOLS (CONTD) 


FORTRAN 

SYMBOL 

MATH 

SYMBOL 

COMMON 

DEFINITION/SUBROUTINE 

IN 

10 

- 

MULTRC 

CURRENT INPUT AND OUPUT FILE FOR 
MULTIPLE RUN - FLIP FLOPS BETWEEN 
3 AND 4/ RE AD AT 

IND(I) 

1=1,25 

— 

0VRLY25 

COMBINATION OF CONSTRAINTS 
SELECTED BASED ON COMBINATORIAL 
PROCEDURES/ COMBINE 

INFF 

- 

INFVC 

A PRINT FLAG WHICH FORCES A PRINT 
WHEN SET NONZERO/ PHZXM 

1NIT 

— 

DYTEM 

A COUNTER IN PREDICTOR-CORRECTOR 
TO START THE ALGORITHM/ DYNXB, 
DYNS2 

INSRT 


RED AT 

FLAG USED TO INDICATE IF A TABLE 
HAS BEEN INSERTED FOR THE CURRENT 
PHASE BEING PROCESSED IN MULTIPLE 
RUN PORTION OF READAT/ READAT 

INTRBL 


SEARC 

A FLAG TO INDICATE THAT THE PROGRAM 
COULD NOT GET TARGETED ON THE 
CURRENT OPTIMIZATION STEP/ TRYIT2 , 
TEST 

INTRY1 

— 

SEARC 

A FLAG TO INDICATE THAT OVERLAY 
(2,5) IS TO BE CALLED WHEN SET 
NON-ZERO/ TRYIT1 

IOPT 

— 

SEARC 

A FLAG WHICH IS SET NON-ZERO TO 
INDICATE THAT THERE WERE NO PREV- 
IOUS OPTIMIZATION STEPS/ TRYIT1 

IPRT 

— 

GENIC 

A FLAG WHICH SUPPRESSES THE TRAJ- 
ECTORY PRINTOUT WHEN SET TO ZERO/ 
MINMYS 

1PRNTB 

- 

INFVC 

THE NUMBER OF FULL PRINT LINES IN 
THE PRINT BLOCK/ INFXM 

IPRN7R 

- 

INFVC 

THE NUMBER OF REMAINING PPINT 
VARIABLES IN THE LAST LINE THAT IS 
NOT FULL/ INFXM 
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DEFINITION OF INTERNAL FORTRAN SYMBOLS CCONTD) 
FORTRAN MATH 

SYMBOL SYMBOL COMMON DEFINITION/SUBROUTINE 


IRANGE 

IREVDT 

IRUNF 

IR1 ( 1 ) 

ISCORECI) - 
1 = 1,2 

ISFLG 

ISTART 

I stem 

1=1,25 
I STEPS 

ISTOP(I) - 
1 = 1 ,4 


IS ZB IK 
ISZEV 


AUXVC A FLAG TO INDICATE THAT THE RANGE 

CALCULATION HAS BEEN INITIALIZED/ 
AUXFMI 

GENIC REVISION DATE PRINTED AT THE 

BEGINNING OF EACH RUN IN READAT/ 
BLKDAT 

REDAT FLAG USED TO INDICATE IF CURRENT 
CASE IS TO BE BYPASSED. IRUNF=0 
IMPLIES CASE WILL BE BYPASSED/ 
READAT 

SERVC A TEMPORARY REUSABLE ARRAY 

MULTRC WORDS IN WHICH SUCCESS OR FAILURE 
OF RUNS ARE PACKED/ RSCORE 

SEARC FLAG WHICH INDICATES THAT THIS 

IS A RESTART ITERATION FOR THE 
DA VI DON OPTION/ R SEARCH 

SEARC A FLAG TO INDICATE THAT THE 

DAVIDON ALGORITHM IS TO BE RE- 
STARTED IF SET NON-ZERO/ TRYIT1 , 
PGM 

OVRLY25 INDICES OF THE SAVED TIGHT 
CONSTRAINTS/ DROP 

DPGVC A FLAG WHICH IS SET NON-ZERO IF A 

STEERING SENSITIVITY MATRIX HAS 
BEEN COMPUTED IN THIS PHASE/ GSA 

SEARC AN ARRAY TO INDICATE HOW THE 

CURRENT PROBLEM TERMINATED/ 

NOMINL , TEST 

ISTOPI 1 )=77B, PROBLEM SOLVED 
ISTOPC 2 )=77B, ITERATION LIMIT 
IST0P(3)=77B, NO CHANGE IN STATE 
IST0P(4)=77B, TIME LIMIT 

MULTRC SIZE OF SEARCH DATA RECORD SAVED 
FOR MULTIPLE RUNS/ READAT 

TGOVC THE NUMBER OF EVENTS BEING 
MONITORED/ TGOEMI 
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DEFINITION 

FORTRAN MATH 
SYMBOL SYMBOL 

ISIZ 

ISV 

ITC(I) 

1 = 1 1 25 

ITERF 

ITRIM 

IVCRT(I) - 
1=1,3 

IVETA 

IVSZ 

IWTFLG 

IXBKT 


OF INTERNAL FORTRAN SYMBOLS (CONTD) 


COMMON DEFINITION/SUBROUTINE 


REDAT THE SIZE OF THE TABLE CURRENTLY 

BEING STORED IN IBKT/ READAT 

MOTVC A POINTER TO THE LAST POINT USED 

IN THE ATMOSPHERE TABLE LOOK-UP/ 
ATMOS2 

SEARC ThE INDICES OF THE TIGHT INEQUAL- 

ITY CONSTRAINTS/ NOMINL, REVISE, 
DROP 

SEARC A FLAG TO INDICATE THE TYPE OF 

ITERATION STEP/ MINMYS 
=0, TARGETING ONLY 
■It OPTIMIZATION ONLY 
=2, TARGETING AND OPTIMIZATION 

LOCAL FLAG TO INDICATE IF TRIM IS BY 

ENGINES (= 0> OR BY FLAPS (= 1)/ 
MOTIAL 


DPGVC ARRAY OF INDICES POINTING TO THE 

EVENT CRITERIA IN IV USED TO 
REACH THE DESIRED ANGLE WHEN 
USING THE GENERALIZED LINEAR 
COMMANDS/ GUIDI, GUID1 

DPGVC INDEX POINTING TO EVENT CRITERIA 

USED TO REACH DESIRED ETA WHEN 
USING GENERALIZED LINEAR COMMANDS/ 
MOTIAL, PROP 

SERVC THE SIZE OF THE COMPUTATIONAL 

COMMON REGION (END-IV+1)/ READAT 


SEARC A FLAG TO INDICATE THAT THE 

DEPENDENT AND INDEPENDENT 
VARIABLES ARE WEIGHTED OR 
UNWEIGHTED/ WUCAL 

REDAT THE NUMBER OF CELLS OCCUPIED BY 

THE CONTENTS OF IBKT/ READAT 
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DEFINITION OF INTERNAL FORTRAN SYMBOLS (CONTD) 


FORTRAN 

SYMBOL 

MATH 

SYMBOL 

COMMON 

DEFINITION/SUBROUTINE 

IXE 

- 

REDAT 

A POINTER TO THE EVENT IN IBKT/ 
PE AD AT 

IXEVG 

- 

REDAT 

A POINTER USED IN THE GENERAL 
DATA ARRAY ( IGEN)/ READAT 

IXEVT 

- 

REDAT 

A POINTER USED IN THE TABLE DATA 
ARRAY (IBKT)/ READAT 

IXEVT 

- 

TGOVC 

AN INDEX ON IEVNT/ TGOEMI 

IXG 


REDAT 

POINTS TO BEGINNING OF CURRENT 
PHASE IN NEW CASE DATA SET/ 
READAT 

IXGEN 

* 

REDAT 

THE NUMBER OF CELLS OCCUPIED BY 
THE CONTENTS OF IGEN/ READAT 

I XT 


REDAT 

AN INDEX USED TO PACK THE TABLES 
INTO IBKT/ READAT 

IXI 

1X2 

— 

REDAT 

INDICES USED TO BUILD THE DATA 
BUFFERS/ READAT 

IXTBL 

— 

REDAT 

THE CORE LOCATION OF THE TABLE 
BEING INSERTED INTO IBKT/ READAT 

jacu) 

1 = 1 ,26 

— 

SEARC 

THE INDICES OF THE ACTIVE 
CONSTRAINTS/ TRYIT1 

JMLT 

— 

REDAT 

POINTS TO END OF LABELED COMMON 
HOLINC IN IV/ READAT 

JMULT 

— 

MULTRC 

SIZE OF TABLE DATA RECORD SAVED 
FOR MULTIPLE RUN/ READAT 

KREEP 


SEARC 

A FLAG WHICH INDICATES IF THE 
ITERATION IS PROGRESSING TOWARDS 
A SOLUTION/ TEST 

LB( I) 
1=1,9 

— 

DPGVC 

THE TRANSFORMATION MATRIX FROM 
THE LAUNCH INERTIAL TO THE BODY 
SYSTEM/ IBMTRX 
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DEFINITION OF INTERNAL FORTRAN SYMBOLS (CONTD) 


FORTRAN 

SYMBOL 

MATH 

SYMBOL 

COMMON 

DEFINITICN/SUBROUTINE 

NFLAG 

- 

SEARC 

A FLAG WHICH IS SET POSITIVE TO 
INDICATE THAT ThE PROBLEM HAS 
CONVERGED/ TEST 

NINDVX 

- 

REDAT 

NUMBER OF INDEPENDENT PARAMETERS 
SAVED FROM LAST RUN/ READAT 

NLDADR(I) 

1=1,25 

- 

REDAT 

ADDRESS OF THE CONTROLS SAVED 
FROM THE LAST SEARCH RUN/ READAT 

NLDPH(I) 

1=1,25 

— 

REDAT 

PHASE NUMBERS ASSOCIATED WITH 
THE CONTROLS SAVED FROM THE LAST 
SEARCH RUN/ READAT 

NOMF 

— 

SEARC 

A FLAG WhICH IS SET NON-ZERO IF 
THE TRAJECTORY BEING RUN IS A 
NOMINAL TRAJECTORY/ NOMINL 

NOO 

- 

SERVC 

FIXED POINT ZERO 

NOl 

m 

- 

SERVC 

FIXED POINT NUMBERS 1,...,15 

m 

m 

N015 




NPAGE 

- 

MULTRC 

PAGE COUNTER/ PAGER 

NPASS 

- 

DYNVC 

AN INTEGRATION PASS FLAG/ DYNXA, 
RUK 

NPC9 

NPC13 

NPC17 

— 

LOCAL 

INPUT VALUES OF NPCC9) , NPC(13), 
AND NPC<17I THAT ARE REQUIRED FOR 
PRINTOUT/ PRINT IC 

NSTEP 

- 

SEARC 

: ■ :i St . 

AN ITERATION COUNTER FOR THE 
CURRENT PROBLEM/ MINMYS 

NTC 

- 

SEARC 

THE NUMBER OF TIGHT CONSTRAINTS/ 
--REVISE 

1* ’ \r :,.2 ■■ - 
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DEFINITION OF INTERNAL FORTRAN SYMBOLS (CONTD) 


FORTRAN 

SYMBOL 

MATH 

SYMBOL 

COMMON 

DEFINITION/SUBROUTINE 

NULL 

- 

SERVC 

A VARIABLE USED TO DETECT IF INPUT 
VARIABLES HAVE BEEN INPUT. THE 
STORED VALUE OF NULL IS 1HU/ 

BLKDAT 

OLDG2 


SfcARC 

G2MAG FROM THE PREVIOUS ITERATION/ 
TEST 

OLDP1 

OLDPZ 


SEARC 

P1NOM AND P2NOM FROM THE PREVIOUS 
ITERATION/ TEST 

OLDU 


SEARC 

UMAG FROM THE PREVIOUS ITERATION/ 
TEST 

OMGSLT 

- 

LOCAL 

OMEGA*SIN(LATC)/ MOTIAL , MOTION 

ONE 

— 

MNMMLT 

THE NAME OF A CELL WHICH CONTAINS 
FLOATING POINT ONE/ NOMHOL 

PCTGC 


TGOVC 

PERCENT OF NOMINAL INTEGRATION 
STEP ALLOWED BY TIME-TO-GO LOGIC 
TO BRACKET THE DESIRED FUNCTION 
VALUE/ TGOEM 

PCTOLD 

' 

SEARC 

THE MAXIMUM PERCENTAGE CHANGE 
ALLOWED ON THE PREVIOUS OPTIM- 
IZATION STEP/ TRYIT1 

PE ( I ) 
1=1,198 


INF IC 

THE CORE ADDRESSES AND NAMES OF THE 
CURRENT PRINT VARIABLES/ INFXMI, 
INFXM 

PGT(I) 
1=1, 45 


TRACKC 

PLANET TO TRACKER GEOGRAPHIC 
TRANSFORMATION MATRIX. (3X3 
MATRIX FOR 5 TRACKERS)/ AUXFMI, 
TRACKER 

PH2F 


PHZVC 

A FLAG TO INDICATE THAT THE TRAJ- 
ECTORY IS TO BE TERMINATED/ PHZXM 
=0, DO NOT TERMINATE THE TRAJECTORY 
NE 0, TERMINATE THE TRAJECTORY 

PI 


SERVC 

PI/ BLKDAT 
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1 

i 

■ 

FORTRAN 

SYMBOL 

MATH 

SYMBOL 

COMMON 

DEF1NITICN/SUBR0UTINE 


LIMGEK 

- 

REDAT 

THE MAXIMUM SIZE OF THE GENERAL 
DATA ARRAY ( IGEN) / READAT 

| 

LIMBKT 

- 

REDAT 

THE MAXIMUM SIZE OF THE TABLE 
DATA ARRAY (IBKT)/ READAT 


LlSTD(l) 
1=1 ,30 

- 

DYTEM 

THE ADDRESSES OF THE CURRENT 
DERIVATIVES/ DYS11, DLGOK 


LISTI(I) 

1=1,30 

- 

DYTEM 

THE ADDRESSES OF THE CURRENT 
INTEGRALS/ DYS11, DLGOK 

- 

LMBKT 

- 

REDAT 

THE CURRENT SIZE CF IEKT/ READAT 

:: 

LMEVT 

- 

REDAT 

THE SIZE OF THE CURRENT EVENT 
BEING INPUT/ READAT 

5 — 

LNGTH 

- 

REDAT 

THE LENGTH OF THE COMPUTATIONAL 
COMMON AREA/ READAT 

H"'. 


LPRNT 

- 

INF VC 

THE LAST PRINT TIME/ INFXMI, 
INFXM 


N 

- 

DYTEM 

THE NUMBER OF INTEGRALS IN THE 
CURRENT PHASE/ DYSI1, DLOOK 


NAGS 

- 

SEARC 

THE NUMBER LF ACTIVE CONSTRAINTS 
INCLUDING THE EC'UALI TY CON- 
STRAINTS/ TRYIT1 

| 1: 

NAMSVR(I) 

1=1,51 

— 

SEARC 

i« r . ■ 

THE HOLLERITH NAMES OF THE 
CONTROL, TARGET, AND OPTIM- 
IZATION VARIABLES/ READAT 

i ; 


NDISC 

— 

DYNVC 

Jt . .S 

* := '% 

A FLAG TO SIGNAL THE INTEGRA- 
TION ALGORITHM THAT THERE IS A 
DISCONTINUITY/ DYNXA 


NEQC 

- 

SEARC " ? 

The number of equality con- 
straints/ MINMYS 

= - 

NETF 

— 

SEARC 

A FLAG WHICH IS SET NON-ZERO IF 
SRCHM=4 AND OPT IS NON-ZERO/ 
MINMYS 
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FORTRAN 

MATH 



SYMBOL 

SYMBOL 

COMMON 

DEFINITION/SUBROUTINE 


PIF - PHZVC A PROGRAM INITIALIZATION FLAG/ 

PHZXMI 

=0, THIS IS THE INITIALIZATION PASS 
*1, THIS IS NOT THE INITIALIZATION 
PASS 


PIC? 


SERVC 

PI OVER TWO/ BLKDAT 

PITMOM 

YAWMCM 

- 

LOCAL 

THE TOTAL MOMENTS IN PITCH AND YAW/ 
TRIM 

PNMCMI) 

1=1,30 


DYTEM 

FIRST DIFFERENCE BETWEEN PREDICTED 
SOLUTION AND CORRECTED SOLUTION 
USED IN FOURTH ORDER PREDICTOR- 
CORRECTOR FORMULA/ FOPC 

PNPl(l) 

1=1,30 

m m 

DYTEM 

PREDICTED SOLUTION IN FOURTH ORDER 
PREDICTOR-CORRECTOR INTEGRATION 
FORMULA/ FOPC 

PROJ(I) 

1 = 1 ,625 

- 

0VRLY25 

THE MATRIX WHICH DETERMINES THE 
PROJECTED GRADIENT/ UPDATS, REVISE 

P1N0M 

P2N0M 

— 

SEARC 

THE VALUES OF PI AND P2 ON THE 
NOMINAL TRAJECTORY/ MINMYS 

P2MIN 

- 

SEARC 

THE LOWER BOUND ON THE WEIGHTED 
ERROR MAGNITUDE/ DELTU 

OLDU ( I ) 
1=1,25 

- 

RED AT 

CONTROLS SAVED FROM LAST SEARCH 
RUN/ READAT 

OSREF 

- 

LOCAL 

DYNAMIC PRESSURE TIMES THE 
AERODYNAMIC REFERENCE AREA/ MOTION 

RERP2 

- 

LOCAL 

(RE/RP )**Z/ MOTIAL 

REVNT 

- 

REDAT 

A VARIABLE USED TO SET A FLAG FOR 
A ROVING EVENT/ READAT 

RPD 

- 

SERVC 

RADIANS PER DEGREE CONVERSION 
FACTOR/ BLKDAT 
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FORTRAN 

SYMBOL 

MATH 

SYMBOL 

COMMON 

DEFINITION/SUBROUTINE 

S(I) 

1=1,25 

- 

SEARC 

THE VALUE OF DU FROM THE PREVIOUS 
ITERATION/ CGM 

SALFA 

CALFA 

- 

LOCAL 

THE SINE AND COSINE OF THE ANGLE 
OF ATTACK/ GUIDI, MOTION 

SAVE ( I ) 
1=1,70 

- 

TGOVC 

THE VALUES OF THE STATE VARIABLES 
AT THE LAST INTEGRATION STEP/ TGOEM 

SAVI T ( I , J ) 
1=1,27 
J= 1 » 5 

— 

SEARC 

AN ARRAY IN WHICH THE RESULTS OF 
THE CURVEFIT OF PI AND P2 ARE 
SAVED/ FGAMA 

SE.ET A 
CBETA 

- 

LOCAL 

THE SINE AND COSINE OF THE ANGLE OF 
SIDESLIP/ GUIDI, MOTION 

SFCT(I) 

1=1,15 

— 

LOCAL 

THE TABLE LOOK-UP VALUE OF THE 
SPECIFIC FUEL CONSUMPTION FOR 
ENGINE 1/ PROP 

SGAM 

CGAM 

— 

LOCAL 

THE SINE AND COSINE OF THF PATH 
ANGLE RELATIVE TO THE ATMOSPHERE/ 
GUIDI, GUID2 

SIDEAL 

— 

AUXVC 

THE SAVED VALUE OF VIDEAL USED 
TO COMPUTE THE REQUIRED VELOCITY 
MARGIN/ AUXFMI 

SINDPH(I) 

1=1,25 

- 

SEARC 

THE VALUES OF INDPH(I) WHICH ARE 
SAVED IN READAT FOR PRNTIC/ READAT 

SIPT 

C1PT 

— 

LOCAL 

THE SINE AND COSINE OF THE THRUST 
INCIDENCE ANGLE IN PITCH REQUIRED 
TO TRIM THE VEHICLE/ TRIM 

SIYT 
Cl YT 

— 

LOCAL 

THE SINE AND COSINE OF THE THRUST 
INCIDENCE ANGLE IN YAW REQUIRED TO 
TRIM THE VEHICLE/ TRIM 

SLAM 

CLAM 

- 

LOCAL 

THE SINE AND COSINE OF THE AZIMUTH 
ANGLE RELATIVE TO THE ATMOSPHERE/ 


GUIDI , GUID2 
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FORTRAN 

SYMBOL 

MATH 

SYMBOL 

COMMON 

DEFINITION/SUBROUTINE - 

SLAMRF 
CLAMP F 

- 

LOCAL 

THE SINE AND COSINE OF THE REFERENCE 
AZIMUTH ANGLE USED IN THE RANGE 
CALCULATIONS/ AUXFMI 

SLATC 

CLATC 

TLATC 


LOCAL 

THE SINE, COSINE, AND TANGENT OF 
THE GEOCENTRIC LATITUDE OF THE 
VEHICLE/ MOTIAL, MOTION 

SLATRF 

CLATRF 


LOCAL 

THE SINE AND COSINE OF THE 
REFERENCE LATITUDE USED IN THE 
RANGE CALCULATIONS/ AUXFMI 

SLG 

CLG 

- 

LOCAL 

SINE AND COSINE OF GEODETIC 
LATITUDE/ MOTIAL 

SLONG 

CLCNG 

- 

LOCAL 

SINE AND COSINE OF INERTIAL 
LONGITUDE/ MOTIAL, MOTION 

SOLVED 

— 

MULTRC 

FLAG SET TO INDICATE HOW THE LAST 
PROBLEM TERMINATED/ SCORE 

SSIGMA 

CS1GMA 

— 

LOCAL 

THE SINE AND COSINE OF THE BANK 
ANGLE/ GUID1, MOTION 

SSTKI ) 

1 = 1 ,625 


OVRLY25 

A MATRIX USED IN COMPUTING THE 
PROJECTION MATRIX AND IS EQUAL 
TO <S*S PRIME )INVERSE/ UPDATE, 
PGM 

S'SVIDL 


AOXVC 

THE ACCUMULATED SUM SQUARED VALUE 
OF VIDEAL FGR EACH STAGE. 

COMPUTED AT THE BEGINNING OF EACH 
PHASE/ AUXFMI 

STEMP( 1) 
1= 1 » 25 

- 

SERVC 

AN ARRAY USED FOR TEMPORARY 
STORAGE 

STEP ( I ) 
1 = 1,6 


SEAPC 

THE TRIAL STEP LENGTH FOR EACH 
TRIAL STEP/ GENMIN 

STKLAT(I) 
CTKLAT ( 1 ) 
1 = 1,5 


TRACKC 

SINE AND COSINE OF GEOCENTRIC 
LATITUDE OF TRACKING STATION 1/ 
AUXFMI, TRACKER 
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FORTRAN 

SYMBOL 

MATH 

SYMBOL 

COMMON 

STKLON(I) 
CTKLON ( I ) 
1 = 1,5 

- 

TRACKC 

STMINP(I) 
1 = 1,2 


SEARC 


STPMAX SEARC 



TABLE ( 1 > 
1=1,1500 

TEMP ( I ) 

1 = 1,£>0 

TGO 

TGRAD 

THRSTT 

TIC I ) 
1=1,15 

TIMIN 

TIMX 

TJD(I) 

1=1,15 


REDAT 

SERVO 
1 GO VC 
SEARC 

LOCAL 
MOT VC 
SEARC 

TGOVC 

COCAL 


FORTRAN SYMBOLS (CONTD) 


DEFIN1TI0N/SUBRGUTINE 


SINE AND COSINE OF THE LONGITUDE 
OF TRACKING STATION 1/ AUXFMI 
TRACKER 

THE MINIMUM STEP SIZE DECREASE 
FFGM THE PREVIOUS STEP WHEN 
GENERATING THE CURVEFIT/ GENM1N 

THE LENGTH OF ThE STEP IN THE 
DIRECTION OF SEAFCH REQUIRED TO 
REACh ThE BOUNDARY FC'R THE NEAREST 
INEQUALITY CONSTRAINT/ TRYIT1, 

TRY IT2 

THE CURRENT TaE LE BEING INPUT. USED 
TO TRANSFER THE TABLE FROM INPUT TO 
THE STORAGE ARRAY I BKT/ RTAB, READAT 

AN ARRAY USED FOR TEMPORARY 
STORAGE 

THE TIME TO GO TO ThE NEXT EVENT/ 

TGO E M 

TIME REQUIRED TO COMPUTE ALL 
TRAJECTORY SENSITIVITIES/ 

MINMYS 

TOTAL THRUST OF THE TRIMMING 
ENGINES/ PROP 

ThE VALUE OF NET THRUST FOR 
ENGINE 1/ PROP 

THE TIME AT WHICH THE EARLIEST 
PHASE IN INDPH(l) OCCURS/ NOMINL , 
SAVIC 

THE LAST VALUE OF TIME/ TGOEM 

The table look-up value of thrust 

FOR ENGINE 1/ PROP 
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FORTRAN 

SYMBOL 

MATH 

SYMBOL 

COMMON 

DEFINITION/SUBROUTINE 

TKLATC(I) 
1 = 1,5 

- 

TRACKC 

GEOCENTRIC LATITUDE OF TRACKING 
STATION 1/ AUXFMI 

TKP/D ( I ) 
1 = ] ,5 

— 

TRACKC 

GEOCENTRIC RADIUS TO TRACKING 
STATION 1/ AUXFMI 

TREF 


CYC VC 

A TIME REFERENCE USED TO COMPUTE 
THE NEXT STEP SIZE (DELT)/ CYCXMI, 
CYC XM 

TREFP 

— 

DYNVC 

ThE TIME REFERENCED TO THE LAST 
PRIMARY EVENT, CYCXMI, TGOEMI 

TREFS 

— 

DYNVC 

THE TIME REFERENCED TO THE LAST 
SECONDARY EVENT/ CYCXMI, TGOEMI 

TRKXRKI) 

1=1.15 


TRACKC 

THE EARTH CENTERED INERTIAL 
COORDINATES OF TRACKING STATION 
I. ( X , Y , Z ) FOR 5 TRACKING 
STATIONS/ AUXFMI, TRACKER 

TTRIAL 


SEARC 

TIME REGUIRED FOR A TRIAL STEP 
IN THE ONE DIMENSIONAL 
MINIMIZATION ROUTINE/ FGAMA 

TWO PI 


S ER VC 

TWO TIMES PI/ ELKDAT 

TVACA 


LOCAL 

TOTAL VACUUM THRUST OF THE 
ENGINES USED TO LIMIT ACCELER- 
ATION/ PRCP 

VXL ( I ) 
1 = 1.3 


MOTVC 

THE INERTIAL VELOCITY VECTOR 
PROPAGATED BY THE LAPLACE EQUATIONS 
IN TWOBQDY/ MOTENC 

VXLO ( I ) 
1 = 1 ,3 


MOT VC 

INITIAL INERTIAL VELOCITY VECTOR 
TO BE PROPAGATED BY THE LAPLACE 
EQUATIONS/ MOTENC 

WGTO 

— * 

LOCAL 

THE CALCULATED VALUE OF INITIAL 
STAGE WEIGHT/ WGT1NI 

WJETTO 

- 

LOCAL 

THE SAVED VALUE OF JETTISON WEIGHT 
FOR PRINTOUT PURPOSES/ WGTINI 
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FORTRAN 

SYMBOL 

MATH 

SYMBOL 

COMMON 

DEFINIT 10N/SUBRCUTINE 

WPRC'PC 

- 

LOCAL 

INITIAL STAGE PROPELLANT WEIGHT/ 
WTG1NI, PROP 

WPUSDh 

- 

LOCAL 

PROPELLANT CCNSUMFD UP TO AND 
INCLUDING LAST PHASE/ WTGINI 

XB ( I ) 

1 = 1,2 

B 

LOCAL 

THRUST TIMES THE PITCH AND YAW 
MOMENT ARMS/ TRIM 

XINF 

- 

SEF.VC 

FLOATING POINT INFINITY (10.E10)/ 
BLKDAT 

XIS AV ( I ) 
1=1,3 


AUXVC 

THE SAVED VALUES OF THE VEHICLE 
POSITION VECTOR RELATIVE TO THE 
EARTH FOR USE IN COMPUTING DPRNG1 
AND DPRNG2/ AUXFMI 

XL ( I ) 
1=1,3 

— 

MCTVC 

THE INERTIAL POSITION VECTOR 
PROPAGATED BY THE LAPLACE 
EQUATIONS/ MOTENC 

XLO( I ) 
1=1,3 

— 

MOTVC 

INITIAL INERTIAL POSITION VECTOR 
TO BE PROPAGATED BY THE LAPLACE 
EQUATIONS/ MOTENC 

XYCMA ( I ) 
XYOME (I ) 
XYOMS ( 1 ) 
1 = 1,3 


DPGVC 

CONSTANTS USED TO COMPUTE COMMANDS 
WHEN USING GENERALIZED LINEAR 
ANGLF OF ATTACK, SIDESLIP, AND BAN 
ANGLE COMMANDS/ GUIDI, GUID1 
XYDMA(1)=1NIT1AL VALUE OF CRITR 
XYGMA(2 >=INITIAL ANGLE VALUE 
XYOMA (3 )=SLGPE 

XYOME(I) 
1 = 1,3 


DPGVC 

CONSTANTS USED TO COMPUTE THROTTLE 
COMMAND WHEN USING GENERALIZED 
LINEAR COMMANDS/ MOT1AL, PROP 
XYCME (1)=IN1TIAL VALUE OF CRITR 
XYOME (2)=INITIAL ETA VALUE 
XYOME (3)=SL0PE 
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FORTRAN 

SYMBOL 

MATH 

SYMBOL 

COMMON 

DEFINITIGN/SU6RC.UTINE 

XYOM1 ( I ) 
XYC M2 ( I ) 
XY0M3 ( I ) 
1=1 ,3 


DPGVC 

CONSTANTS USED TO COMPUTE COMMANDS 
WHEN USING GENERALIZED LINEAR 
INERTIAL OR RELATIVE EULER ANGLE 
COMMANDS/ GU1DI, GU1D1 
XY0M1(1)= INITIAL VALUE OF CRITR 
XY0M1 (2)=INITIAL ANGLE VALUE 
XY0M1 ( 3 )= SLOPE 

YAWRC 

- 

LOCAL 

VALUE OF YAWR SAVED AT END OF 
A PHASE FOR PRINTOUT/ GUIDI 

YAWRH 


MOTVC 

HISTORY VALUE OF YAW ATTITUDE. 
USED WHEN RE-INITIALIZING BODY 
ATTITUDES BETWEEN EVENTS/ GUIDI 

YES ( 1 ) 
1 = 1, A 


SEARC 

THE VALUES OF THE POLYNOMIALS 
GENERATED DURING THE ONE 
DIMENSIONAL MINIMIZATION 
AT THEIR PREDICTED MINIMUMS/ 
TRYIT1 

YP ( I » J ) 
1 = 1,4 
J= 1 ,30 


DYTEM 

BACK DERIVATIVES USED IN FOURTH 
ORDER PREDICTOR -CORRECTOR 
FORMULA/ FOPC 

YPVAVU) 

1=1,30 

*“ 

DYTEM 

INITIAL VALUES OF INTEGRALS USED 
IN PREDICTOR-CORRECTOR FORMULA/ 
FOPC 

Yim 

Y2 ( I ) 
1=1,30 


DYTFM 

WORKING STORAGE USED BY THE RUNGE- 
KUTTA INTEGRATION ALGORITHM/ RUK 

ZlTRY(l) 

1=1,6 

- 

SEARC 

SAVED VALUES OF PITRY(I)/ TRYIT1 

Z2TRY ( I ) 
1 = 1,6 

- 

SEARC 

SAVED CALUES OF P2TRY<I)/ TRYIT1 

ZMAG 

- 

SEARC 

SAVED VALUE OF UMAG/ TRYIT1 

ZMAX 

— 

SEARC 

SAVED VALUE OF STPMAX/ TRYIT1 


V-22 



DEFINITION OF INTERNAL FORTRAN SYMBOLS (CCJNTD) 


FORTRAN 

SYMBOL 

MATH 

SYMBOL 

COMMON 

ZP1DS 

- 

SEARC 

ZP2DS 

- 

SEARC 

ZSTEPt I) 
1 = 1,6 

- 

SEARC 

ZYESII ) 
1=1,5 

- 

SEARC 


DEFINITION/SUBROUTINE 
SAVED VALUE OF 0P1DS/ TRY1T1 

saved value of dp2ds/ tryiti 

SAVED VALUES OF THE TRIAL STEP 
SIZES/ TRYITI 

SAVED VALUES OF YES Cl)/ TRYITI 
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VI. POST SUBROUTINE INDEX 


SUBROUTINE 


OCTAL 

NAME 

OVERLAY 

SIZE 


ABT 

(0,0) 

74 

ADDREL 

(2,6) 

116 

AERO 

(2,3) 

252 

AERCHI 

(2,3) 

155 

AERrA 

(2,3) 

106 

ANGLE2 

(0,0) 

115 

ANMPT 

(2,3) 

263 

ARTBR 

(2,0) 

71 

ATANH 

(0,0) 

32 

ATAN3 

(0,0) 

27 

ATMOS 

(2,0) 

34 

ATMOS 1 

(2,0) 

43 

ATMOS2 

(2,0) 

275 

ATMCS3 

(2,0) 

507 

AUXFM 

(2,3) 

614 

AUXFMI 

(2,2) 

374 

BACKCI 

(2,0) 

66 

EACKOR 

(2,0) 

66 

BLKDAT 

(0,0) 

7 

BTW 

(0,0) 

105 

EUCKET 

(0,0) 

123 

CALE 

(2,0) 

120 

CALES 

(2,3) 

26 

CALSPE 

(2,3) 

21 

CGM 

(2,5) 

43 

CLGM 

(2,3) 

30 

CLSPFL 

(2,0) 

114 

COMB IN 

(2,5) 

56 

CONIC 

(2,3) 

544 

CONICT 

(2,3) 

547 

CON VO 

(2,3) 

226 

COSH 

(0,0) 

32 

CUBMIN 

(0,0) 

131 

CYCXM 

(2,3) 

220 

CYCXMI 

(2,2) 

116 

CYCYM 

(2,3) 

67 

DATA 

(2,0) 

7 

DELTU 

(2,5) 

45 

DERIV 

(2,3) 

21 

DERVI 

(2,2) 

15 
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POST SUBROUTINE INDEX (CONTD ) 


SUBROUTINE 

NAME 


GRAV 

GS A 

CSAI 

GSENSR 

GU1DI 

Gumx 

GU1D1 

GUID2 

HSWGT 

IPMTRX 

INFXM 

INFXMI 

INPUTX 

INTGPL 

INVM 

1RMTFX 

IRTBR 

ITERO 

LEVFL 

MADD 

MASTER 

MATPY 

MINMYS 

MOMENT 

MCN1TR 

MOTENC 

MCTIAL 

MOTION 

M7RXM 

MTRXT 

MTRXTV 

MTRXV 

NMLTER 

NOMHf’L 

NCM1NL 

OLGM 

OLGCM 

ORETR 

ORBTRT 

PAD 


OCTAL 

SIZE 


114 

225 

16 

53 

523 
212 
600 

524 
47 

525 


(2,3) 

302/333 

(2,2) 

177 

( 1,0) 

65t> 

(2,2) 

37 

(0,0) 

547 

(2,0) 

63 

(2,0) 

54 

(2,6) 

1,250/1 ,545 

(2,0) 

20 

(0,0) 

35 

(0,0) 

7,343/11,^31 

(0,0) 

77 

(2,0) 

125 

(2,3) 

167 

(2,3) 

46 

(2,3) 

152 

(2,2) 

1 ,400 

(2,3) 

554/725 

(0,0) 

66 

(0,0) 

66 

(0,0) 

55 

(0,0) 

54 

(1,0) 

22 

(1,0) 

117 

(2,0) 

220 

(2,0) 

10 

(2,0) 

15/35 

(2,2) 

241 

(2,2) 

241 

(2,0) 

166 


OVERLAY 


2,3) 

2,3) 

2 , 2 ) 

2,3) 

2 , 2 ) 

2 , 0 ) 

2 , 0 ) 

2,3) 

2,3) 

2 , 0 ) 



T 
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SUBROUTINE 


OCTAL 

NAME 

OVERLAY 

SIZE 


PAGER 

(0,0) 

55/104 

P BLOCK 

(2,3) 

1,043 

PGM 

(2,5) 

235 

PHZXM 

(2,3) 

172 

PHZXMI 

(2,2) 

45 

POLY 

(0,0) 

46 

PGLY1 

(2,0) 

22 

PCLY2 

( ? , C ) 

2o 

PGLY3 

(2,0 

33 

PPT 

(0,0) 

107 

PRNTIC 

(2,2) 

4 ,627/ *,,726 

PROP 

(2,3) 

614 

GMULT 

(2,2) 

74 

GRTATE 

(2,2 ) 

72 

QUAT1 

(2,2) 

60 

GUAT2 

(2,2) 

26 

QUAT3 

(2,2) 

52 

READAT 

(1,0) 

3,112 

PF VDA7 

(0,0) 

7 

REVISE 

(2,5) 

136 

R GENOA 

(1,0) 

773 

RRTRR 

(2,0) 

101 

RSCGRE 

( 1,0) 

100 

RSE ARC 

( 1,0) 

367 

RTAE 

(1,0) 

74 

RTELML 

(1,0) 

1 ,015 

RUK 

(2,3) 

147 

SAVIC 

(2,0) 

101 

SCORE 

(1,0) 

31 

SGM 

(2,5) 

20 

SEPCH 

(1,0) 

42 

SETESN 

(2,0) 

74 

SETIC 

(2,0) 

140 

SETIV 

(2,0) 

76 

SHRINK 

< 1,0) 

37 

SINH 

(0,0) 

32 

SOLVE 

(2,0) 

150 

SP 

(0,0) 

30 

SPEL 

(1,0) 

33 
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POST SUBROUTINE INDEX (CONTD ) 


SUBROUTINE 


OCTAL 

name 

OVERLAY 

SIZE 

SUN 

( 2,2) 

263 

SVDQ 

(2,3) 

2,067 

SVDQI 

(2,3) 

147 

TAB 

(2,0) 

510 

TEST 

(2,0) 

154 

TGOEM 

(2,3) 

376 

TGOEMI 

(2,2) 

331 

THPM 

(0,0) 

125 

THPOSM 

(0,0) 

160 

TMOTM 

(2,3) 

122 

TPOSM 

(0,0) 

64 

TRACKM 

(2,3) 

275 

TRAJ 

(2,0) 

202 

TRAJX 

(2,0) 

120 

TRIM 

(2,3) 

217 

TRITAB 

(2,0) 

350 

TRYIT1 

(2,0) 

323 

TRYIT2 

(2,0) 

160 

TSPXM 

(2,0) 

45 

TWOBDY 

(2,3) 

356 

UNIT 

(0,0) 

46 

UNITDU 

(2,0) 

40 

UPDATS 

(2,5) 

164 

UPNOM 

(2,0) 

22 

UPNOMS 

(2,3) 

25 

VCRCSS 

(0,0) 

41 

VDOT 

(0,0) 

30 

WGTINI 

(2,2) 

265 

WINDS 

(2,0) 

74 

WUCAL 

( 2,5) 

272 

X ITER 

(2,3) 

201 

XRNGE1 

( 2,3) 

114 

XRNGE2 

(2,3) 

217 

ZERO 

(0,0) 

27 
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